Экспериментальное определение сечений продуктов глубокого расщепления и деления в мишенях из natW и 181Ta при облучении протонами в диапазоне энергий 40 - 2600 МэВ тема автореферата и диссертации по физике, 01.04.16 ВАК РФ

Федеральное государственное унитарное предприятие «Государственный Научный Центр Российской Федерации Институт Теоретической и Экспериментальной Физики»

На правах рукописи

4803

У

Флоря Сергей Николаевич

Экспериментальное определение сечений продуктов глубокого расщепления и деления в мишенях из па'\¥ и шТа

при облучении протонами в диапазоне энергий 40 - 2600 МэВ

Специальность: 01.04.16 - Физика атомного ядра и элементарных частиц

АВТОРЕФЕРАТ

диссертации на соискание ученой степени кандидата физико-математических наук

1 5 СЕН 2011

Москва 2011 г.

4853163

УДК 621.039.584

Работа выполнена в ФГУП «ГНЦ РФ ИТЭФ», г. Москва.

Научный руководитель: доктор физ,- мат. наук Титаренко Ю.Е.

(ФГУП «ГНЦ РФ - ИТЭФ», г. Москва)

Официальные оппоненты: доктор физ.- мат. наук Голубев А.А.

(ФГУП «ГНЦ РФ - ИТЭФ», г. Москва)

доктор физ.- мат. наук Говердовский A.A. (ФГУП «ГНЦ РФ - Физико-энергетический институт им. А.И. Лейпунского», г. Обнинск)

Ведущая организация: ФГУП «Научно производственное

предприятие «Радиевый институт им. В.Г. Хлопина», (г. Санкт - Петербург)

Защита диссертации состоится 11 октября 2011 г. в 11 часов на заседании диссертационного совета Д 201.002.01 в конференц-зале ИТЭФ по адресу: г. Москва, ул. Б. Черемушкинская, д. 25.

С диссертацией можно ознакомиться в библиотеке ИТЭФ. Также диссертация и автореферат доступны по запросу через электронную почту Sergey.Phlorya@itep.ru.

Автореферат разослан « 15 » августа 2011 г.

Ученый секретарь диссертационного совета кандидат физ.-мат. наук

Васильев В.В.

Общая характеристика работы

По состоянию на конец 2010 года общий объем накопленного отработанного ядерного топлива (ОЯТ) на объектах Госкорпорации «Росатом» составляет ~ 21 ООО тонн. К этому количеству ежегодно добавляется ~ 650 тонн ОЯТ, выгружаемого из 32 гражданских реакторов, которые работают на территории России. Исходя из этого и учитывая 33 новых энергоблока, запланированных в генеральной схеме строительства АЭС, прогноз накопления ОЯТ (без его переработки) соответствует рис. 1. Видно, что количество ОЯТ будет увеличиваться на ~ 10 ООО тонн каждые 10 лет, при этом средняя суммарная активность 1 тонны ОЯТ составляет ~

Дата

Рис. 1. Прогноз наработки ОЯТ на реакторах России Таким образом, наряду с проблемами обеспечения безопасной эксплуатации атомных электростанций (АЭС) (техническими и террористическими) в качестве первоочередной стоит так же проблема обращения с ОЯТ.

В концепции обращения с ОЯТ, утвержденной в Госкорпорации «Росатом» в 2008 году, подтвержден базовый принцип обращения с ОЯТ - его

переработка, а в качестве стратегических направлений выдвинуты следующие:

1) окончательная изоляция (захоронение) образующихся при переработке радиоактивных отходов (РАО) (будет регламентироваться в федеральном законе о РАО);

2) создание системы долговременного контролируемого хранения ОЯТ (будет регламентироваться в федеральном законе о ОЯТ);

3)сбалансированное вовлечение продуктов регенерации ОЯТ в ядерный топливный цикл;

4) развитие технологий переработки ОЯТ (в России из выгружаемого ежегодно ОЯТ перерабатывается не более 15 %).

В то же время признается, что поскольку ОЯТ и РАО, образующиеся в процессе переработки ОЯТ, содержат ряд радионуклидов с большими периодами полураспада (табл. 1), невозможно полностью гарантировать их надежную изоляцию в хранилищах и захоронениях, обеспечивающую радиационную безопасность человека и окружающей среды на весь период их потенциальной опасности (по истечении которого удельная активность радионуклидов, содержащихся в ОЯТ и РАО, снизится до значений, позволяющих освободить их от регламентации норм радиационной безопасности).

Таблица 1. Радионуклиды, определяющие активность и токсичность ОЯТ

Временной интервал, лет Определяющие радионуклиды

до 100 55Ре, 58Со, 59№, 908г, 106Ки, Ш8Ь, 137',34Сз, 144Се, ,47Рш, ,55'154Еи

100- 1000 63-39М,60Со,,37С3, 1М8ш

1000- 10000 240-239Ри,241Ат

104 - 1 о5 14С, 59№, 937г, 94НЪ, 240'239Ри, 237Кр, 243Ат

>105 99Тс, п91, 239Ри

Два последних стратегических направления должны позволить исключить риски, связанные с длительным хранением.

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

материалов и перевода их в относительно короткоживущие или стабильные нуклиды.

Одной из перспективных трансмутационных концепций является электроядерная технология, которая и рассматривается в настоящей работе. Ее физической основой является использование подкритического энерговыделяющего бланкета с интенсивным внешним источником нейтронов. В качестве такого источника служит мишень из тяжелого материала, облучаемая протонным пучком сильноточного ускорителя (ток пучка ~ 10 мА или более, энергия ~ 1 ГэВ).

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

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

1) высокая общая наведенная активность МУ, что потребует регламентных ограничений по срокам эксплуатации и проведению ремонтно-профилактических работ;

2) образование в МУ радиоактивных продуктов реакций (РПР) с высокой радиотоксичностью (например, а-излучатели, изотопы йода, тритий), попадание которых во внешнюю среду создало бы значительную радиационную угрозу, которая может определить порядок и сроки утилизации

облученных МУ;

3) накопление в конструкции МУ вторичных нуклидов (как радиоактивных, так и стабильных), имеющих значительные величины сечений поглощения нейтронов и снижающих нейтронную эмиссионную эффективность («изотопное отравление» МУ);

4) образование газообразных элементов, которые, скапливаясь в пустотных дефектах твердых МУ, могут приводить к охрупчиванию и даже разрушению кристаллической решетки;

