Радиохимия, т.42, №2 (2000) 180-187

УДК 539.3: 621.039.8

КОМПЛЕКСНАЯ ДИФФУЗИОННО-ЗОНДОВАЯ ТОМОГРАФИЯ ГЕТЕРОГЕННЫХ СРЕД

И.Н.Бекман

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

Перспективы применения микротомографии на радиоактивных диффузионных газовых зондах для диагностики неоднородных материалов связаны с различием в механизмах и интенсивности процессов растворения и диффузии радиоактивных газов в отдельных компонентах (структурных образованиях) исследуемой гетерогенной среды [1-4]. В неразрушающих вариантах диффузионно-зондовой томографии (ДиТ) за развитием процесса миграции радионуклида в твердом теле следят по изменению состояния радиационного поля вокруг образца. Из пространственного распределения интенсивности излучения, испускаемого диффундирующей примесью, его энергии и вариаций этих величин во времени рассчитывают нестационарные концентрационные профили зонда, из которых, в свою очередь, определяют механизм массо-переноса и значения коэффициентов диффузии и растворимости зонда в отдельных компонентах гетерогенной структуры и на границах их раздела. Из спектров локальных коэффициентов диффузии и констант растворимости находят пространственное расположение компонентов структуры, транспортные связи между ними и их пропускную способность.

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

В качестве примера рассмотрим проблему применения ДиТ для пространственной локализации двух адсорбционно активных центров (точечных или в виде включений в виде сфер или цилиндров), обнаружения наличия транспортных связей между ними и оценки их пропускной способности.

Прежде всего попытаемся использовать для структурных целей метод динамической авторадиографии, позволяющий снимать последовательные кадры развития процесса диффузии [5, 6]. Как известно, для повышения разрешающей способности АРГ следует использовать мягкое излучение, например, a-частицы. Но тогда возникает проблема возможной глубины анализа, поскольку невозможно обнаружить источники излучения, расположенные на глубине большей, чем пробег a-частиц в исследуемом материале.

В качестве примера рассмотрим диагностику двух точечных адсорбционно активных центров, между которыми возможно существование транспортной связи. Пусть в нашем распоряжении имеются образцы в виде пластины, в каждом из которых расположены по два точечных дефекта, способные активно сорбировать радиоактивный газ (адсорбционно активные центры). Растворением и диффузией зонда в материале изделия в первом приближении пренебрегаем. Возможны две ситуации: 1) дефекты не связаны друг с другом, но каждый имеет свой транспортный путь на внешнюю поверхность образца, по которому возможно его насыщение или дегазация; 2) один дефект (№1) имеет транспортную связь со вторым (№2), от которого диффузионный путь ведет на поверхность (дефект №1 не имеет транспортной связи с поверхностью). Каждый дефект характеризуется собственной емкостью по отношению к тестовому газу. Пропускная способность каждого транспортного пути определяется собственным значением константы скорости диффузии:

где D - коэффициент диффузии зонда по транспортному пути, L -длина транспортного пути (При последовательной диффузии L= l₁+ l₂, где l₁-длина диффузионного пути, связывающего дефекты №1 и №2, l₂ - длина транспортного пути, связывающего дефект №2 с поверхностью образца).

Из-за ограниченности пробега a-частиц в материале нам придется рассмотреть набор из шести образцов. В образцах 1 и 2 дефекты не связаны друг с другом, а в образцах 3-6 связаны друг с другом. Во всех образцах дефект №2 расположен на глубине (h₂=L₂=15 мкм) существенно меньшей, чем величина пробега R_a=50 мкм, и a-излучение хорошо регистрируется детектором. Дефект №1 в образцах 1, 3 и 5 расположен на глубине hR_a (см. рис.1а). Полагали, что емкости дефектов (максимально возможная концентрация зонда в адсорбционно активном центре) и коэффициенты диффузии зонда по обоим транспортным каналам одинаковы (D=10^-7 см²/с). Схемы расположения дефектов представлены на рис.1а, а количественные характеристики собраны в табл. 1.

