Безуглов М.В., Малышевский В.С., Малыхина Т.В., Торговкин А.В., Фомин Г.В., Шраменко Б.И. Фоторождение космогенного бериллия-7 в земной атмосфере
Научная статья
Электронный научный журнал "Исследовано в России"а 589 атмосфере
Безуглов М.В. (1), Малышевский B.C. (vsmalyshevsky@sfedu.ru) (1), Малыхина Т.В. (3), Торговкин А.В. (2), Фомин Г.В. (1), Шраменко
Б.И.(2)
(1) Южный федеральный университет, 344090, Ростов-на-Дону, Россия (2) Национальный научный центр Харьковский физико-технический институт,
61108, Харьков, Украина (3) Харьковский национальный университет имени В.Н.Каразина, 61077, Харьков,
Украина
Введение
На сегодняшний день мониторинг содержания радионуклидов в приземном слое атмосферы позволяет сделать вывод о том, что существенный вклад в радиоактивность приземного воздуха дает короткоживущий изотоп Be космогенного происхождения. Вариации содержаний Be в воздухе связаны с солнечной активностью и имеют характерный сезонный ход и широтную зависимость [1,2]. Благодаря быстрому распаду его активность изменяется в растениях в зависимости от синоптических условий. Поэтому Be представляет интерес не только с точки зрения радиоактивного воздействия на биологические системы, а также может являться индикатором скоростей обмена в растениях, и, как следствие, показателем накопления природными средами загрязняющих веществ, поступающих из атмосферы [3]. Именно это делает его удобным индикатором для быстрой оценки возможного атмосферного загрязнения и воздушного обмена в окружающей природной среде. Поэтому изучение механизмов и закономерностей процессов возникновения, переноса и миграции радионуклида Be в объектах экосферы и на границах их раздела весьма актуально и своевременно.
Считается, что основные реакции, приводящие к образованию изотопов бериллия в атмосфере Земли, протекают при взаимодействии космических лучей с ядрами азота и кислорода [4], которые являются основными составляющими атмосферного воздуха. Это
А."71йаа "7аа 1 .аа "7
так называемые реакции скалывания 7N (р,Х) 4Ве , gO (р,Х) 4Ве , 7N (п,Х) 4Ве и 80 (п,Х) 4Ве . Другим возможным механизмом образования изотопа бериллия Be в верхних слоях атмосферы могут быть фотоядерные реакции, которые до настоящего времени не учитывался.
Для расчета скорости образования изотопов Be в земной атмосфере благодаря этому механизму необходимо знать сечения фотоядерных реакций. К сожалению, в литературе очень мало данных по реакциям многочастичного фоторасщепления ядер с
Электронный научный журнал "Исследовано в России"а 590 7 17аа 7а Д7 1а 7
Сечение фоторождения Beа в реакциях 6Са (у,Х)4Ве , 7Nа (у,Х)4Ве , sOа (у,Х)4Ве
1аа 19а 7а 14а 1
Для определения сечения фоторождения Be в реакциях бС (у,Х)4Ве , 7N (у,Х)4Ве , sO (у,Х)4Ве в области гигантского дипольного резонанса был проведен эксперимент по облучению тормозным излучением на линейном ускорителе электронов ЛУ-40 ряда мишеней, включающих кислород, азот и углерод. Энергия ускоренных электронов составляла 90 МэВ, ток 3,26 мкА. Химическими веществами, содержащими данные элементы, были выбраны: нитрид бора (BN), нитрид алюминия (A1N), корунд ( АЬОз), пирографит (С). В качестве тормозного конвертора использовались четыре танталовые (Та) пластины толщиной 1 мм, разделенные промежутком 1 мм. Для вычисления сечения фоторождения изотопа 4В е необходимо знание плотности потока тормозных фотонов в месте расположения мишеней. Для этой цели было проведено компьютерное моделирование прохождения первичных электронов с энергией 90 МэВ через модель мишенной сборки, имеющую параметры, которые соответствуют параметрам экспериментальной мишенной сборки. Для проведения расчетов была разработана компьютерная программа "КГРТ". Программа разработана на языке С с использованием библиотеки классов Geant 4 версии 9.4 (декабрь 2010).
Диаметр пучка электронов задавался равным 5 мм, геометрия пучка электронов задавалась с помощью класса G4UniformRand. Параметры модели мишенной сборки описаны с помощью методов класса G4DetectorConstruction (определены составные части модели мишенной сборки, в частности, определены геометрические параметры и материалы, параметры визуализации и др.). В таблице 1 представлены основные геометрические параметры составных частей модели мишенной сборки (все составные части имеют цилиндрическую форму). На рис. 1 представлен схема облучения и устройство мишенной сборки
Электронный научный журнал "Исследовано в России"а 591 сборки
№ |
Материал |
Толщина, мм |
Диаметр, мм |
Комментарий |
0 |
Ti |
0.05 |
19 |
p = 4.5 г/см3 |
1 |
Air |
40.0 |
19 |
p = 1.29 мг/см3 |
2 |
Та |
1.0 |
19 |
p = 16.62 г/см3 |
3 |
Air |
1.0 |
19 |
p = 1.29 мг/см3 |
4 |
Та |
1.0 |
19 |
p = 16.62 г/см3 |
5 |
Air |
1.0 |
19 |
p = 1.29 мг/см3 |
6 |
Та |
1.0 |
19 |
p = 16.62 г/см3 |
7 |
Air |
1.0 |
19 |
p = 1.29 мг/см3 |
8 |
Та |
1.0 |
19 |
p = 16.62 г/см3 |
9 |
Air |
2.0 |
19 |
p = 1.29 мг/см3 |
10 |
Al |
0.05 |
19 |
p = 2.7 г/см3 |
11 |
Air |
0.95 |
19 |
p = 1.29 мг/см3 |
12 |
Pb |
0.02 |
19 |
p = 11.35 г/см3 |
13 |
Air |
0.98 |
19 |
p = 1.29 мг/см3 |
14 |
Al |
0.05 |
19 |
p = 2.7 г/см3 |
15 |
Al |
0.5 |
19 |
p = 2.7 г/см3 |
16 |
A1N |
3.0 |
19 |
порошок р= 1 г/см3 |
17 |
Al |
2.5 |
19 |
р = 2.7 г/см3 |
18 |
BN |
3.0 |
19 |
порошок р= 1 г/см3 |
19 |
Al |
2.0 |
19 |
р = 2.7 г/см3 |
20 |
A1203 |
1.0 |
19 |
р = 3.97 г/см3 |
21 |
С |
3.0 |
16 |
р = 2.265 г/см3 |
22 |
Air |
1.0 |
16 |
р= 1.29 мг/см3 |
Программа KIPT имеет два режима работы - пакетный и интерактивный. Пакетный
режим работы программы необходим для моделирования прохождения через установку
большого количества первичныхаа электронов (с целью набора статистики).
Интерактивный режим программы имеет модуль визуализации для наиболее полного представления о взаимном расположении составных частей установки. При визуализации прохождения частиц через установку можно использовать режим визуализации поверхности или контурный режим. Визуализация установки программно реализована с использованием библиотеки OpenGL.
Электронный научный журнал "Исследовано в России"а 592
Гаа аа в cd
Рис. 1 Схема облучения мишенной сборки на ускорителе ЛУ-40. Та-конеертер из 4- х пластин тантала;А,В -контейнеры из А1, заполненные порошкообразными нитридами AIN и BN; С - мишень из корунда (А1203);
D- мишень из углерода.
На рис. 2 представлены скриншоты программы в интерактивном режиме, на которых можно видеть модель мишенной сборки (все компоненты отображаются в контурном режиме), траекторию первичного электрона с энергией 90 МэВ, тормозные гамма-кванты (ярко-зеленый цвет) и т.д.
О |
Fri Feb 25.а 7:54 PM |
ewer-0 (OpenGLIniniediaieX) |
Hie Edit View Terminal nextа majorаа release.аа И *** Thisа isа justа a wa rningа message. C4V1sManager:аа Usingа С 4TrajectoryDrawByChar geа asаа default traject oryа model. Seeа commandsаа inа /vis/ modeli ng/traj ectories / forа otherа options. Idle>а /run/beamOnаа 1 Idle>а /vis/viewer/ref resh Idle:- /run/beamOn 1 Idle> /run/beamOn 1 Idle> /run/beamOn 1 Idle> /run/beamOn 1 Idle> Q |
'ЦApplicationsаа Actions V5<SQ^ 0
\*$\а \ШtetG49+^localho5t:-/g494/geа i [m: - -/g494/geant4.9.4/myj|_] viewer-0 (OpenGLIrnmediateX!аа J_J|
Рис. 2. Интерактивный режим программы.
Для описания физических процессов использовалась модель низких энергий электромагнитных процессов (используются данные EPDL97, EEDL97). При моделировании отслеживались гамма-кванты с энергией более 20 МэВ, прошедшие через всю сборку, включая мишень АЬОз (корунд), и достигшие С (углерод). Энергии гамма-
Электронный научный журнал "Исследовано в России"а 593 0аа 10аа 20аа 30аа 40аа 50аа 60аа 70аа 80аа 90
-_ 0,050 '-_ 0,045 0,040 0,035 0,030 -_ 0,025 '-_ 0,020 0,015 0,010 |
geant4 9.4 simulation |
7x105
Energy of primary electrons E =90 MeV number of primary electrons: ft , =2.3*10 1аа 'аа pupnrc |
6x105 5x105 4x105
чЧ'
------ f(x) = 1/x |Ч ХЧ Ngammal |
I 3x105
о
2x105
1х105
0
ПППЧ |
жI..................... I...................... I....................... I...................... I...................... I.................... |
20аа 30аа 40аа 50аа 60аа 70аа 80аа 9сГ"
Еу, MeV
Рис. 3. Смоделированный спектр гамма-квантов, прошедших через сборку, использовалась модель низких
энергий электромагнитных процессов в Geant4
На данном графике (рис. 3) представлены спектры (синие точки на графике) гамма-квантов с энергией от 20 до 90 МэВ, а также график функции f(x)=l/x (красным цветом). Энергияаа первичныхаа электроноваа Ееаа =аа 90аа МэВ,аа количествоаа первичныхаа электронов
NeventS=2.3xl07.
Также было проведено моделирование прохождение электронов с энергией Ее = 90 МэВ через данную установку, но физические процессы для электронов, позитронов и гамма-квантов описывались с помощью классов библиотеки Geant4, использующих модель Penelope. В таблице 2 представлен перечень классов библиотеки Geant4, используемых для первого и второго способов моделирования.
На рис. 4 представлен смоделированный спектр гамма-квантов с энергией выше 20 МэВ, прошедших через сборку (и вошедших в С ). Энергия первичных электронов Ее=90 МэВ, число первичных электронов NeVents=5661660. Диаметр пучка = 5 мм.
Из результатов моделирования видно, что модель низких энергий электромагнитных процессов (используются данные EPDL97, EEDL97) и модель Penelope, используемые в библиотеке классов Geant4, дают несколько различный выход тормозных фотонов в области 60-90 МэВ. Для дальнейших расчетов использовались спектры, смоделированные при использовании модели низких энергий.
Электронный научный журнал "Исследовано в России"а 594 двух способах моделирования
Тип |
Модель 1 ("LE") |
Модель 2 ("Penelope") |
У |
G4LowEnergyRayleigh G4LowEnergyPhotoEl ectri с G4LowEnergyCompton G4LowEnergyGammaConversion |
G4PenelopeRayleighModel G4PenelopePhotoElectricModel G4PenelopeComptonModel G4Penel opeGammaC onver si onModel |
е- |
G4eMultipleScattering G4LowEnergyIonisation G4LowEnergyBremsstrahlung |
G4UrbanMscModel93 G4PenelopeIonisationModel G4PenelopeBremsstrahlungModel |
е+ |
G4eMultipleScattering G4LowEnergyIonisation G4LowEnergyBremsstrahlung G4eplusAnnihilation |
G4UrbanMscModel93 G4PenelopeIonisationModel G4PenelopeBremsstrahlungModel G4PenelopeAnnihilati onModel |
110000-100000-90000 -80000Х 70000 -60000 -50000 -40000 -30000Х 20000 -10000-0- |
15 20 25 30 35 40 45 50 55 60 65 70 75 80 85 90 95
I ' I |
I ' I |
тЧIЧIЧIЧr- |
1ЧIЧIЧIЧr |
Frequency Counts of A| |
I ' I ' I ' I
Ее = 90 MeV Nevents = 5661660 |
0,05
Geant 4 9.4 Penelope |
- 0,04
- 0,03
0,02
0,01
Чi i | i | i | i | i | i | i | i | i | i | i | i | i [Ч
15 20 25 30 35 40 45 50 55 60 65 70 75 80 85 90 95
Bin Center
Рис. 4. Смоделированный спектр гамма-квантов, прошедших через сборку, использовалась модель Penelope
в Geant4
Расчет сечений образования изотопаВеаа из различных мишеней производился по формуле:
<7
А -А
0аа т
O0-m-NAV-(\-e~At^)-\0~
24
(1)
где а - сечение, (барн), Фо- плотность потока у-квантов, (1/см ), Ао - активность мишени (Бк), Ат - атомная масса изотопа мишени, Nav - число Авогадро, m - масса изотопа в мишени, (г), X - постоянная распада Be .
Измерение активности каждой мишени после облучения проводилось на спектрометрическом комплексе CANBERRA InSpector 2000 с энергетическим разрешением по линии 1332 кэВ не ниже 1,74 кэВ и относительной погрешностью определения активности не более 6 %.
Электронный научный журнал "Исследовано в России"а 595 Результаты расчета средних значений (от порога до 90МеВ) сечений образования Be на различных мишенях представлены в таблице 2.
Таблица 2. Сечение образования Be.
Ядро -мишень |
масса, а.е. |
Масса мишени, г |
Поток гамма 1/см2*с 1мкА |
Активность ,Бк |
Сечение (барн) |
С |
12 |
0,91 |
8.77Е+11 |
1.827Е+04 |
2,582Е-04 |
N (из AIN) |
14 |
0,3 |
1.00Е+11 |
1.315Е+03 |
5.769Е-04 |
О |
16 |
0,48 |
5,60Е+10 |
6J80E+02 |
3.794Е-04 |
B+N |
0,561 |
3,140Е+03 |
|||
В (из BN) |
11 |
0,247 |
8J7E+11 |
1J64E+03 |
8,419Е-05 |
N (из BN) |
14 |
0,314 |
1.00Е+11 |
1.376Е+03 |
5.768Е-04 |
В литературе содержится весьма ограниченное количество экспериментальных данных по фоторождению Be в различных мишенях, но полученное нами усредненное значение образования Be на углероде достаточно хорошо согласуется с результатами , полученными в работе [5]. Следует отметить, что расчетное значение сечения образования Be значительно зависит от правильности выбора порога фотоядерной реакции, который определяется не только количеством выбитых из ядра нуклонов, но и типом вылетающих частиц : n, d, а и др. Поскольку в данном эксперименте невозможно было зафиксировать тип вылетающих из ядра-мишени частиц, то пороги фотоядерных реакций на азоте и кислороде рассчитывались по числу выбитых нуклонов, и составляют 56 и 72 МэВ, соответственно. При расчете сечения реакции бС (у, па)4Ве на углероде значение порога, согласно [5], принималось равным 25 МэВ. В дальнейшем, могут потребоваться дополнительные исследования для уточнения типа канала реакций 7N(y, Х)4Ве и 80(у, Х)4Ве , что в свою очередь, приведет к уточнению вычисленного сечения фоторождения Be на указанных ядрах.
Моделирование потока гамма квантов в атмосфере
Для проведения численных расчетов плотности потока гамма квантов в требуемом энергетическом диапазоне,а т.е.а от десяти доа сотни МэВ,аа использована программный
Электронный научный журнал "Исследовано в России"а 596
100 100
40 60
Е (MeV)
Рис. 5. Восходящая и нисходящая компоненты потока гамма-квантов на широте 40 и на высоте 5 км (а) и
10 км (Ъ) при умеренной солнечной активности.
При расчете каскада в стандартной модели атмосферы распределение по энергии первичного потока галактических протонов выбиралось в виде [8]:
Электронный научный журнал "Исследовано в России"а 597 Е (Е +2т с1)(Е+х + Ф)
2 |
(?+Ф)(?+2т/+Ф)
(2)
гдеа Еаа - кинетическая энергия протона (МэВ),а тсаа - энергия покоя протона (938 МэВ),
Ф - так называемый модулирующий множитель, учитывающий солнечную активность, который изменяется в диапазоне от 400 МэВ при минимальной солнечной активности до 1000 МэВ при максимальной.а Эмпирические константыа СДаа иаа Xаа равны следующим
величинамаа Ср = 1.244-106 1/(см2 сек МэВ), х = 780ехр(-0.00025?^) .
Полный поток гамма квантов разбивался на нисходящую (Downward) и восходящую (Upward) компоненты. Каждая из компонент имеет свою зависимость от глубины атмосферы и энергетическое распределение. Некоторые результаты моделирования потока гамма-квантов при умеренной солнечной активности показаны на рис.5. На рис.6 представлено сравнение полученных нами результатов с данными экспериментальных измерений и компьютерного моделирования, проведенного в работе [9]. Как видно программный модуль Planetocosmic пакета библиотек Geant4 дает удовлетворительное согласие с результатами измерений.
10: |
О СО N Е о |
1(Г |
10" |
10"! |
ю-
: |
ж ч |
|||
-аа |_ |
1 |
L+tf |
||
I |
тг>. |
|||
: |
1 Iа 1 |
|||
Дп( п |
||||
=- |
1 1 1аа 1 |
г |
||
1 |
||||
100 |
10
Е, MeV
Рис. 6. Сравнение результатов моделирования потока гамма-квантов на высоте 5 км (сплошная линия) с данными экспериментальных измерений и компьютерного моделирования (пунктирная линия), проведенного
в работе [9].
Электронный научный журнал "Исследовано в России"а 598 Фотоядерное производство Be в атмосфере
Расчет скорости образования Pt(h) изотопа Be в z'-ом фотоядерном канале на заданной глубине атмосферы h можно провести по следующей формуле:
Pt(h) = Nt(A)/dEai(E)J(E,h) Nt(A) < ог > /dEJ(E,h),(3)
где, iV. (И) - число атомова мишени в одном грамме воздуха на высотеа /г,аа а. (Е) - сечение
соответствующей фотоядерной реакции,а <а.>- среднее значение сечения в интервале
(Е1,Е2), J(E,H) - плотность потока фотонов. Если плотность потока фотонов задается как
[J(E,h) ] = 1/(см2 сек МэВ), то формула (3) даст число изотопов Be7, рождаемых в 1 сек в
1 грамме воздуха на заданной высоте, т.е. \_Pt(h) ] = 1/(г сек).
Интервала интегрированияа выбираетсяа ота порогаа фотообразованияаа Е1 = ЕЩ~таа в
соответствующем канале реакцииаа до максимальногоаа значения Е2. Для канала
80 (у,Х)4Ве интервал интегрирования может быть распространен вплоть до энергии 1 ГэВ если принять во внимание результаты работы [6]. Отсутствие аналогичных данных для канала 7N (у,Х)4Ве не позволяет для него распространить интегрирование в высокоэнергетичную область. Однако можно сделать верхнюю оценку, экстраполируя по аналогии с реакцией gOа (у,Х)4Ве , среднее значение сечения до энергий 1 ГэВ.
На рис. 7 представлены результаты численных расчетов скорости фотопроизводства изотопов Be в зависимости от высоты на экваторе с учетом суммарного вклада от восходящей и нисходящей компонент потока гамма квантов. Максимум производства образования Be достигается на высоте 15-17 км. На рис. 7 приведены также результаты расчетов [10] производства изотопов Be за счет реакций скалывания на протонах и нейтронах. Как видно из рисунков, вклад фотоядерного механизма в производство Be в экваториальной области земной атмосферы не мал и сравним с вкладом от протонного и нейтронного механизмов.
Электронный научный журнал "Исследовано в России"а 599 Equatorialа Be production
10а
'га
10
- 10s
10"
ЧжЧ Photoproduction
--ж-- Proton+Neutron Production
10-
10
15аа 20
h, km
25
30
35
40
Рис. 7. Скорость фотопроизводства изотопов Be7 в земной атмосфере в зависимости от высоты на экваторе (красная линия). Синяя линия - вклад от протонного и нейтронного механизма [10]. Пунктиром показана прогнозная оценка при экстраполяции данных в азотном канале до энергии гамма квантов 1 ГэВ.
На рис. 8 представлены результаты численных расчетов широтной зависимости скорости фотопроизводства изотопов Be также с учетом суммарного вклада от восходящей и нисходящей компонент потока гамма квантов. Как следует из приведенных результатов, фотоядерное производство Be в отличие от протонного и нейтронного механизмов [10] обнаруживает довольно слабую широтную зависимость.
40 Latitude
Рис. 8. Широтная зависимость скорости фотопроизводства изотопов Be7
Электронный научный журнал "Исследовано в России"а 600 В работе на пучке тормозного излучения с верхней границей энергиии 90 МэВ линейного ускорителя электронов ННЦ ХФТИ проведены исследования с целью оценки средних значений сечений (от порога до 90 МэВ) реакций А(у,Х)Ве на элементах В, С, N и О природного изотопного состава. На основании существующих данных и измеренных сечений фотообразования Be , а также моделирования ядерно-электромагнитного каскада в атмосфере впервые исследован вклад фотоядерного механизма в производство космогенного радиоизотопа Be в земной атмосфере и его широтная зависимость. Показано, что вклад фотоядерного механизма сравним с вкладом протонного и нейтронного механизмов (см., например, [4,8,10]) образования Be в атмосфере. Отсутствие данных для азотного фотоядерного канала в области высоких энергий гамма-квантов не позволяет получить более детальные оценки. Поэтому полученные в работе результаты следует рассматривать как нижнюю границу фотопроизводства Be в атмосфере. Фотоядерное производство Be в отличие от протонного и нейтронного механизмов обнаруживает довольно слабую широтную зависимость. Таким образом, проведенное исследование показывает, что вклад фотоядерного механизма в производство космогенного Be необходимо учитывать при анализе процессов накопления и переноса Be в приземном слое атмосферы. итература
- Бураева Е.А., Давыдов М.Г., Зорина Л.В., Малышевский B.C., Стасов В.В. Содержание космогенного Ве-7 в приземном слое воздуха умеренных широт.// Атомная энергия, 2007, Т. 102, В.6, С. 370-374.
- Бураева Е.А., Давыдов М.Г., Зорина Л.В., Малышевский B.C., Стасов В.В. Содержание Ве-7 в приземном слое воздуха г. Ростова-на-Дону.// АНРИ, 2007, № 1, С.63-67.
- Петрова Т.Б., Охрименко С.Е., Власов В.К., Микляев П.С. Содержание бериллия-7 в атмосферном воздухе г. Москвы. //АНРИ, 2003, № 3, С.22-29.
- Yoshimori М. Production and behavior of beryllium-7 radionuclide in the upper atmosphere.//Advances in Space Research, 2005, V.36, P.922-926.
- Кириченко В.В., Ходячих А.Ф., Вацет П.И., Догюст И.В., Золенко В.А. Исследование реакций 12C(y,pa)7Li и 12C(y,na)7Be при Еутах=120 МэВ. //Ядерная физика, 1979, Т.29, В.З, С. 572-581.
- Matsumura Н., Washiyama К., Haba Н., Miyamoto Y., Oura Y., Sakamoto К., Shibata S., Furukawa M., Fujiwara I., Nagai H., Kobayashi Т., Kobayashi K. Target-dependence of light
Электронный научный журнал "Исследовано в России"а 601
Usoskin I.G., Kovaltsov G.A. Production of cosmogenic 7Be isotope in the atmosphere: Full
Все научные статьи