5) образование металлоидных элементов с высокой коррозионной способностью (фтор, хлор, кислород), вызывающих коррозию конструкционных элементов МУ.

Поэтому важнейшими исходными данными для количественных оценок значимости и последствий перечисленных выше эффектов, наряду с конструкционными и эксплуатационными характеристиками электроядерных установок (ЭлЯУ), являются сечения образования вторичных нуклидов при взаимодействии ускоренных протонов с МУ и конструкционными материалами ускорителя. Значительное число вторичных нуклидов и разнообразие возможных условий эксплуатации ЭлЯУ объективно вынуждают рассматривать в качестве основы для формирования константного обеспечения по протонным сечениям расчетные методы. Они могут быть основаны на использовании (альтернативном или совместном) существующих моделирующих компьютерных программ: Bertini, ISABEL, СЕМ03.02, INCL4.2+ABLA, INCL4.5+ABLA07, PHITS, CASCADE07 и т.д. Однако выявляются серьезные расхождения между расчетными и экспериментальными результатами, которые, в некоторых случаях, могут достигать нескольких порядков. Данное обстоятельство требует проводить непрерывную работу по их верификации и совершенствованию. Для этого необходимы опорные значения сечений РПР, полученные экспериментальным путем. Кроме того, для ряда РПР, имеющих первостепенную важность при оценке отмеченных выше последствий эксплуатации МУ и нуждающихся, в силу этого, в повышенной степени точности и достоверности, экспериментальные значения сечений образования будут являться, возможно, безальтернативной исходной информацией.

Экспериментальное определение сечений РПР при взаимодействии протонов промежуточных и высоких энергий с конструкционными материалами мишенного устройства ЭлЯУ и верификация высокоэнергетических транспортных программ, используемых для его расчета, являются предметами исследования настоящей диссертации, чем, принимая во внимание выше сказанное, определяется ее актуальность.

Цель диссертации состоит в формировании наборов экспериментальных значений сечений РПР, образующихся при взаимодействии протонов промежуточных и высоких энергий с тонкими мишенями из па1\¥ и шТа (перспективными материалами для использования в качестве мишенных узлов ЭлЯУ). Верификация высокоэнергетических транспортных программ, широко применяемых в мировой практике, с использованием полученных экспериментальных данных.

Исходя из этого, можно выделить основные задачи диссертации:

1) целевой анализ имеющейся экспериментальной информации по протонным сечениям дляпа'\У и |8|Та;

2) выбор диапазона энергий протонов для облучения экспериментальных образцов;

3) проведение облучений образцов, физических измерений и обработки их результатов;

4) определение сечений образования РПР для ""'ДУ и Та в широком диапазоне энергий протонов;

5) сравнение полученных экспериментальных значений сечений РПР, где это возможно, с опубликованными ранее результатами;

6) проведение расчетов по различным высокоэнергетическим транспортным программам для определения сечений РПР ви 181Та;

7) сравнение экспериментальных значений сечений РПР с соответствующими расчетными данными и качественная оценка отмеченных расхождений;

8) верификация высокоэнергетических транспортных программ, с

использованием полученных экспериментальных данных.

Методы исследований. Диссертационная работа использует

следующие основные методические приемы:

1) а-спектрометрия облученных протонами образцов исследуемых материалов без предварительной химической сепарации;

2) у-спектрометрия облученных протонами образцов исследуемых материалов без предварительной химической сепарации;

3) сравнение полученных результатов с имеющейся экспериментальной информацией и расчетными моделями.

Научная новизна выполненных исследований заключается в следующем: определено 1895 значений сечений, что позволило получить 192 функции возбуждения для natW и 173 для 181Та. Из них 99 и 76 функций возбуждения, соответственно, для natW и 181Та измерены впервые. В уже изученных функциях возбуждения для natW и 18|Та получено 723 и 497 значений сечений в других, ранее не изученных, диапазонах энергий. С использованием полученных данных оценена предсказательная способность и осуществлена верификация семи ВЭТ программ: MCNPX (Bertini, ISABEL), СЕМ03.02, INCL4.2, INCL4.5, PHITS, CASCADE07. На защиту выносятся:

1) результаты измерений сечений РПР в протонных реакциях для natW и 181Та при энергии протонов 40, 70, 100, 150, 250, 400, 600, 800, 1200, 1600, 2600 МэВ;

2) результаты сравнения полученных значений сечений с имеющейся экспериментальной информацией;

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

Практическое значение работы определяется возможностью использования ее результатов:

1) для верификации и совершенствования расчетных программ, используемых при создании и эксплуатации ЭлЯУ;

2) в качестве самостоятельных ядерных данных;

3) для дополнения мировых баз данных по протонным сечениям.

Структура и объем диссертации

Диссертация состоит из введения, пяти глав, заключения (общих выводов по диссертации) и приложения, оформленного отдельным томом; имеет объем 118 страниц печатного текста; содержит 46 рисунков, 29 таблиц и библиографию (159 названий).

Краткое содержание работы

Во Введении раскрыта актуальность темы диссертации, определенны цель, задачи и методология исследования. Обоснована научная новизна и практическая значимость диссертации, рассмотрены вопросы стратегии обращения с отработанным ядерным топливом и возможность создания замкнутого топливного цикла на базе ЭлЯУ.

В Главе 1 диссертации выполнено обзорное рассмотрение ЭлЯУ. Вольфрам и тантал рассматриваются как главные кандидаты на роль мишенного узла «твердого типа» в ЭлЯУ. Их преимущество, по отношению к мишенным узлам из свинца, висмута или свинцово-висмутовой эвтектики, заключается в отсутствии наработки изотопов полония, которые в основном и определяют радиотоксичность облученной мишени ЭлЯУ. Однако стоит отметить, что в мишенях из тантала и вольфрама не удаётся исключить полностью радиотоксичность, обусловленную а-эммиторами, из-за образования в них usGd, который образуется во всех тяжелых мишенных материалах при взаимодействии с ними протонов с энергией выше ~ 0.5ГэВ. Отметим, что к настоящему времени выполнена всего одна работа, в которой проведены систематические измерения образования I48Gd в nalW, 181Та и 197Au для протонов с энергией 0.6 и 0.8 ГэВ, что недостаточно для практических приложений, в которых исходная энергия протонов, как правило, не ниже ~ 1 ГэВ.