Табл.1 Параметры адсорбционно активных центров

Регистрацию излучения, покидающего пластину, будем осуществлять с помощью динамического варианта трековой авторадиографии. Зондирование структуры проводим с использованием a-излучающего радиоактивного газа. Рассмотрим режим сорбции. Будем полагать, что образец в виде пластины насыщается радиоактивным газом из постоянного источника. В некоторый момент времени диффузию останавливают, образец извлекают из резервуара, и проводят авторадиографирование со входной по отношению к диффузанту стороны пластины. (Для регистрации распределения треков a-частиц по поверхности образца используется фотоматериал со слоем фотоэмульсии заведомо превышающий пробег в нем регистрируемого излучения). Затем образец вновь помещают в резервуар с радиоактивным газом и продолжают процесс диффузии. Последовательные "кадры" дифузионного фильма снимают до полного насыщения адсорбционных центров.

В ходе математического эксперимента расчеты трековых авторадиограмм от точечного источника a-излучения проводили по формулам:

где h - глубина залегания точечного источника, x₁, y₁ - координаты начала проекции на АРГ трека, x₂, y₂ - координаты конца проекции трека (предполагается, что сам источник находится в точке с координатами x = y = 0), R_a - пробег a-частиц (пробеги в материале образца и фотоэмульсии считали одинаковыми), q -угол вылета a-частицы, отсчитываемый от перпендикуляра к плоскости АРГ, a - угол направления трека в плоскости наблюдения. В ходе расчета АРГ углы q и a меняли по закону случайных чисел, имеющих равномерное распределение.

Трековые авторадиограммы от всех образцов, зарегистрированные при равновесной сорбции, приведены на рис.1б. Видно, что АРГ фиксирует пространственное расположение источников в плоскости x, y, но не всегда достаточно однозначно и не всех: из шести образцов, только на двух авторадиограммах 1б и 3б представлены две "звезды", что соответствует реальности. "Звезды", связанные с дефектом №1 на АРГ 2б, 4б и 6б отсутствуют, т.к. адсорбционно активный центр расположен на глубине, больше пробега a-частиц и излучение от адсорбированного на нем газа не достигает фотопластинки. АРГ образца 5 фиксирует оба дефекта в виде одной "звезды" (сумма звезд, представленных на рис.1б), т.к. источники излучения находятся точно друг под другом. Таким образом, визуальное рассмотрение авторадиограмм позволило только в двух случаях из шести однозначно и правильно определить число дефектов в материале и провести их пространственную локализацию.

Рис. 1. Комплексная диффузионно-зондовая томография точечных адсорбционно активных центров (дефектов).
а) Схема расположения адсорбционно активных центров и транспортных путей в образцах.
б) Трековые авторадиограммы распределения a-излучающего радиоактивного зонда, адсорбированного в дефектах (Режим равновесной сорбции, авторадиографирование входной относительно тестового газа поверхности пластины, одинаковая сорбционная емкость всех дефектов).
в) Спектры a-излучения, покидающего образец.
г) Кинетические кривые сорбции (пунктир) и десорбции (сплошные линии) радиоактивного зонда (сорбция/десорбция на отдельных дефектах, расположенных в зоне пробега, и суммарная сорбция/десорбция образцом).

Дополнительную информацию можно получить из снятия энергетического спектра a-излучения, исходящего из входной по отношению к диффундирующему зонду поверхности пластины (или из анализа распределения проекций треков по длине, что требует трудоемких расчетов пространственных углов расположения треков и длин проекций). Как известно, вид спектра a-частиц, зарегистрированный от некоторой концентрации радиоактивного газа на глубине h, С(t) изменяется во времени (в процессе диффузии) по закону:

