Все научные статьи

Безуглов М.В., Малышевский В.С., Малыхина Т.В., Торговкин А.В., Фомин Г.В., Шраменко Б.И. Фоторождение космогенного бериллия-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 в приземном слое атмосферы. итература

  1. Бураева Е.А., Давыдов М.Г., Зорина Л.В., Малышевский B.C., Стасов В.В. Содержание космогенного Ве-7 в приземном слое воздуха умеренных широт.// Атомная энергия, 2007, Т. 102, В.6, С. 370-374.
  2. Бураева Е.А., Давыдов М.Г., Зорина Л.В., Малышевский B.C., Стасов В.В. Содержание Ве-7 в приземном слое воздуха г. Ростова-на-Дону.// АНРИ, 2007, № 1, С.63-67.
  3. Петрова Т.Б., Охрименко С.Е., Власов В.К., Микляев П.С. Содержание бериллия-7 в атмосферном воздухе г. Москвы. //АНРИ, 2003, № 3, С.22-29.
  4. Yoshimori М. Production and behavior of beryllium-7 radionuclide in the upper atmosphere.//Advances in Space Research, 2005, V.36, P.922-926.
  5. Кириченко В.В., Ходячих А.Ф., Вацет П.И., Догюст И.В., Золенко В.А. Исследование реакций 12C(y,pa)7Li и 12C(y,na)7Be при Еутах=120 МэВ. //Ядерная физика, 1979, Т.29, В.З, С. 572-581.
  6. Matsumura Н., Washiyama К., Haba Н., Miyamoto Y., Oura Y., Sakamoto К., Shibata S., Furukawa M., Fujiwara I., Nagai H., Kobayashi Т., Kobayashi K. Target-dependence of light

Электронный научный журнал "Исследовано в России"а 601

  1. Desorgher L., Fluckiger E.O., Gurtner M., Moser Ml, Butikofer R. ATMOCOSMICS: a GEANT4 code for computing the interaction of cosmic rays with the Earth's atmosphere. //International J. Modern Phys. A, 2005, V.20, P.6802-6804.
  2. Masarik J., Beer J. Simulation of particle fluxes and cosmogenic nuclide production in the Earth's atmosphere. //J. Geophys. Res., 1999, V.104, P.12099-12111.
  3. Cecchini S., Chiarusi Т., Pagliarin A., Giovannini G., On the altitude dependence of gamma-rays spectra in the Earth's atmosphere. //29th International Cosmic Ray Conference Pune (2005).
Usoskin I.G., Kovaltsov G.A. Production of cosmogenic 7Be isotope in the atmosphere: Full      Все научные статьи