Приведены результаты целевого информационного анализа имеющихся баз ядерных данных по сечениям взаимодействия протонов с ядрами "atW и 181Та. Обращено внимание на существенный недостаток требуемых для создания ЭлЯУ ядерных данных при энергиях протонов выше 1 ГэВ, в связи с чем сделан вывод о необходимости значительного увеличения объема соответствующей экспериментальной информации.

В Главе 2 описана методика экспериментального определения

независимых скоростей реакций образования нуклидов Щ*1, когда продукт

реакции (7Ь А1) образуется в ходе ядерной реакции, и кумулятивных я,"" , когда продукт реакции А2) образуется в ходе всех процессов (как непосредственно в ходе реакции, так и распада независимо образовавшихся генетических предшественников). Приведены погрешности определения скоростей реакций, а также значения использованных сечений мониторной реакции 27А1(р,х)22На.

Метод а-спектрометрии применялся только для определения скорости реакции образования 148Сс1 (Т1/2=74.6 лет, Ец=3.183 МэВ, Т1=100%) без его химической сепарации. Метод у-спектрометрии применялся для определения скоростей реакций /? /?+ и 1Т нуклидов. В измерениях использовались: а-спектрометр на базе детектора АЗ 00-17-АМ с разрешением 17 кэВ на а-линии с энергией 5275.3 кэВ 243Аш и два у-спектрометра с НР детекторами ОС-2518, с разрешением 1.8 кэВ на у-линии с энергией 1332 кэВ 60Со.

Сечение образования любого РПР и его погрешность вычислялись по формулам:

где: ]?итШ - кумулятивная или независимая скорость реакции образования продукта реакции, Ф - плотность потока протонов.

Так как химическая сепарация |48С(1 из-за незначительного количества его образования в образце не проводилась, то определить напрямую скорость счета а-частиц в пике затруднительно. Поэтому методика её определения была построена на расчете удельной активности шСс1 в образце.

Удельная активность ,48Сс1 в образцах па'\У и 181Та определялась путем сравнения расчетных и экспериментальных а-спектров, полученных в идентичных геометриях измерения исходя из предположения, что для образцов одинаковой геометрической формы и размеров при равномерном распределении швс1 в исследуемом образце скорость счета в заданном

а

псгшгНп4 ,ситИпЫ ___

(1)

Ф

(2)

энергетическом диапазоне спектра пропорциональна удельной активности 148С(1. Пробеги а-частицы с энергией 3183 кэВ в па1\У и 181Та равны 10.09 и 9.88 мг/см2 соответственно. Расчетная скорость счета для образцов цилиндрической формы в выбранном энергетическом диапазоне а-спектра от их толщины, представленная на рис. 2, показывает, что их толщина много больше длины свободного пробега а-частиц в этих материалах.

При таком подходе необходимо учитывать, что помимо а-частиц детектор реагирует на р- и характеристическое у-излучение радиоактивных ядер, образовавшихся в образцах вместе с 148Сс1. Поэтому необходимо ввести граничную энергию Есиыг Яля выделения области регистрации а- «пика» свободной от р- и у- составляющей. Для определения граничной энергии Еси10й-проводились дополнительные эксперименты, в которых каждый экспериментальный образец экранировался алюминиевой фольгой толщиной 15 мкм. В качестве примера на рис. 3 приведен спектршТа (Ер=800 МэВ) с экраном и без него.

ш

^ 2 о ^

СО

о. о Ё «и

ш ч:

(0

I 04

5 о 5 10 15 20

толщина образца, мг/см2 Рис.2. Расчетная скорость счета а-детектора в зависимости от удельной толщины образцов па|'\Л' и шТа

т

т

-l-1-1-1—r~

1500 2000

энергия, кэВ

2500

' I 1 1

3000

Рис.3. Спектры |81Та (Ер= 0.8 ГэВ), измеренные с экраном (А1 фольга толщиной 15мкм) и без него

Для моделирования спектра а-частиц шОс1 использовалась программа МСИРХ, позволяющая рассчитать энергетическое распределение импульсов А(Е) в объеме а-детектора с учетом всех ядерно-физических процессов взаимодействия а-частиц со средой, учитываемых этой программой. С учетом выше изложенного, формула для определения скорости реакции образования

148,

ОсЗ имеет вид

j^cum _

1

cutoff

1

N,ag 4Mjr 1-е"

(3)

где Ат,п1} - скорость счета в расчетном спектре в диапазоне энергий от ЕсЫо$до 3183 кэВ, АшоЯ - скорость счета в измеренном спектре в том же диапазоне энергий, - время облучения образца, ыт - число ядер в образце.

Формулы для определения кумулятивных и независимых скоростей реакций при использовании метода у-снсктрометрии имели вид:

А 1

j^cumfmd

Ктаг-Ъ-е, Fy

ЯГ

МТае"Пг-е2УП

2 ^

(4)

(5)

Яр =

'41

/С"

13

Л Г, Л2]

ситНгиЛ №

К Тч-Пг

(6)

(7)

N таг'Ч г£ 7

где Л, = Л А\ = А\ , А\ = А гг-кИг - параметры, определяемые путем фитирования экспериментальных точек кривых распада «материнского» и «дочернего» ядер методом наименьших квадратов («материнский нуклид -индекс 1», «дочерний нуклид - индекс 2»); Л^ - число ядер в облучаемом образце; щ, г]2 - абсолютные выходы у- квантов; Х\, кг - постоянные распада; Б[, - абсолютная эффективность спектрометра при энергиях у- квантов Ех (первый нуклид) и Ег (второй нуклид); у^ - коэффициент ветвления, т.е. вероятность перехода материнского нуклида в дочерний; ^ - время облучения; Р1, - функционалы, учитывающие не стационарность облучения.

Поправки ки, учитывающие поглощение у-квантов в материале экспериментального образца, определялись с использованием полного сечения взаимодействия у-квантов у-ой энергии с веществом (барн/атом) и толщины экспериментального образца (г/см2).

В ряде случаев, когда условия измерения не позволяли определить коэффициент А\, определялись супракумулятивная скорость реакции и соответствующее сечение

]2

---~-(8)

' ^Ч у цтЛ __^ 2

л,-л2' 1 ИТч-п 1-е2-Р2

Плотность потока протонов и погрешность определялись из формул аналогичным (1) и (2), при условии, что сечение используемой мониторной реакции 27А1(р,х)22На известно.

Глава 3 посвящена описанию постановки эксперимента. Описывается процесс компоновки облучаемого стека из образцов и А1 мониторов. Приведены характеристики облучаемых образцов исследуемых материалов и их изотопный состав.