где Е - энергия a-частицы, покидающей образец, R_e - пробег a-частицы в фотоэмульсии (полный пробег R_a = R_e+R_m, где R_m пробег a-частицы в исследуемом материале) а - коэффициент связи между энергией излучения и его пробегом в материале (Е = аR_e), V- объем дефекта, С - концентрация находящегося в дефекте газа, А - удельная радиоактивность газа, h - глубина залегания источника. Диапазоны изменения параметров: 0 < h <R_a, 0 £ q £ arcCos(h/R) или

На рис.1в приведены энергетические спектры всех рассматриваемых здесь шести образцов, зарегистрированные после окончания диффузии (при равновесной сорбции), когда установлено равновесие между газовой и твердой фазами. При расчетах полагали, что a-частица с энергией 5 Мэв приводит к пробегу в исследуемом материале 50 мкм, т.е. а = 0.1 МэВ/мкм. Из рисунка видно, что из-за ограничений, связанных с пробегом, спектры 2в, 4в и 6в дают одну линию, что совершенно не соответствует действительности. Однако, упомянутые выше трудности в интерпретации образца 5, возникшие из-за расположения источников излучения точно друг под другом, на спектре 5в однозначно разрешаются: на "экспериментальном" спектре две линии четко демонстрируют наличие двух источников излучения. Спектры 1в, 3в и 5в позволяют по максимальной энергии определить глубину залегания каждого дефекта, а по высоте пика - концентрацию зонда в каждом адсорбционном центре.

Таким образом, простейший пример двух адсорбционно активных центров с различным расположением в пространстве и с разной эффективностью транспортной связи между ними показывает, что одновременное использование трековой авторадиографии и ядерной спектроскопии для режима равновесной сорбции в трех ситуациях из шести (там, где дефекты расположены в зоне пробега) позволяет однозначно определить число дефектов, найти все их пространственные координаты, емкости по тестовому газу, а также длины диффузионных путей L₁ и L₂. Дефекты, залегающие глубже зоны пробега здесь себя никак не проявляют и догадаться об их существовании невозможно.

Перейдем теперь к рассмотрению перспектив использования процессов сорбции и десорбции для сортировки дефектов по степени их связности друг с другом.

Кинетика сорбции и десорбции газа в двух связанных друг с другом дефектах в предположении кинетического механизма последовательной диффузии описывается решением системы уравнений:

	(3а)
	(3б)

Здесь С₁ и С₂ концентрации радиоактивного газа, адсорбированные в дефектах 1 и 2, соответственно, _t - время диффузии.

Система уравнений (3) решается при условиях:

Решение можно представить в виде:

	(4а)
	(4б)

где  ;

Суммарная кинетика сорбции/десорбции задается формулой:

(5)

Следствия:

В случае сорбции в первоначально пустые ловушки:

В случае десорбции до полного опустошения ловушек

где

Если дефекты не связаны друг с другом, то концентрации зонда в них в ходе сорбции (кинетический режим параллельной диффузии) описываются формулами:

	(6а)
	(6б)

а в ходе десорбции

	(7а)
	(7б)

Рис. 2. Нестационарная сорбция и десорбция тестового радиоактивного газа образцом, содержащим два адсорбционно активных центра. (Одинаковая сорбционная емкость всех дефектов).
а) Обыкновенный масштаб
б) Логарифмический масштаб
Пунктир - сорбция, сплошные линии - десорбция.
Цифры на кривых - номера образцов (см. рис. 1а).

Примеры изменения концентрации зонда в каждом дефекте (измеряется по площади соответствующего пика в a-спектре) и суммарной концентрации зонда (измеряется по всей площади спектра) во времени в ходе экспериментов сорбции и десорбции для всех образцов приведены на рис.1г. Видно, что в одних случаях (рис.1г, 1, 3, 5) число кривых соответствует числу адсорбционно активных центров, в других (рис.1г, 2, 4, 6) - кинетика диффузии регистрируется только для дефектов, расположенных в зоне пробега. По форме кинетической кривой легко различить случаи параллельной и последовательной диффузии.

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