Облучения тонких мишеней из и шТа, в диапазоне энергий

(0.04 ГэВ - 2.6) ГэВ, выполнялись на быстром выводе синхротрона ФГУП «ГНЦ РФ ИТЭФ» У-10. Ускоренный до заданной энергии протонный пучок, состоящий из 4-х банчей длительностью ~ 250 не каждый, из ускоряющего кольца переводился в транспортный канал, обеспечивающий вывод протонов со следующими параметрами: интенсивность - ~ 2-Ю11 протон/сброс; поперечное сечение - эллипс с осями ~ 10x15 мм; полная длительность ~ 1 мке; частота следования импульсов ~ 15 имп./мин. Схема транспортного канала для облучения образцов вместе с элементами быстрого вывода протонного пучка представлена на рис. 4.

В связи с тем, что предложенная методика определения сечения 148Gd накладывает жесткие требования на определение энергии дискриминации Ecutofr, с особой тщательностью были проведены измерения энергетической калибровки а-спектрометра.

Калибровка а-спектрометра проводилась с помощью набора образцовых спектрометрических а-источников (ОСАИ) (239Pu, 238Pu, 226Ra и 233U). В качестве примера на рис. 5 приведен а-спектр 233U+238Pu+239Pu. После выполнения калибровочных измерений проводились измерения фона а-камеры, результаты представлены на рис. 6 (время набора составляло 105 часов). На рис. 7 представлена измеренная энергетическая калибровка, из которой видно, что положения центров пиков 226Ra смещены вниз примерно на 10 каналов относительно положений пиков остальных источников. По этой причине энергетическая калибровка проведена отдельно для 226Ra и остальных источников. Смещение линий 226Ra обусловлено различной толщиной их активных слоев.

пучка:

обозначения: 1 - подставка для размещения облучаемой мишени; 2 - токовый трансформатор; 3 - выходной фланец вакуумного протонопровода; 4 -поворотный магнит; 5 и 6 - дуплет квадрупольных линз, обеспечивающих фокусировку пучка; 7 - септум магнит; 8 и 9 - магнитные блоки ускоряющего кольца; 10 - кикер магнит с углом отклонения 15 мрад; 11 - корректирующий магнит (с углом отклонения 20-30 мрад)

На рис. 8 приведена зависимость эффективности регистрации спектрометра, полученная с использованием 239Ри (линия 5156.6 кэВ), как функции дистанции между источником и поверхностью детектора, измеренной без коллиматора. Измеренная эффективность регистрации хорошо согласуется с расчетной, полученной посредством программы для моделирования энергетического спектра в а-спектрометре ААБ!. Данная величина так же рассчитывалась с использованием программы МС№Х, где задавалась реальная геометрия образец- детектор, а-частицы задавались с начальной энергией 3183 кэВ, с изотропным угловым и однородным пространственным распределениями по всему объему каждого экспериментального образца. Значения эффективности регистрации для значения дистанции между источником и поверхностью детектора 4 мм, использованной в измерениях, представлены в табл. 3.

Таблица 3. Эффективности регистрации а-спектрометра

Дистанция от источника до детектора, мм Эффективность регистрации измеренная23'Ри Эффективность регистрации рассчитана по АА81 Эффективность регистрации рассчитанная по МС№Х£„

4 0.294±0.018 0.2942±0.0015 0.29485±0.00015

10 1 с; Ё го л га i 1 О : 0.01 233и гз9Ри 11 2'Ч228тч1 233и « жГи Ж \Ж ■ и1 Ш 24Ra 212в Thl 2,2& 1 1 226Th ¡2,гв1 21:25Т I 2,6Ро ia г- ^ t 1+rAf

1300 1400 1500 1600 1700 1800 1900 2000 канал

Рис. 5. Измеренный спектр ОСАИ 233U+238Pu+239Pu

О 200 400 600 800 1000 1200 1400 1600 1800 2000 канал

Рис. 6. Фоновый счет а-спектрометра