На рис.2а приведены временные зависимости количества газового зонда, прошедшие через поверхность пластины в режимах сорбции и десорбции для всех шести образцов (Рис.1а и Табл.1). Расчеты проводили в предположении равенства адсорбционной емкости дефектов. На рис.2б те же результаты представлены на рис.2б в логарифмическом масштабе. Для шести образцов мы получили пять кривых сорбции и пять кривых десорбции (кривые для образцов 3 и 5 совпали, т.к. здесь 1= 2). Кинетика диффузии в образце 2 для дефекта N2 описывается одной экспонентой. Здесь присутствие дефекта N1, лежащего вне зоны пробега, никак не проявляется. Остальные кривые имеют аномальную форму (результаты не спрямляются в логарифмическом масштабе). Математический анализ формы интегральной кинетической кривой позволяет отсортировать пары связанных и не связанных дефектов. Особенно эффективным является сравнение кинетических кривых сорбции и десорбции. С этой целью построим на одном графике кривые десорбции и кривые сорбции, пересчитанные по формуле:

Из рис.3 видно, что для образцов с несвязанными друг с другом дефектами (параллельная диффузия) кривые сорбции и десорбции полностью совпадают. Для "гирлянд" связанных друг с другом дефектов (последовательная диффузия) кривые сорбции и десорбции не совпадают друг с другом (десорбция опережает сорбцию). Данный эффект является надежным критерием наличия в образце сети связанных дефектов. Разложение сложной кинетической кривой и нахождение отдельных вкладов центров адсорбции и диффузионных путей для образцов 1 и 2 представляет собой достаточно простую задачу (чем сильнее различаются длины диффузионных путей, тем выше точность расчетов). Для последовательной диффузии задача усложняется, особенно при малых отличиях в длинах диффузионных путей. Отметим, что метод десорбции в этом случае существенно более информативен, чем метод сорбции (см. рис.3б где для образцов 3-6 кривые десорбции имеют явно выраженные изгибы. тогда как сорбция описывается практически линейными зависимостями. Важно подчеркнуть, что в случае связанных друг с другом дефектов удается обнаружить наличие центров адсорбции на любой глубине (в том числе - вне зоны пробега). Количественная обработка кинетических кривых позволяет рассчитать константы диффузии, а затем (при известном коэффициенте диффузии) найти длины соответствующих транспортных путей.

Рис. 3. Сравнение кинетических кривых сорбции и десорбции в предположении одинаковой емкости всех дефектов.
а) Обыкновенный масштаб
б) Логарифмический масштаб
Пунктир - сорбция, сплошные линии - десорбция.
Цифры на кривых - номера образцов (см. рис. 1а).

Коротко остановимся на осложнениях, связанных с анализом дефектов с существенно различной емкостью. Пусть емкость более близкого к поверхности образца дефекта №2 много меньшей емкости дефекта №1 (С10=200, С20=20). Кинетики сорбции и десорбции для этого случая приведены на рис.4. Видно, что построенные в логарифмическом масштабе зависимости (рис.4б) почти прямолинейны, разложить их на отдельные составляющие практически невозможно и приходится ограничиваться расчетом некоторых "эффективных" констант скоростей диффузии, что существенно снижает эффективность методики. Существенно лучшие результаты можно получить, если кинетики сорбции и десорбции регистрировать раздельно для каждого дефекта (см. ниже).

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