: о 2Ии | 1 П (ГП 1 Г| гм 1 1 1МЧ.Г 1 1 1 ... .................

- □ "Ъа

: V 238Ри

: Л 239Ри

Л энергия[кэВ]=А*канал+В

-А=3,434±0.005, В=21 ±7, х2=1.1 ■

- -----А=3.41210.006, В=86±10, х=0Л [

ЧЭШ ..............................I......... ...........................

1300 1400 1500 1600 1700 1800 1900 2000 канал

Рис. 7. Энергетическая калибровочная зависимость а-спектрометра

Рис. 8. Экспериментальная (239Ри) и расчетная зависимость (АА81, МС№Х) эффективности регистрации детектора от расстояния до образца. Погрешности измерений в пределах точек

Примеры измеренных а-спектров облученных образцов MW и 18'Та приведены на рис. 9. Как видно они существенно отличаются от спектров «тонких» стандартных а-источников на рис. 5. Их отличие обусловлено торможением а-частиц испущенных ,48Gd с энергией 3183 кэВ в материале экспериментального образца, что приводит к полному размытию монолинии 148Gd и регистрации детектором а-частиц от энергии 3183 кэВ до нижней энергии дискриминации спектрометра (~ 80 кэВ). Присутствие в облученных образцах большого количества ß-активных ядер приводит также к регистрации спектрометром их суммарного спектра. Различие их функций возбуждения, разная длительность облучения и время высвечивания обуславливают отсутствие постоянной энергетической границы ß- и а-спектров.

Как указывалось выше, в работе применялись два у-спектрометра, калибровка которых проводилась с помощью набора образцовых спектрометрических у-источников. Результаты и детали калибровки, в том числе и их абсолютные эффективности регистрации, приведены в диссертации. В качестве примера на рис. 10 приведен у-спектр natW.

Обработка измеренных у-спектров осуществлялась с помощью пакета программ GENIE-2000. Идентификация у-линии и расчет сечений образования ядер-продуктов производился с помощью разработанной в ИТЭФ программой SIGMA и базой ядерных данных PCNUDAT.

В Главе 4 приведены результаты измеренных экспериментальных сечений для исследуемых в настоящей работе natW и шТа. В итоге, методом у-спектрометрии было измерено 1013 для natW и 882 для ШТ сечений. Методом а-спектометрии было измерено по 6 сечений образования ,48Gd в natW и 181 Та. Количество полученных сечений образования РПР каждого типа по всем энергиям подробно представлено в табл. 4. Средняя точность определения сечений образования продуктов реакций составляет: natW - 14.8 %, 181Та -14.4 % (рис. 13 и 14). Полученные данные сравнивались с данными полученными в ZSR (Zentrum für Strahlenschutz und Radioökologie, Hannover) (табл. 5).

энергия, кэВ

Рис. 9. Пример измеренных а-спектров, облученных NatW и l8,Ta

Е, кэВ

Рис. 10. у -спектр natW (№5, Ер = 2.6 ГэВ, время набора 900 с)

s 150

i 100

50

Q>. знач. 11.151

Ср. знач. 14.7891

О 25 50 75 100 0 25 50 75 100

Погрешность, %

Рис. 11. Распределения погрешностей скоростей реакций и сечений natW

i 150

1100 я

ra о а

g 50

S

s У

Ср. знач. 10.26%

Ср. знач. |4Ж№

25

50

75 100 О 25 Погрешность, %

50

75

100

Рис. 12. Распределения погрешностей скоростей реакций и сечений шТа

Таблица 4. Общее количество измеренных выходов РПР различных типов по каждому эксперименту

Нуклид Энергия, ГэВ Тип выхода Всего

с+с* i i(m)

0.04 3 4 - 2 9

0.07 8 7 - 2 17

0.1 17 10 2 2 31

0.15 24 10 3 3 40

0.25 34 11 4 4 53

,8,Та 0.4 58 14 6 4 81

0.6 68 19 8 5 99

0.8 70 19 9 3 100

1.2 95 25 11 13 143

1.6 99 23 9 16 146

2.6 111 23 13 17 163

Нуклид Энергия, Тип выхода Всего

ГэВ с+с* 1 ¡(ш)

0.04 8 6 2 3 19

0.07 14 11 2 4 31

0.1 23 13 3 5 44

0.15 30 14 4 5 53

0.25 40 19 5 5 69

0.4 58 14 6 5 82

0.6 74 18 9 5 105

0.8 79 17 12 6 ИЗ

1.2 106 24 15 11 155

1.6 110 25 14 15 163

2.6 121 26 15 18 179

ВСЕГО 1250 353 152 153 1907

Таблица 5. Сравнение результатов ИТЭФ и гКС для шТа и па1\У

Энергия ИТЭФ, МэВ Энергия гэя, МэВ ИТЭФ N(/N/N3(7 гея К Нуклид Тип 1 Оетэф Рэо

'"Та

66 ± 1 70.8±2 17/6/3 8/6/3 2.31 Та ,75НГ с с с 224 ± 18 234 ± 19 2.28 ± 0.20 491 ±621 407 ± 18 13.7 ± 1.4 1.17

97 ± 1 97.3±2 118±2 31/13/11 18/13/11 1.89 ,й1и с с 7 ±12 1.48 ±0.16 52.8 ±7.3 3.25 ±0.43 1.25

148 ± I 148±1 156 ±1 40/18/16 23/18/16 2.25 |йЬи ,67УЪ ,иТт с с с 79.5 ±7.1 22 ± 2.6 2.91 ± 0.44 127 ±14 7.71 ±0.82 6.16 ±0.74 2.25

248 ± 1 256±1 53/22/21 29/22/21 1.57 171Ьи "ль ¡(т+§) ¡(т+я) 14.6 ±3.3 0.13 ±0.01 46 ±8 0.06 ±0.02 1.43

339±2 337±2 81/33/28 39/33/28 1.7 153ТЬ Н9Сс1 ,46а<1 '"Ей "5Еи с* с с с с 6.48 ± 0.92 2.83 ±0.29 1.29 ±0.13 1.78 ±0.23 0.55 ± 0.07 2.34 ± 0.4 1.14 ± 0.17 0.41 ±0.05 0.52 ±0.14 0.17 ±0.03 1.35

559 ±2 561±1 99/44/42 58/44/42 1.61 н6Еи шРт 139Се с с с 14.3 ±2 5.16 ± 0.71 2.16 ±0.26 1.34 ± 1.5 1.85 ±0.45 1.15 ± 0.15 1.31

799 ±2 102/47/47 61/48/47 1.18

1198 ±3 1180±1 143/55/54 63/55/54 1.24 8!у с 1.13 ± 0.1 0.57 ±0.14 1.22

1598 ±4 1570±1 150/56/55 68/56/55 1.29 !«у с 1.83 ±0.17 0.7 ±0.2 1.25

2605 ± 8 2580±1 165/73/70 86/73/69 1.39 |731.ц ,!|С(1 88у с с с 27.4 ±1.7 12.7± 1.3 3.87 ±0.36 49.2 ±6.2 26.1± 3.7 1.35 ±0.53 1.34

«•V

68 ± 1 68.6 ± 1.9 31/8/7 10/8/7 1.35 1»,н{. с 0.08 ±0.01 0.15 ±0.02 1.21

99 ± 1 99.6±1.12 95.5± 0.8 44/11/10 16/11/10 1.27 т1м с 18.9 ± 1.7 35.1 ±3.7 1.17

Энергия ИТЭФ, МэВ Энергия МэВ ИЛ ИТЭФ N</N/N,0 Р Нуклид Тип Фитэф Рз»

149 ± 1 147 ± 1.4 154 ± 1.3 53/15/14 19/15/14 1.18 гаЬи с 83.7 ±7.9 139.6 ± 14.6 1.П

249 ± 1 254 ± 1.7 260 ±1.1 69/22/22 27/22/22 1.37 1.37

400 ±2 334 ± 1.7 83/25/20 32/25/20 1.67 гаШ ,боЕг ШТЬ 149С<1 ,16оа с с с* с с 83 ± 8.4 19.3 ±2.1 2.43 ± 0.35 1.16 ± 0.13 0.48 ±0.05 134.3± 13.6 10.2 ± 1.7 0.63 ±0.14 0.51 ±0.06 0.13 ±0.02 1.36

600 ±2 559 ± 1.5 105/38/35 43/38/35 1.30 шЕи 96Тс "'■¿г с ¡(т+в) с 4.09 ±0.38 0.16 ±0.02 0.24 ± 0.03 2.44 ± 0.33 0.28 ± 0.04 0.49 ± 0.06 1.23

799 ±2 764 ± 1.5 114/40/37 49/40/37 1.40 "3Ьи 131Ва "Тс с с ¡(т+Й 55.4 ±5.6 1.47±0.14 0.19 ±0.03 100.2 ± 10.4 0.8 ± 0.1 0.39 ±0.06 1.26

1199 ± 3 1160±1.4 155/47/45 51/47/45 1.36 тЬи "й с с 43.1 ±4.3 0.26 ± 0.02 80.9 ± 8.2 0.55 ±0.06 1.31

1599 ±4 1570±1.4 163/49/47 55/49/47 1.29 ,73Ьи »8 у с с 36.4 ±3.7 1.8 ± 0.17 68.1 ±6.9 0.82 ±0.15 1.18

2605 ± 8 2590±1 179/55/50 56/55/50 1.4 ,ббуь и6са |4!Еи 8!у 56Со с с с с с 24.9 ±2.3 20.0 ± 1.9 18.1 ±1.8 4.27 ± 0.40 0.213 ±0.021 15.4 ± 1.8 8.56 ±0.87 11.8± 1.3 0.99 ±0.18 0.11 ±0.02 1.26

Результаты сравнения показывают, что число сечений РПР, измеренных в

ИТЭФ (Ы0) в ~ 4 раза превосходит число сечений, измеренных в гБЯ. Это обусловлено возможностью измерения в ИТЭФ короткоживущих продуктов реакций. Так как в составах измеренных продуктов реакций наблюдаются различия, то в качестве критерия согласованности экспериментальных результатов ИТЭФ и гБИ использовался среднеквадратичный фактор отклонения <Р>, который вычислялся по формуле (9). Фактор рассчитывался как для полной совокупности данных (<?>), так и для совокупности данных в которой значения, вышедшие за предел За, были отброшены

(Ц ~а1\чАсг1 <рЗа>). Выполненный анализ показал, что значение

фактора в наборах экспериментальных данных, полученных в ИТЭФ и 2811, и не выходящих за предел За составляет для шТа ~ 30 % и ~ 25 %, что намного ниже различия экспериментальных и расчетных результатов. Таким

образом, за счет более широкого и унифицированного диапазона энергий, в котором проводились измерения, и большего массива экспериментальных данных за счет измерения короткоживущих нуклидов данные ИТЭФ имеют преимущество по отношению к данным ZSR.

В Главе 5 полученные результаты были сравнены с расчетными моделями BERTIN, ISABEL, СЕМ03.0, INCL4.2+ABLA, INCL4.5+ABLA07, PHITS, CASCADE07. Вычисления проводилось при 25 значениях энергии протонов от 0.03 до 3.5 ГэВ. Предсказательная способность программ оценивалась с использованием среднеквадратичного фактора отклонения <F>:

5

где стсхр - сечение, полученное в данной работе, аса1 - расчетное сечение.

Поскольку программы вычисляют только независимые сечения продуктов реакций, а измеряются в основном коммулятивные (см, табл. 4), то для увеличения статистики сравнений, в соответствии с цепочками радиоактивных переходов, вычислялись соответствующие значения расчетных кумулятивных сечений. Необходимые формулы представлены в диссертации. В качестве примера на рис. 13 представлена схема образования |480с1. Результаты сравнения экспериментальных и расчетных сечений, полученные по формуле (9), приведены в табл. 6. Таблица 6. Значения фактора Г для па'\У и шТа

Модель Образец Энергия, МэВ

40 70 100 150 250 400 600 800 1200 1600 2600-

Bertini W 5.75 3.02 2.49 2.26 2.78 2.30 2.28 1.94 2.34 2.82 2.96

Та 7.43 2.07 1.83 1.54 2.40 2.40 2.29 2.14 2.94 2.74 3.22

ISABEL W 6.66 3.42 2.53 2.11 2.76 2.47 2.51 2.20 2.34 2.82 2.96

Та 7.14 2.24 1.88 2.29 2.16 2.54 2.79 2.56 2.94 2.74 3.22

СЕМ03.02 W 5.52 3.22 2.88 2.50 2.32 1.99 2.13 2.29 2.87 2.86 3.31

Та 1.61 1.85 2.21 1.59 1.42 2.86 4.17 4.19 4.30 3.43 3.33

INCL4.2+ABLA W 6.99 3.80 2.78 2.57 2.36 2.93 3.46 3.06 4.09 3.78 3.80

Та- 5.83 3.09 2.31 2.17 2.02 2.90 3.84 3.35 3.65 3.68 3.94

INCL4.5+ABLA07 W 6.05 4.34 2.66 2.14 1.73 1.84 2.11 1.76 2.03 1.89 2.01

Та 3.89 2.09 2.33 1.81 1.48 1.59 1.85 1.75 1.80 1.73 2.01

PHITS W 8.41 4.30 3.19 2.61 3.52 2.23 2.58 2.20 2.84 2.88 2.90

Та 2.44 1.91 1.83 2.49 3.25 2.27 2.37 2.08 2.64 2.51 2.81

CASCADE07 W 6.26 3.30 3.24 2.80 2.97 2.45 2.42 2.48 3.25 3.20 3.34

Та 9.43 1.70 2.79 1.75 2.94 3.16 2.65 2.42 3.11 3.17 2.96

>5р,

72Ли 173Аи <1™- 177Т1

,6.05 „2р( ^- ,76Н§ ^- ,!0рЬ

161^__13% ) 160Тд 66% ) 160Н[ 99.3°/^ 1б0^и

81 »1 «1 11

,56Ш |56Ьи ,56УЬ -**» 156Тгп 156Ег

-5 I ^

-4 -1 =4 -4

1вЫ _> 1!2УЬ_> "2Тт-> 153Ег |52Но _«*> |520у

И * г §4 =4

148Тт_> ,48Ег_> |4»Но-> |48Оу-> "*ТЬ-■> ,4,Сс1

£1 А

'Чи

'"Ьи ШУЬ

11 «1

»'Та '"Ш^- "5Оз<- "9Р1<-

Рис. 13. Цепочка распада 1480с1

3Нё

Полученные результаты показывают, что расчетные данные не согласуются с экспериментальными с требуемой точностью ~ 30% (значение фактора <Б> меняется от 1.52 до 9.43), особенно это заметно в области низких энергий.

Примеры измеренных расчетных функций возбуждения и сечений, полученных другими авторами, приведены на рис. 14-16.

к ю

о" ю

к

в

о

¡г

о

О

г - г

/> л- /'" ,1 / .*' """Г"-'

/ * шТа

Энергия протонов, МэВ

Рис. 14. Экспериментальные и расчетные функции возбуждения образования 148Сс1впаЧУи181Та

.'*йо, ......' : \\ г : : 1 1 : ( М1| 1 1 . . М >■) . . ......1 ........... • ■' -' 1/6та» ........' ' ; ** ".'..,-,. I ........1 • • '

.......... ¡V г <ё:' о.«1' / ; У • ...... 1 ....... 1 ' • Рес [ : ' ; • / ■ ...... ' ......1 ' ■ 'ЕЦ ' ...... / Г / • : .Р «' / -• • • / -......... • ' .....1 ' ' "

ю

и ю 2 1

О

10 10 10 10 ю 10

Энергия протонов, МэВ

Рис. 15. Экспериментальные и расчетные функции возбуждения для па'\У(р,х)

I'tfw;" ........1 ' ': ' ч Г | 1 : 1 ■ I Л 'I -». :; \ • . ; \ \ • : ......1 ......., \ ! - Т\ ; ? : f ; i ! < 175l|cl ......1 ' ': г ; • ^ f -O-f*^ о : - ■! : : ................

!'tezrc"' ........■ / ; . L &: Ч ; • г •/ ; 1 : 84RbKm.g)........ '1 Г / ! : 1 Л 1 1 " 1 .......1 , 1 '.: : ВС, ... * Л »' / * / ' * ......1 . '.....1 . ■

2 3 2 3 2 3

10 10 10 10 10 , 10

Энергия протонов, МэВ Рис. 16. Экспериментальные и расчетные функции возбуждения для 181Та(р,х)

Результаты определения фактора для 181 Та и па1\¥, вместе с результатами определения фактора ^ для этих же этих программ, но для других мишеней: ла,Сг, 56Ре, па':№, 93№>, па1РЬ и 209В1 использовались для определения их предсказательной способности в широком диапазоне масс мишеней и энергий налетающих протонов.

Результаты глобального сравнения представлены на рис. 17-23.

40 7(1 1 00 150 250 400 600 000 1200 1600 2601) РиНаи

Рис. 17. Фактор Р для ВЕЯНМ в диапазоне энергий (0.04 - 2.6) ГэВ и диапазоне масс от 50 до 209

40 70 100 1 50 250 400 600 600 1200 ftotort ЕII«'ay (M»V>

Рис. 18. Фактор F для ISABEL в диапазоне энергий (0.04 - 2.6) ГэВ и диапазоне масс от 50 до 209

□ 1-2 »2-3 иЗ-5 «5-10

ft

<

40 70 100 150 250 40Q 600 800 1200 1600 2600 Proton Enerjy

Рис. 19. Фактор F для CEM03.02 в диапазоне энергий (0.04 - 2.6) ГэВ и диапазоне масс от 50 до 209

40 ?0 100 150 250 400 600 f ProtonEiiei{y|MeVl

1200 1600 2600

Рис.20. Фактор Б для РШТБ в диапазоне энергий (0.04 - 2.6) ГэВ и диапазоне масс от 50 до 209

J 1 „ РЬ \ j <

—S | tx / /

J т- ( Г /

ш z KI) J- № fe J

Ш1 / 4

N R1 Г £ — / 4 * J \

ЬО 70 ЮС 150 250 400 600 I Proton Enet gy (Ме\Л

1200 1600 2600

Proton Energy (MeVl

Рис.21. Фактор Б для ШСЬ4.5.5 + АВЬА07 в диапазоне энергий (0.04 - 2.6) ГэВ и диапазоне масс от 50 до 209

Рис. 22. Фактор Б для ШСЬ4+АВЬА в диапазоне энергий (0.04 - 2.6) ГэВ и диапазоне масс от 50 до 209

40 70 100 150 250 400 600 Proton Energy (MeV)

1200 1600 2600

Рис. 23. Фактор Б для СА5САБЕ.07 в диапазоне энергий (0.04 - 2.6) ГэВ и диапазоне масс от 50 до 209

Полученные результаты позволили выделить диапазоны энергий и масс

с недостаточной (меньше фактора 2) предсказательной способностью для каждой из рассмотренных программ, а именно: BERTINI: ~ 2600 МэВ (W, Та), 40 - 100 МэВ (Nb, Fe, Сг), 40 - 70 МэВ (Bi, Pb, W, Та, Nb); ISABEL: 2600 МэВ (W, Ta), 250 - 1200 МэВ (Ta, Nb, Ni, Fe, Cr), 70- 150 МэВ (Nb, Ni, Fe, Cr), 40 -70МэВ (W, Ta, Nb, Ni, Fe, Cr), 40- 100 МэВ (Bi, Pb); CEM03.02: 500 - 2600 (W, Ta), 40 - 100 (Fe, Cr), 40(Bi, Nb), 40 - 70(W); INCL4+ABLA: 400 - 2600 (W,Ta); 150 - 400 (Nb, Ni, Fe, Cr); 40 - 150 (Ni, Fe, Cr); 40-100 (Bi, Pb, W, Ta, Nb); INCL4.4.5 +ABLA07: 40 - 70(W, Ta, Nb, Ni), 40 (Bi); PHITS: 40 - 400(Fe, Cr), 40 - 150 (Nb, Ni), 100 - 250 (W, Ta, Nb), 40-100 (Bi, Pb, W); CASCADE07: 1600 -2600 (Bi, Pb, W, Ta, Nb, Ni, Fe, Cr), 1200 - 1600 (W, Ta, Nb, Ni, Fe, Cr), 250-400 (Fe, Cr), 40 - 250 (Ni, Fe, Cr), 40 - 150 (Nb, Ni, Fe, Cr), 70 - 150 (Bi, Pb, W), 4070 (Bi, Pb, W, Ta, Nb, Ni, Fe, Cr).

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

2) Обосновано использование метода a-спектрометрии для экспериментального определения сечений РПР без предварительной химической сепарации из облученных образцов.

3) Определены сечения РПР для мишенных/конструкционных материалов ЭлЯУ "atW и шТа при энергии протонов 40, 70, 100, 150, 250, 400, 600, 800, 1200, 1600, 2600 МэВ. Всего получено 1013 для natW и 882 сечения для 181Та, которые представлены в виде 192 и 173 функций возбуждения для nalW и 18íTa, соответственно. Из них 99 функции возбуждения для natW и 76 функций возбуждения для 18'Та представлены впервые. Все данные переданы для включения в библиотеку ядерных данных EXFOR. Значения сечений сравнены с данными аналогичных экспериментов, выполненных в других исследовательских группах. В большинстве случаев в ходе сравнения наблюдается хорошая согласованность данных, что свидетельствует о достоверности полученных в диссертации результатов. В то же время результаты сравнения демонстрируют новизну, поскольку количество представленных в ней значений сечений значительно превышает объем

данных других групп.

4) Для natW и шТа проведено сравнение полученных результатов с расчетными, полученными с использованием высокоэнергетических транспортных программ: MCNPX(Bertini, ISABEL), СЕМ03.02, INCL4.2 + ABLA, INCL4.5 + ABLA07, PHITS, CASCADE.07. Сравнения расчетов с экспериментальными данными показывает, что значения фактора <F> лежат в диапазоне от 1.5 до 9.4 для различных программ. Эти значения соответствуют отклонению расчетов от эксперимента на (50 - 800) %. Таким образом, для тяжелых ядер даже для более совершенных программ отклонение превышает требуемую точность 30%. Расхождения особенно велики в области низких энергий.

5) Используя данные сравнения экспериментальных результатов с расчетными, полученными в других работах (natCr, 56Fe, natNi, Nb, " Pb, Bi) для каждой вышеуказанной программы выполнен анализ ее предсказательной способности путем построения двумерных диаграмм для всего энергетического диапазона и всех наборов ядер-мишеней, что дает возможность распространить результаты сравнения на ядра мишени, для которых эксперименты не проводились.

6) Результаты такого глобального анализа подтверждают вывод о том, что все программы на основе модели внутриядерного каскада нуждаются в дальнейшем развитии.

Автор выражает глубокую благодарность:

1) коллегам по лаборатории Фундаментальных ядерно-физических исследований ФГУП «ГНЦ РФ ИТЭФ» за помощь при выполнении экспериментальных и расчетных исследований, описанных в диссертации;

2) сотрудникам синхротрона ФГУП «ГНЦ РФ ИТЭФ» У-Ю за обеспечение 11 типов режимов работы ускорителя, используемых при облучении экспериментальных мишеней;

3) Международному ' научно-техническому центру и Госкорпорации «Росатом» за финансовую поддержку представленных в диссертации исследований.

Личное участие автора. Автор принимал непосредственное участие в качестве полноправного участника научной группы в исследованиях,

описанных в диссертации, на всех этапах их выполнения: анализе литературных данных, создании экспериментальной методики определения а-сечений, изготовлении экспериментальных образцов и их компоновкой для облучения, калибровке используемой аппаратуры, измерения а- и у-спектров облученных образцов и мониторов, обработке полученных а- и у-спектров и вычислении значений сечений и потоков. Проведение расчетов с использованием программ NCNPX(Bertini, ISABEL), INCL4.2, INCL4.5. Анализ полученных экспериментальных и расчетных результатов и верификация высокоэнергетических транспортных программ были проведены автором совместно с научным руководителем.

Апробация работы. Результаты работы докладывались на: 1) Сессии-конференции секции ЯФ ОФН РАН 2007 «Физика фундаментальных взаимодействий», ФГУП «ГНЦ РФ ИТЭФ», октябрь, 2007, Москва, Россия;

2. Международной конференции «NUFRA 2007 - International Conference on Nuclear Fragmentation», 24 сентября - 1 октября, 2007, Анталия (Кемер), Турция;

3) Первом международном рабочем совещании «Accelerator Radiation Induced Activation ARIA 2008», Paul Scherrer Institute,13-17 октября, 2008, Виллиген, Швейцария;

4) Международной конференции «Shielding Aspect of Accelerators, targets and irradiation Facilities (SATIF-9)», Oak Ridge National laboratory, 21-23 апреля, 2008, Окридж, США;

5) Научной сессии НИЯУ МИФИ - 2010, январь, 2010,Обнинск, Россия;

6) Рабочей встрече участников Координационного Исследовательского Проекта «Анализ Ядерных Моделей и Программ», СБА, февраль, 2010, Saclay, Париж, Франция;

7) Научной сессии НИЯУ МИФИ - 2011, февраль, 2011,Обнинск, Россия.

Публикации автора по теме диссертации 1)Ю.Е. Титаренко, В.Ф. Батяев, ... С.Н Флоря и др. Измерение и моделирование сечений образования продуктов реакций при облучении natW и шТа протонами энергиями (0.04 - 2.6) ГэВ, Ядерная физика, 2011, т. 74, №4, с. 574-595.

2)Ю.Е. Титаренко, В.Ф. Батяев, ... С.Н Флоря и др. Сечения мониторных реакций 27А1(р, x)24Na, 27А1(р, x)22Na, 27А1(р, х)7Ве при энергиях протонов (0.04 - 2.6) ГэВ, Ддерная физика, 2011, т. 74, №4, с. 531-547.

3) Ю.Е. Титаренко, В.Ф. Батяев, ... С.Н Флоря и др. Измерение и моделирование сечений образования I48Gd в тонких n3tW и шТа мишенях при воздействии протонов с энергиями в диапазоне (0.4 - 2.6) ГэВ, Ядерная физика, 2011, т. 74, №4, с. 596-604.

4)Titarenko Yu. Е., Batyaev V. F., ... Florya S. N. et al. 148Gd formation in thin na,W and 181Ta targets induced with (0.4 - 2.6) GeV. First International Workshop on Accelerator Radiation Induced Activation ARIA 2008, October 13-17, 2008, Paul Scherrer Institute, Villigen, Switzerland, PSI Proceedings 09-01, January 2009, ISSN 1019-0643, p. 68-74.

5) Titarenko Yu. E., Batyaev V. F..... Floiya S. N. et al. Experimental and

theoretical study of 148Gd formation in thin natW targets induced with (0.4 - 2.6) GeV proton. Shielding Aspect of Accelerators, targets and irradiation Facilities (SATIF-6). 21-23 April 2008, Spallation Neutron Source, SNS, 701 Scarboro Road, Oak Ridge, TN, USA Proceedings NEA/NSC/DOC (2010)7, p. 181-186.

6)Titarenko Yu. E., Batyaev V. F., ... Florya S. N. et al. Nuclear Data Requirements for Accelerator-Driven Systems. NUFRA 2007 - International Conference on Nuclear Fragmentation, Antalya (Kemer), September 24-October 1, 2007, http://fias.uni-frankfurt.de/nufra2007/talks/titarenko.pdf.

Подписано к печати 20.07.11 г. Формат 60x90 1/16

Усл. печ. л. 4,0 Уч.-изд. л. 0,68 Тираж 100 экз. Заказ 578

Отпечатано в ИТЭФ, 117218, Москва, Б.Черемушкинская, 25