Очевидно, что временная развертка позволяет существенно повысить разрешение энергетических спектров испускаемых из образца a-частиц. С точки зрения определения числа дефектов важное значение имеют временные изменения каждой линии a-спектра. Возникающие здесь проблемы иллюстрирует рис.1г, на котором приведены кинетики локальной и полной сорбции/десорбции из всех образцов. Особое внимание следует уделить анализу образцов 2, 4 и 6, когда в наличии имеется два источника излучения, но в спектре - лишь одна линия (связанная с излучением только от адсорбционного центра №2, расположенного в зоне пробега излучения, по изменению во времени интенсивности которого и построены кривые рис.4). Видно, что кинетические кривые для образцов 4 и 6 имеют аномальную форму, анализ которой позволяет обнаружить присутствие глубоко расположенного (невидимого) дефекта. Дефекты на глубинах вне зоны пробега обнаруживают себя, если они связаны с поверхностью через дефект в зоне пробега. В этом случае, из анализа формы локальной кинетической кривой для центра №2 можно рассчитать параметры kD₁ и kD₂ и, зная коэффициент диффузии D, определить длины диффузионных путей L₁ и L₂ (истинную глубину залегания дефекта №1, h₁, найти невозможно). Если дефекты не связаны друг с другом и дефект №1 имеет прямой выход на поверхность, то обнаружить его присутствие ни методом динамической авторадиографии, ни разверткой энергетических спектров во времени нельзя. Тем не менее, рассмотренные здесь методики позволили достаточно подробно разобраться в пяти ситуациях из шести

Рис. 4. Нестационарная сорбция и десорбция тестового радиоактивного газа образцом, содержащим два адсорбционно активных центра различной емкости. (Сорбционная емкость дефекта №2 (рис.1а) в десять раз меньше емкости дефекта №1).
а) Обыкновенный масштаб
б) Логарифмический масштаб
Пунктир - сорбция, сплошные линии - десорбция.
Цифры на кривых - номера образцов (см. рис. 1а).

Адекватный анализ образца 2 (два несвязанных друг с другом дефекта, один из которых расположен вне зоны пробега) можно провести, если обратиться к суммарной кинетической кривой десорбции, измеренной по увеличению активности газовой фазы. Сравнение ситуаций образцов 2 и 6 (топология дефектов одинакова, но в образце 2 они не связаны, а в образце 6 связаны друг с другом) показывает, что кинетические кривые отличаются как друг от друга, так и от ситуации с одним дефектом. Анализ подобных кривых позволяет определить число дефектов (вне зависимости от глубины залегания) и найти параметры k₁ и k₂

Рис. 5. Нестационарная сорбция и десорбция радиоактивного зонда на адсорбционно активном центре, ближайшем к поверхности (дефект №2, рис.1а).
а) Одинаковая сорбционная емкость дефектов
б) Сорбционная емкость дефекта №2 (рис.1а) в десять раз меньше емкости дефекта №1
Пунктир - сорбция, сплошные линии - десорбция.
Цифры на кривых - номера образцов (см. рис. 1а).

Как уже упоминалось, при произвольном отношении емкостей адсорбционных центров, задача количественной обработки экспериментальных данных существенно усложняется. В этой связи рассмотрим перспективы использования для этих целей локализующей спектроскопии. На рис. 5а приведены кривые сорбции и десорбции только для дефекта N2. Видно, что при одинаковой емкости ловушек (рис.5а) кинетические кривые мало информативны и для анализа предпочтительнее использовать кинетики, измеренные для всего образца (рис.2). При различной емкости ловушек наоборот, суммарные кинетические кривые неинформативны (рис.4б), но зато весьма наглядны кривые десорбции для одного центра N2 (рис.5). Обнаружить механизм последовательной диффузии можно уже визуально: экспоненциально ниспадающие зависимости сменяются на пикообразные. Это прекрасно видно, когда все дефекты расположены в зоне пробега. Но даже если дефект N1 лежит за пределами пробега, то его присутствие существенно и наглядно искажает форму кинетической зависимости, измеренной для дефекта N2 (см. кривую 6 на рис.5б).

Подводя итоги можно сказать, что комплекс методик позволил для всех образцов найти такие параметры, как число центров адсорбции, их предельные емкости по тестовому газу, константы и коэффициенты диффузии по транспортным путям, а также длины транспортных путей. Для случаев залегания всех дефектов в зоне пробега удалось рассчитать все три координаты их расположения. Для трех образцов определена локализация лишь источника №2, тогда как местоположение источника №1 (залегающего вне зоны пробега) найти не удалось. Отметим, что ситуация с глубоко залегающими дефектами не безнадежна. Поскольку по косвенным данным доказано их существование и даже оценены некоторые важные параметры, то для полной диагностики образцов с такими дефектами нужно провести аналогичные эксперименты с радионуклидом, обладающим более жестким излучением. Погрешность локализации дефекта 1 новым зондом будет выше, чем ошибка локализации дефекта 2 старым зондом с низкоэнергетическим излучением, но все же координаты расположения глубокого дефекта будут определены.

Аналогичный подход можно применить и к адсорбционно активным неоднородностям структуры, имеющим определенную форму и размеры и сложным образом расположенные в пространстве. Некоторые примеры приведены на рис.6. Видно, что авторадиограммы для a-излучающего зонда и равновесного режима сорбции (рис.6б) позволяют хорошо различить слоистые среды от отдельных включений. Координаты центров включений в плоскости x,y определяются достаточно надежно. Соосные включения (образец 5 рис.6) регистрируются как одно скопление треков. Ядерная спектроскопия (рис.6в) позволяет с трудом отличить слоистые среды от включений типа цилиндра (вертикальный цилиндр и толстый слой дают практически одинаковую форму спектра). Спектры полезны в определении глубины залегания отдельных включений, их размеров, формы и емкости по тестовому газу. Энергетические спектры особенно эффективны для измерения глубины залеганий скоплений и разграничения соосно расположенных центров адсорбции. Очевидно, что при увеличении плотности неоднородностей структуры и большом разнообразии их форм, задача интерпретации результатов томографии существенно усложняется. Одним из примеров может служить образец 4 рис.6.

Рис. 6. Комплексная диффузионно-зондовая диагностика материалов с различными типами неоднородности структуры.
а) Топология неоднородностей структуры (форма источника излучения): вертикальный цилиндр (1), слой, прилегающий к поверхности пластины (2), сфера (3), слой, сфера и цилиндр (4), соосно расположенные сферы (5).
б) Трековые авторадиограммы распределения a-излучающего радиоактивного зонда, адсорбированного в дефектах (Режим равновесной сорбции, авторадиографирование входной относительно тестового газа поверхности пластины, одинаковая сорбционная емкость всех дефектов).
в) Спектры a-излучения, покидающего образец

Аналогичный подход можно применить и к адсорбционно активным неоднородностям структуры, имеющим определенную форму и размеры и сложным образом расположенные в пространстве. Некоторые примеры приведены на рис.6. Видно, что авторадиограммы для a-излучающего зонда и равновесного режима сорбции (рис.6б) позволяют хорошо различить слоистые среды от отдельных включений. Координаты центров включений в плоскости x, y определяются достаточно надежно. Соосные включения (образец 5 рис.6) регистрируются как одно скопление треков. Ядерная спектроскопия (рис.6в) позволяет с трудом отличить слоистые среды от включений типа цилиндра (вертикальный цилиндр и толстый слой дают практически одинаковую форму спектра). Спектры полезны в определении глубины залегания отдельных включений, их размеров, формы и емкости по тестовому газу. Энергетические спектры особенно эффективны для измерения глубины залеганий скоплений и разграничения соосно расположенных центров адсорбции. Очевидно, что при увеличении плотности неоднородностей структуры и большом разнообразии их форм, задача интерпретации результатов томографии существенно усложняется. Одним из примеров может служить образец 4 рис.6.

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

Литература.

1. Бекман И.Н. // Вестн. МГУ. Сер. 2: Химия. 1994. Т.35, N6, С. 547.
2. Бекман И.Н. // Радиохимия. 1998. Т.40. N2. С. 150.
3. Бекман И.Н., Бунцева И.М. // Радиохимия. Т.40. N2. С. 160.
4. Бекман И.Н., Бунцева И.М. // Радиохимия. Т.40, N2. С. 168.
5. Бекман И.Н. // Радиохимия. 1981. Т. 23, N 5. 750.
6. Бекман И.Н., Швыряев А.А. // Радиохимия. 1982. Т. 24. N 1. С.126.