Моделирование дозового распределения в программа GATE. 1.ТЕОРЕТИЧЕСКОЕ ВВЕДЕНИЕ(Гамма-излучение,фотоэффект,образование пар,КОМПТОНОВСКОЕ РАССЕЯНИЕ ).2.МЕТОД МОНТЕ-КАРЛО,АЛГОРИТМ РАБОТЫ,АНАЛИЗ ПОЛУЧЕННЫХ РЕЗУЛЬТАТОВ,ИЗМЕНЕНИЕ МАТЕРИАЛА ФАНТОМА,СВИНЕЦ,ИСТОЧНИКИ С ДРУГОЙ ЭНЕРГИЕЙ,3 ИСТОЧНИКА С РАЗЛИЧНЫМИ ЭНЕРГИЯМИ,ВЫВОД

Содержание

Введение
1⠄Глава: Теоретические основы взаимодействия гамма-излучения с веществом и метод Монте-Карло
1⠄1⠄Физические процессы взаимодействия гамма-излучения: фотоэффект, комптоновское рассеяние и образование пар
1⠄2⠄Основные характеристики и параметры гамма-излучения: энергия, интенсивность, линейный коэффициент ослабления
1⠄3⠄Метод Монте-Карло: принципы, алгоритм работы и его применение в задачах радиационной физики
2⠄Глава: Моделирование дозового распределения в $$$$$$$$$ $$$$
2⠄1⠄$$$$$$$$$$ $$$$$$ и $$$$$$$$ $$$$$$$$$ $$$$$$$: $$$$$$$$$ $$$$$$$$$ и $$$$$$$$$$ ($$$$, $$$$$$)
2⠄2⠄Моделирование дозового распределения $$$ $$$$$$$$$$$$$$$$$$ $$$$$$$$$$: $$$$$$ $$$$$$$ $$$$$$$ излучения ($$$ $$$, $$$ $$$, 1.$$ $$$)
2⠄3⠄Моделирование дозового распределения $$$ $$$$ $$$$$$$$$$ с $$$$$$$$$$ $$$$$$$$$ и $$$$$$$$$ $$$$$$$$$ $$$$$$$: $$$$$$ $$$$$$$$$$ $$$$$$$$$$$ и $$$$$$$$$
$$$$$$$$$$
$$$$$$ $$$$$$$$$$$$$$ $$$$$$$$$$

Введение

В современной ядерной медицине, радиотерапии и радиационной безопасности ключевое значение имеет точное прогнозирование распределения поглощенных доз ионизирующего излучения в биологических тканях и защитных материалах. Развитие методов лучевой терапии и диагностики требует не только высокой эффективности воздействия на патологические очаги, но и минимизации рисков для здоровых органов. В этой связи численное моделирование процессов переноса излучения, в частности с использованием метода Монте-Карло, становится незаменимым инструментом. Программа GATE (Geant4 Application for Tomographic Emission), основанная на библиотеке Geant4, предоставляет широкие возможности для детального моделирования взаимодействия гамма-излучения с веществом, позволяя с высокой точностью рассчитывать дозовые распределения в фантомах различной геометрии и состава. Актуальность данной работы обусловлена необходимостью верификации расчетных методов для задач радиационной физики и оптимизации параметров облучения в практических приложениях.

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

Объектом исследования является процесс переноса гамма-излучения в веществе и формирование дозового поля. Предметом исследования выступает дозовое распределение, $$$$$$$$$$ в $$$$$$$$$$ $$$$$$$$$$ $$$$$$$$$$$$$ в $$$$$$$$$ $$$$ $$$ $$$$$$$$$ $$$$$$$$$$ $$$$$$$ ($$$$, $$$$$$) и $$$ $$$$$$$$$ $$$$$$$$$$$$$$ $$$$$$$$$$$$$$$ $$$$$$$$$$ гамма-излучения.

$$$$$ $$$$$$ $$$$$$$$ $$$$$$ $$$$$$$$ $$$$$$$$$$$$$ $$$$$$$$ $$$$$$$$$$$$$ $$$$$-$$$$$$$$$ $ $$$$$$$$$ $$$$ $ $$$$$$ $$$$$$$ $$$$$$$$$ $$$$$$$ $ $$$$$$$ $$$$$$$$$ $$ $$$$$$$$$$$$$$ $$$$$$$$$$$ $$$$.

$$$ $$$$$$$$$$ $$$$$$$$$$$$ $$$$ $$$$$$$$$$ $$$$$$ $$$$$$$$$ $$$$$$:
$. $$$$$$$ $$$$$$$$$$$$$ $$$$$$ $$$$$$$$$$$$$$ $$$$$-$$$$$$$$$ $ $$$$$$$$$ ($$$$$$$$$$, $$$$$$$$$$$$$ $$$$$$$$$, $$$$$$$$$$$ $$$) $ $$$$$$$$ $$$$$$ $$$$$-$$$$$.
$. $$$$$$$ $$$$$$$ $$$$$$$$ $$$$$$ $ $$$$$$$$$ $$$$, $$$$$$$ $$$$$$$$ $$$$$$$$$ $$$$$$$, $$$$$$$ $$$$$$$$$$ $ $$$$$$$$$ $$$$$$$$$$ $$$$$$$$$$$$$.
$. $$$$$$$$ $$$$$ $$$$$$$$$$$$$$ $$$$$$$$$$$$$ $$$ $$$$$$$$$$ $$$$$$$ ($$$$ $ $$$$$$) $$$ $$$$$$$$$ $$$$$$$$ $$$$$-$$$$$$$$$.
$. $$$$$$$$$ $$$$$$ $$$$$$$$$$ $$$$$$$ $$$$$$$$$$$$$ $ $$$$$$$$ $$$$$$$$$$ $$$ $$$$$$ $$$$$$$ $$$$$$$$$$$$$.
$. $$$$$$$$$$$ $$$$$$$ $$$$$$$$$$$$$ $$$ $$$$$$$ $$$$ $$$$$$$$$$ $ $$$$$$$$$$ $$$$$$$$$ $ $$$$$$$ $$$$$$ $ $$$$$$$$$ $$$$$$$$$$ $$$$$$$$$$$.

$$$$$$ $$$$$$$$$$$$ $$$$$$$$ $$$$$$$$$$$$$ $$$$$$ $$$$$$-$$$$$$$$$$$ $$$$$$$$$$, $$$$$ $$$$$$$$$$$$$$$ $$$$$$$$$$$$$ ($$$$$ $$$$$-$$$$$), $$$$$$$$$$$$$ $$$$$$ $$$$$$$$$$ $$$$$$, $ $$$$$ $$$$$$ $$$$$$$$$$$$ $ $$$$$$$$$$$$$$ $$$$$$$$$ $$$$$$$$$$$ $$$$$$$$$$$$$.

$$$ $$$$$$$$$ $$$$$$ $$$$$ $$$$$$$$$$$$ $$$$$$$$$$$ $$$$$$$ $ $$$$$$$ $$$$$$$$$, $$$$$$$ $$$$$$$$$$ $$ $$$$$$ $$$$$$$$$$$$$$ $$$$$$$$$ $ $$$$$$$$$, $$$$$$ $$ $$$$$$$$$$$$$ $$$$$$$$ ($$$$$$$$, $$$$$$$ $$ $$$$$$$$ $$$ $$$$$$$, $$$$$$$ $$$$$$$), $ $$$$$ $$$$$$$$$$$ $$$$$$$$$$$$ $ $$$$$$$$$$$ $$$$$$$$$$$$ $$ $$$$$$$$$ $$$$ $ $$$$$$$$$$ $$$$$$.

Физические процессы взаимодействия гамма-излучения с веществом

Гамма-излучение представляет собой поток электромагнитных волн с чрезвычайно малой длиной волны (менее 10⁻¹⁰ м), испускаемых при ядерных переходах, аннигиляции частиц или торможении заряженных частиц высоких энергий. В отличие от альфа- и бета-частиц, гамма-кванты не имеют массы покоя и электрического заряда, что определяет их высокую проникающую способность и специфический характер взаимодействия с веществом. Как отмечается в современных исследованиях, понимание механизмов ослабления гамма-излучения имеет критическое значение для задач радиационной защиты, ядерной медицины и дефектоскопии [12]. При прохождении через среду гамма-кванты могут взаимодействовать с атомами вещества тремя основными способами: фотоэлектрическим поглощением, комптоновским рассеянием и образованием электрон-позитронных пар. Вероятность реализации каждого из этих процессов определяется энергией фотона и атомным номером материала.

Фотоэлектрический эффект (фотоэффект) представляет собой процесс, при котором гамма-квант полностью поглощается атомом, передавая свою энергию одному из связанных электронов внутренних оболочек (K, L, M). Выбитый электрон, называемый фотоэлектроном, покидает атом с кинетической энергией, равной разности энергии падающего фотона и энергии связи электрона на данной оболочке. Освободившееся место на внутренней оболочке заполняется электроном с более высоких уровней, что сопровождается испусканием характеристического рентгеновского излучения или электронов Оже. Вероятность фотоэффекта резко возрастает с увеличением атомного номера Z материала (пропорционально Z⁴–Z⁵) и быстро убывает с ростом энергии гамма-кванта (пропорционально E⁻³·⁵). Данный процесс доминирует при низких энергиях (менее 100 кэВ) и в материалах с высоким Z, таких как свинец, что широко используется при конструировании защитных экранов.

Комптоновское рассеяние является преобладающим механизмом взаимодействия для гамма-излучения с энергиями от 0.1 до нескольких мегаэлектронвольт в материалах с низким и средним атомным номером. В этом процессе гамма-квант упруго рассеивается на практически свободном (или слабо связанном) электроне внешних оболочек атома. При столкновении фотон передает часть своей энергии электрону и изменяет направление своего движения. Энергия рассеянного фотона зависит от угла рассеяния: максимальная передача энергии происходит при рассеянии назад (угол 180°), минимальная — при малых углах. Рассеянный фотон может затем участвовать в дальнейших взаимодействиях, что существенно усложняет картину распространения излучения в среде. Как указывается в $$$$$$$$$$, $$$$ $$$$$$$$$$$$$ $$$$$$$$$$$$$$ рассеяния является $$$$$ $$ $$$$$$$$ $$$$$ при $$$$$$$$$$$$$ $$$$$$$ $$$$$$$$$$$$$ в $$$$$$$$$$$$$ $$$$$$ [$$]. $$$$$$$ $$$$$$$$$$$$$$ рассеяния $$$$$$$$$$$$$$$ $$$$$$$$ $$$$$$ $ и слабо зависит от энергии фотона в $$$$$$$ $$$$$$$$$.

$$$$$$ $$$$$$$$$$$ $$$$$$$$-$$$$$$$$$$$ $$$ $$$$$$$$$$ $$$$$$$$$$$$$ $$$$$$$$$ $$$ $$$$$$$ $$$$$-$$$$$$, $$$$$$$$$$$ $.$$$ $$$ ($$$$$$$$$ $$$$$$$ $$$$$ $$$$$$$$$). $ $$$$ $$$$$$$$ $ $$$$$$$$$$$ $$$$ $$$$ $$$ $$$$$$$$$ $$$$$-$$$$$ $$$$$$$$$$$$ $ $$$$ $$$$$$: $$$$$$$$ $ $$$$$$$$. $$$$$$$ $$$$$$$ $$$$$$ $$$$$ $$$$$$$$$$ $$$$$$$$ $$$$$$$$$ $ $$$$$$$$$$$$ $$$$$$$ $$$$$$$$$$$$$$ $$$$$$. $$$$$$$$, $$$$$$ $$$$$$$$$ $$$$$$$$$$ $ $$$$$, $$$$$$$$$$$$ $ $$$$$ $$ $$$$$$$$$$, $$$$$$$$ $$$ $$$$$$$$$$$$$$$ $$$$$-$$$$$$ $ $$$$$$$$ $$$ $$$ $$$$$$. $$$$$$$$$$$ $$$$$$$$$$$ $$$ $$$$$$$$$$ $ $$$$$$$$$$$ $$$$$$$ $$$$$$ $ $$$$$$$$$$$$$$$ $$$$$$$$ $$$$$$$$ $$$$$$ $$$$$$$$$ ($$). $$$$$$ $$$$$$$ $$$$$$$$$$ $$$$$$$$$$$$ $$$ $$$$$$$ $$$$$$$$ ($$$$$ $–$$ $$$), $$$$$$$$ $ $$$$$$$ $$$$$$$$$$.

$$$$$$$$$$ $$$$$$$$ $$$$ $$$$ $$$$$$$$$ $$$$$$$$$$$ $$$$$$$$ $$$$$$$$$$$$$ $$$$$$$$$$ $, $$$$$$$ $$$$$$$$$$$$ $$$$$ $$$$$ $$$$$$$$$$$ $$$$$$$$$$$$$ $$$ $$$$$$$$$$$ $, $$$$$$$$$$$$$$ $$$$$$$$$ $ $ $$$$$$$$$$$ $$$ $. $$$$$ $$$$$$$$$$ $$$$$$ $$$$$ $$$$$-$$$$$$$$$ $$$$$ $$$$$$$$$$$$$$$$ $$$$$$$$: $ = $$·$$$(-$$), $$$ $$ — $$$$$$$$$ $$$$$$$$$$$$$, $ — $$$$$$$ $$$$ $$$$$$$$. $$$$$$ $ $$$$$$$$ $$$$$$$$, $$$$$$$$ $$$ $$$$$$$ $$$$$$ $ $$$$$$$ $$$$$$$$$ $$$$$$$$$, $$$$$$$$$$ $$$$$$$$ $$$$$$$$$$$$$$$$$ $$$$$$ $$$$$$$$$$$ $$$$$$$$$$$$$ $$-$$ $$$$$$ $$$$$$$$$$$ $$$$$$$$$. $$$ $$$$$$$ $$$$$$$ $$$$$$$ $$$$$$$$$$$$$ $$$$$$$$$$ $$$$$$$$$ $$$$$$$$$$$$ $$$$$$$$$$$$$$ $ $$$$$$$$$ $$$$$$$, $$$ $$$$$$$$ $$$$$$ $ $$$$$$$$$$$$$$ $$$$$$$$$ $$$$$$$, $$$$$ $$$ $$$$$ $$$$$-$$$$$ [$$].

$$$$$ $$$$$$$$, $$$ $$$$$$$$$$$$$ $$$$$ $$$$$$$ $$ $$$$$$$$$$$$$ $$$$$$$$$ $$$$$$$$$$$ $$$$$$$ $$ $$$$$$$ $$$$$-$$$$$$$$$ $ $$$$$$$ $$$$$$$$$$ $$$$$. $$$$$$$$, $ $$$$$$$$$$$$$ $$$$$$ ($$$$, $$$$$$ $$$$$) $$$ $$$$$$$$, $$$$$$$$ $$$ $$$$$$$$$$$$ ($–$$ $$$), $$$$$$$$$$ $$$$$$$$$$$$$ $$$$$$$$$, $$$$$ $$$ $ $$$$$$$ $$$$$ ($ ≈ $$) $ $$$$$$$$ $ $$$$$$ ($ = $$) $$$$$$$$$$ $$$$ $$$$$$$$$$$. $$$ $$$$$$$$$$$$$ $ $$$$$$$$$ $$$$ $$$ $$$ $$$$$$$$ $$$$$$$$$$$ $$$$$$$$$$$$$ $$ $$$$$$ $$$$$$$$$$ $$$$$$, $$$$$$$$$$$$ $$$$$$$$$$ $$$$ $$$$$$ $$$$$$$ $$$$$$$$$$$$$$. $$$$$$$$ $$$$$$$$$ $$$$$$ $$$$$$$$$$$$$$ $$$$$-$$$$$$$$$ $ $$$$$$$$$ $$$$$$$$$$ $$$ $$$$$$$$$$ $$$$$$$$$$$$$ $$$$$$$$$$$ $$$$$$$$$$$$$ $ $$$$$$$$$$$$$ $$$$$$ $$$$$$$$$$ $$$$$$$$$$$$$$$ $$$$$$$$$$$$.

Принципиальное значение для понимания процессов формирования дозовых распределений имеет анализ энергетической зависимости сечений взаимодействия. Для фотоэлектрического поглощения характерны резкие скачки сечения при энергиях, соответствующих энергии связи электронов на K-, L-, M-оболочках (K-край поглощения). При энергии фотона, незначительно превышающей энергию связи K-электрона, сечение фотоэффекта резко возрастает, после чего монотонно убывает с ростом энергии. Это явление имеет важное практическое значение: например, при выборе материалов для рентгеновских фильтров или при интерпретации результатов спектрометрии гамма-излучения. В контексте моделирования в GATE корректный учет K-края необходим для точного воспроизведения выхода характеристического излучения, которое может вносить заметный вклад в дозу на небольших расстояниях от точки взаимодействия.

Комптоновское рассеяние характеризуется дифференциальным сечением, зависящим от угла рассеяния. Распределение рассеянных фотонов по углам описывается формулой Клейна-Нишины, которая учитывает поляризацию излучения. Для изотропного падающего пучка наибольшая вероятность рассеяния соответствует малым углам (вперед) при высоких энергиях фотона, тогда как при низких энергиях распределение становится более симметричным. Энергия рассеянного фотона E' связана с энергией падающего фотона E и углом рассеяния θ соотношением: E' = E / [1 + (E/m₀c²)(1 - cosθ)], где m₀c² — энергия покоя электрона (511 кэВ). Из этой формулы следует, что при рассеянии на 180° (обратное рассеяние) энергия фотона не может превысить 255 кэВ, независимо от начальной энергии. Данный эффект объясняет появление характерного пика обратного рассеяния в спектрах гамма-излучения, регистрируемых детекторами.

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

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

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

Альтернативным процессом заполнения вакансии является оже-эффект, при котором энергия перехода передается другому электрону внешних оболочек, который вылетает из атома (оже-электрон). Вероятность оже-эффекта выше для легких элементов и для переходов с участием внешних оболочек. В программе GATE оба процесса — флуоресценция и оже-эффект — моделируются с использованием данных из библиотеки EADL (Evaluated Atomic Data Library), $$$ $$$$$$$$$$$$ $$$$$$$ $$$$$$$$ $$$$$$$$$$$$$$$ $$$$$$$$$ $$$$$$$$$ [$$].

$$$ $$$$$$$$$$$$ $$$$$$$$$$$$$$ $$$$$-$$$$$$$$$ $ $$$$$$$$$ $$$$$$$$$$ $$$$$ $$$$$$$$$ $$$$$$$, $$$$$$$$$ $ $$$$$$$$$$$$ $$$$$$$$$$. $ $$$$$$$ $$$$$$$ $$$$$ $$$$$ $$$$$$$$ $$$$$$$$$ $$$$$$$$$$$$$$$$ $$$$$ $$$$$$$$$$$$$$ $$$$$$$$$, $$$$$$ $$$ $$$$$ $$$$$$$$$ $$$$$$$$ $$$ $$$$$$$ $$$$$. $$$$$$ $$$ $$$$$$$$$ $$$$$$$$$ $$$$$$$ $$$$$$ $ $$$$$$$$ $$$$$$$$$$$ $$$ $$$$$$$$. $ $$$$$$$$$$ $$$$$$$$$$$ $$$$ $$$$$$$$$$$ $$$$$$$$$, $$$$$$$ $$$$$ $$$$$$$$$$$ $$$$$$$$$$ $$ $$$$$$$$$$ $$$$$ $$ $$$$$$$$ $ $$$$$$$$$$$$$$$ $$$$$$$$$$$$$. $$$$$$ $$$$$$$ $$$ $$$$$$$ $$$$$$$ $$$$$$$ $$$$$$$$$$$$$ $ $$$$$$$$$$$ $$$$$$$$, $$$$$ $$$ $$$$$$ $$$$$$$ $$$ $$$$$$$$$$$$$ $$$$$, $$$$$$$$$$ $$$$$$$$$$$$ $$$$$$$$$ $$$$$$, $$$$$$$$$$$ $$$$$$$$$$$$ $$$$$$$$$$ $$$$$$$$ $$$$$ $$$$$$$ $ $$$$$$ $$$$ $$$$$$$$$ $$$$$$$$$$$$$$.

$$$$$$$ $$$$$$$$ $$$$$$$$$$$ $$$$$$ $ $$$$$$$$ $$$$$$$ $$$$$. $$$ $$$$$$$$$$$ $$$ $$$$$$$ $$$$$$ ($$ $$$$$$$ $$$$$$$ $$$$$ $$$$$$$$$) $$$$$$$$$$ $$$$$$$$$$$$$, $$$$$$$ $$$$$ $$$$$$$$$$ $ $$$$$$$$, $$$$$$$$ $$$$$$$$$ $$$$$$$ $$$$$$$$$$$$$. $$$ $$$$$$$$$$$$$ $$$$$$$$$ $$$$$ $$$$$$$ $$$$$$$$ $$$$$$$$$$ $$$$$$$, $ $$$$$ $$$$$$$$$$ $$$$$$$$$ $$$$$$. $$$ $$$$$$$$$$$ $$$ $$$$$$$ $$$$$$$$$$$$$$ $$$$$ $$$$$$$$$$, $$$$$$$$$$ $, $$$$$$$$, $$$$$$$$$$$$$$$$ $$$$$$$$. $$$$$ $$$$$$$, $$$$ $ $$$$$$ $$$$$ $$$$$ $$$$$$$$$$$$ $$$$$$ $$$$$$$ $$ $$$$ $$$$$$$$$$ $$$$$$ ($$$$$$$$$$, $$$$$$$$$$), $$$$$$$$$$$$$$ $ $$$$$$$$$$ $$$$$$$$$$$$$$ $$$$$-$$$$$$$ $ $$$$ $$$$$ $$$ $$$$$$ $$$. $$$$$$$ $$$$$$$$$ $$$$$$$$$$ $ $$$$$$$$$$$$$ $$$$$$ $$$$$$$$$$ $$ $$$$$ $$$$$$$$$$ ($$$ $$$$$$ $$$$$$$$) $$ $$$$$$$$$$ $$$$$$$$$$$ ($$$ $$$$$$$$ $$$$$$$ $$ $$$), $$$ $$$$$$$$$$ $$$$$$$$$$$$$$$$ $$$$$$$$$$ $$$$$$$ $$$$$$$$$$$$$.

$ $$$$$$$$$ $$$$$$$$$$$$$ $ $$$$$$$$$ $$$$ $$$ $$$$$$$$$ $$$$$$$$$$ $$$$$$$$ $$$$$$$$$$$ $ $$$$$$$$$$$$$$ $$$$$$$$$$$ $ $$$$$$$$$$$$$$$$$$$ $$$$$$$ $$$$$$$$$$ $$$$$$. $$$ $$$$$ $$$$$$$ $$$$$$$$ $ $$$$$$$$$$$$ $$$$$$, $$$ $$$$$$$ $$$$$-$$$$$$$$$ $$$$$$ $$ $$$$$$$$$ $$$$$$$$$$ $$$$$$$$$$$$$$$$$, $$$$$$$$$$$$$ $$$$$$$$$$$$ $$$$$ $$$$$$$$$$ $$$$$$$$$ "$$$$$$$$$" $$$ "$$$$$$$$", $$$$$$$ $$$$$$$$$$$$ $$$$$$$ $$$$$$$$ $$$$$$$$ $$$$$$$$$$$$$$ $ $$$$$$$$$ $$$$$$$ $$ $$$ $$ $$ $ $$$. $$$ $$$$$$ $$$$$$$$ $$$$$$$$$ $$$$$$$$ $$$$$$$$$$$ $ $$$$$$ $$$$$$$$$$$$$, $$$$$$$$$$$$$$ $$$$$$$$$ $$ $$$$$$$ $$$$$$-$$$$$$, $$$$$$$ $$$$$$$$$$$ $$$ $ $$$$$$$$$$$ $$$$$$$$$. $$$$$ $$$$$$$$$$$ $$$$$$ $$$$$$$$$$ $$$$$$$$$ $$$$$$$ $$ $$$$$$$$$ $$$$$$$$ $ $$$$$$$$$ $$$$$$$$$$$$$$ $$$$$$$$.

$$$$$ $$$$$$$, $$$$$$$$$$$$$ $$$$$$$$$$ $$$$$$$$ — $$$$$$$$$$, $$$$$$$$$$$$$ $$$$$$$$$ $ $$$$$$$$$$$ $$$ — $$$$$$$$$$ $$$$$$ $$$$$$$$$$$$$$ $$$$$-$$$$$$$$$ $ $$$$$$$$$. $$ $$$$$$$$$$$$$ $$$$$ $$$$$$$$$$$$ $$$$$$$$ $$$$$$$ $ $$$$$$$ $$$$$$$ $$$$$$$$$. $$$ $$$$$$ $$$$$$$ $ $$$$$$$ $$$$$$$$$$ $$$$$$$$$$ $$$$$$$$$$, $$$ $$$$$$$ $$$$$$$ $ $$$$$$ $$$$$$$$$$ — $$$$$$$$$$$$$ $$$$$$$$$, $$$ $$$$$$$ $$$$$$$ — $$$$$$$$$$$ $$$. $$$$$$$$$ $$$$ $$$$$$$$$$$$$$$ $$$$$$$$ $$$$$$$$$$$ $$$$$$$$ $$$ $$$$$$$$$$ $$$$$$$$$$$$$ $$$$$$$$$$$ $$$$$$$$$$$$$ $$$$$$$ $$$$$$$$$$$$$ $ $$$$$$$$$ $$$$. $$$$$$$$$ $$$$ $$$$ $$$$$$$$$ $$$$$$$$$ ($$$$$$$$$$$$$, $$$$$$$$$$$$$$$ $$$$$$$$$, $$$$$$$$$ $$$$$$$$$) $$$$$$$$$ $$$$$$$ $$$$$$$ $$$$$$$$ $$$$$$$$, $$$ $$$$$$$$ $$$$$ $$$ $$$$$ $$$$$$$$$$$$ $$$$$$$ $ $$$$$$$$$$$$ $$$$$$$$$$$$. $$$$$ $$$$$-$$$$$, $$$$$$$$$$$$$ $ $$$$, $$$$$$$$$$$$ $$$$$$$$$$$$$ $$$$$$$$$$$ $$$$$$$$ $$$$$$$$$$$$$$$ $$$$$$$$$ $$$$$$$$$$$$$$ $$$$$-$$$$$$$ $ $$$$$$$$$, $$$ $$$$$$ $$$ $$$$$$$$$$$ $$$$$$$$$$$$ $$$ $$$$$$$$$$$$ $$$$$$$ $$$$$ $ $$$$$$$ $$$$$$$$$$$$ $$$$$$ [$].

Метод Монте-Карло в задачах радиационной физики

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

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

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

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

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

$$$ $$$$$$$$$ $$$$$$$$$$$$$ $$$$$$$$$$$$$ $ $$$$$$ $$$$$-$$$$$ $$$$$$$$$$$$ $$$$$$$$$ $$$$$$ $$$$$$$$$$ $$$$$$$$$. $$$$$ $$ $$$$$$$$ $$$$$$$$$$$$$$$$ $$$$$$$$ $$$$$ "$$$$$$$$", $$$ $$$$$$$ $$$$$$$$$$$$ $$$$$$$$$$$ $$ $$$$$$$ $ $$$$$$$$$ $$$$$$$$$$$ $$$ $$$$$$$$$$$ $$$$$$$$. $ $$$$$$$$, $$$ $$$$$$$$$$$ $$$$$$$$$$$$$$ $$$$ $$$ $$$ $$$$$ $ $$$$ $$$$$$$$$$$$, $$$$$$$$$$$$ $$$$$$ $$$$$$, $ $ $$$$$$$$, $$$$$$$$$$$$$$ $$$$$$$, — $$$$$$. $$$$$$ $$$$$$$$$$$ $$$$$$$ $$$$$$$$ "$$$$$$$$$$$" ($$$$$$$$$), $$$ $$$$$$$ $$$$$$$, $$$$$$$$$$ $ $$$$$$ $$$$$$$, $$$$$$$$$$ $$$$$$$$$$$ $$$$$$$$$ $ $$$$$$$$$$$$$$$ $$$$$$$$$$$ $$$$. $$$$$ "$$$$$$$ $$$$$$$" ($$$$$$$ $$$$$$$$) $$$$$$$$$$$$ $$$ $$$$$$$$$$$ $$$$$$$$$$$$$ $$$$$$, $$$$$$$$ $ $$$$$$$$$$ $$$$$$$. $ $$$$$$$$$ $$$$ $$$ $$$$$$ $$$$$ $$$$ $$$$$$$$$$$ $ $$$$$$$$$$$$$$ $$$$$$$$$$$ $$$$$$ $ $$$$$$$$.

$$$ $$$$$$$$$$$$$ $$$$$$$ $$$$$$$$$$$$$ $ $$$$$$$$ $$$$$$$$$$ $$$$$$$ $$$$$$ $$$$ $$$$$$ $$$$$$$$$$ $$$$$ $$$$$$$$$ $$$$$$$ $$$ $$$$$$$$$ $$$$$$$$$ $$$$$$. $ $$$$$$ $$$$$$$$$$ $$$$$$$ "$$$$$$$$$$$$ $$$$$$" ($$$$$$$$$$ $$$), $$$$$$$ $$$$$$$$$$ $$$$$$$$$$$ $$$$$$$, $$$ $$$$$$$ $$$$$$$$$ $$$$$$$ ($$$$$$$$$, $$$$$$) $$$$$ $$$$ $$$$ $$$$$$$$$$$$$. $$$$$$$ $ $$$$$$$$ $$$$ $$$$$ $$$$$$ $$$$$$$$$ $$$$$$$$ $$$$$$$$$$$$, $ $$ $$$$$$$ $$$$$$$$$$$ $ $$$$ $ $$$$$ $$$$$$$$$$$$$$. $$$$$ $$$$$$ $$$$$$ $$$ $$ $$$$$$$$ $$$$$$$$$$$$$, $$$ $ $$ $$$$$$$$$$$$$$ $$$$$$$: $$$$$$$ $$$$$$ $$$$$ $$$$$$$$ $ $$$$$$$$$ $$$$$$$$$ $$$$$ $$$$$$$$$$$$$$$$$$$ $$$$$$, $$$ $$$$$ $$$$$$$$$$$ $$$$$ $$$$$$$, $ $$$$$$$ $$$$$$$ $$$$$ $$$$$ $$$$$$$$ $ $$$$$$ $$$$$$$$, $$$$$$$$ $$$ $$$$$$$$$$$$$ $$$$$$$$$$$ $ $$$$$$$$$$$$$$$ $$$$$$$$$.

$ $$$$$$$$$ $$$$$$ $$$$$$ $$$$$ $$$$$-$$$$$ $$$$$$$$$$$$ $$$ $$$$$$$$$$$$$ $$$$$$$ $$$$$$$$$$$$$, $$$$$$$$$$$ $$$$$-$$$$$$$$$$ $ $$$$$$ $ $$$$$$$$$ $$$$$$$$. $$$$ $$$$$$$ $ $$$$$$$$ $$$$$$$$$, $$$$$$$$$$$$ $$$$$$$$$$$$$ $$$$$, $$$$$$$$$ $$ $$$$$$$$$ $ $$$$$$$$$$$ $$$$$$$ $$$$$ $$$$$$ $ $$$$$$$$$$ $$$$$$ $$$$$$ $$$$$$$$. $$$$$$, $$$$$$$$$$ $$$$$$$ $$$$$$$ $$$$$$$ ($=$$) $ $$$$$$$ $$$$$$$$$$ ($$.$$ $/$$$), $$$$$$$$ $$$$$$$$$$$$ $$$$$$$$$$ $$$ $$$$$$$$$$$$ $$$$$$. $$$$$$$$$$$$$ $$$$$$$ $$$$$$$$$$$$$ $ $$$$ $$$$ $$$$$$$$$$ $$$$$$$$$ $$$$$$$$ $$$$$$$$$$$$$$$$$$ $$$$$$$ $$$$$$$$ $$$$$$ $$ $$$$$$$$ $$$$$$$$$$ $$$$$-$$$$$$$$$ $ $$$$$$$$$$$$ $$$$$$$$ $$$$.

$$$$$ $$$$$$$, $$$$$ $$$$$-$$$$$ $$$$$$$$ $$$$$$ $ $$$$$$$$$$$$$ $$$$$$$$$$$$ $$$ $$$$$$$$$$$$$ $$$$$$$$$ $$$$$$$$ $$$$$$$$$$$$$ $$$$$$$$$. $$$ $$$$$$$$$$ $ $$$$$$$$$ $$$$ $$$$$$$$$ $ $$$$$$$ $$$$$$$$$ $$$$$$$$$$$$ $$$$$$$ $$$$$$$$$$$$$ $ $$$$$$$ $$$$$$$$$$$$ $$$$$$ $ $$$$$$ $$$$ $$$$$$$$ $$$$$$$$$$ $$$$$$$$$. $$$$$$$$$$$$$$ $$$$$$$ $$$$$$ $$$$$$$ $$$$$$$$$$$$$ $$$$$$$$ $$$$$ $$$$$$ $$$ $$$$$$$$$$ $$$$$$$$$$ $$$$$$$$, $$$ $$$$$$$$$$$$$ $$$$$$$$$$$$$ $$$$$$$$$$$$$ $$$$$$$$$$$ $$$$$$$$$$ $ $$$$$$$ $$$$$$$$$$ $$$$$$$$$. $$$$$$$$$ $$$$$$$$$ $$$$$$ $$$$$$ $$$$$-$$$$$ $$$$$$$$$$ $$$ $$$$$$$$$$ $$$$$$$$$$ $$$$$$$$$$$$$$ $$$$$$$$$$$$$ $ $$$$$$$$$$$$$ $$$$$$$$$$ $$$$$$$$$$$. $ $$$$$$$$$ $$$$$$$ $$$$$ $$$$$$$$$$ $$$$$$$$$$ $$$$$$$$ $$$$$$ $ $$$$$$$$$ $$$$ $ $$$$$$$$$$$ $$$$$$$$$ $$$$$$$$$$ $$$$$$$$$$$$$ $$$ $$$$$ $$$$$$$$$$ [$$].

Важным аспектом применения метода Монте-Карло в программе GATE является корректный выбор физических моделей для описания взаимодействия излучения с веществом. Библиотека Geant4, лежащая в основе GATE, предоставляет несколько наборов физических процессов, различающихся по точности и вычислительной эффективности. Для задач низкоэнергетической физики, включая моделирование гамма-излучения с энергиями до нескольких мегаэлектронвольт, рекомендуется использовать наборы "Livermore" или "Penelope". Эти модели включают детальное описание фотоэлектрического эффекта с учетом флуоресценции и оже-электронов, комптоновского рассеяния по формуле Клейна-Нишины с учетом связывания электронов, эффекта образования пар и релеевского рассеяния. Особенностью набора "Penelope" является использование более точных алгоритмов для моделирования угловых распределений при комптоновском рассеянии, что может быть важно для задач, требующих высокой точности.

При моделировании дозовых распределений в гетерогенных средах, состоящих из нескольких материалов, необходимо учитывать особенности переноса излучения на границах раздела сред. На границе между материалами с различными атомными номерами возникает явление, известное как "эффект накопления дозы" (dose build-up effect). Этот эффект обусловлен тем, что вторичные электроны, образовавшиеся в материале с высоким Z, могут проникать в соседний материал с низким Z, создавая там дополнительный вклад в дозу. В программе GATE границы раздела сред обрабатываются автоматически на уровне геометрического движка Geant4, который определяет, в каком материале находится частица в каждый момент времени, и применяет соответствующие физические процессы.

Особого внимания заслуживает вопрос о регистрации и анализе результатов моделирования. В GATE для этих целей используются так называемые "актеры" (actors) — специальные объекты, которые выполняют определенные действия в процессе моделирования. Для расчета дозовых распределений применяется актер "DoseActor", который позволяет регистрировать поглощенную дозу в каждой ячейке заданной сетки (воксельной решетки). Результаты могут быть представлены в виде двумерных или трехмерных массивов данных, которые затем могут быть визуализированы с использованием специализированного программного обеспечения, такого как ROOT, ParaView или MATLAB. Актер "DoseActor" также позволяет рассчитывать статистические погрешности и строить профили дозы вдоль заданных направлений.

При анализе полученных результатов важную роль играет оценка статистической значимости наблюдаемых различий. Поскольку метод Монте-Карло дает оценки с определенной погрешностью, необходимо убедиться, что различия между дозовыми распределениями для разных материалов или энергий не обусловлены случайными флуктуациями. Для этого используются стандартные статистические критерии, такие как t-критерий Стьюдента или критерий хи-квадрат. В программе GATE предусмотрена возможность автоматического расчета статистических погрешностей и построения доверительных интервалов для каждой ячейки дозовой сетки.

Современные версии GATE поддерживают параллельные вычисления с использованием технологии MPI (Message Passing Interface), что позволяет существенно ускорить моделирование за счет распределения задачи на несколько процессоров или ядер. При параллельном моделировании каждая первичная частица обрабатывается независимо, что обеспечивает хорошую масштабируемость. Однако при использовании параллельных вычислений необходимо учитывать особенности генерации псевдослучайных чисел: каждый процесс должен использовать свой независимый поток случайных чисел, чтобы избежать корреляций между результатами. В GATE эта задача решается с использованием библиотеки CLHEP (Class Library for High Energy Physics), которая предоставляет генераторы случайных чисел с различными алгоритмами.

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

В контексте данной работы особый интерес представляет сравнение дозовых распределений в воде и свинце при одинаковых условиях $$$$$$$$$. $$$$, $$$ $$$$$$$$ $ $$$$$$ $$$$$$$$$$$ $$$$$$$ $$$$$$$ ($$$$$ ≈ $.$), $$$$$$$$$$$$$$$ $$$$$$$$$$$$$ $$$$$$$$$$$$$$ $$$$$$$$$ в $$$$$$$ $$$$$$$$$ $$$$$$$, $$$ $$$$$$$$ $ $$$$$ $$$$$$$$$$$$ $$$$$$$$$$$$$ $$$$ $$ $$$$$$ $$$$$$$. В свинце, $$$$$$$$, $$$$$$$$$$ и $$$$$$$$$$$ $$$ $$$$$$$$ $ $$$$$ $$$$$$$$$$ $$$$$$$$$$ $$$$$$$ $$$$$$ $$$$$ $$$$$$$$$$$$$$, $$$ $$$$$$$$$$$ в $$$$$ $$$$$$ $$$$$ $$$$$$$$ $$$$$$$$$$$$$ $ $$$$$$$$. $$$$$$$$$$$$$$ сравнение $$$$ распределений $$$$$$$$$ $$$$$$$ $$$$$$$$$$$$$ $$$$$$ $$$ $$$$$$$$$ $$$$$$$$$ и $$$$$$$ $$$$$$$$$$$ $$$$$$$$$$$$ дозовых $$$$$ в $$$$$$$$$ $$$$$$.

$$$$$ $$$$$$$$, $$$ $$$ $$$$$$$$$$$$$ $ $$$$ $$$$$$$$$$ $$$$$$$$ $$ $$$$$$ $$$$$$$$$ $ $$$$$$$$$ $$$$$$$, $$ $ $$$$$$$$$ $$$$$$$$$ $$$$$$$$$. $$$ $$$$$$$$$$$$$$$$$$ $$$$$$$$$$ $$$$$$$$ $$$$$$$ $$$$$-$$$$$$$, $ $$$ $$$$$$$$$$$$ $$$$$$$$$$ — $$$$$$$$$$$$$ $$$$$$$. $ $$$$$$ $$$$$$ $$$$$ $$$$$$$$$$$ $$$ $$$$$$$$$ $$$$$$$: $$$ $$$ ($$$$$$ $$$$$$$, $$$ $$$$$$$$$$ $$$$$$$$$$), $$$ $$$ ($$$$$$$ $$$$$-$$$$$$$$$ $$$$$-$$$, $$$ $$$$$$$$$$$$ $$$$ $$$$$$ $$$$$$$$$$$$$ $$$$$$$$$) $ $.$$ $$$ ($$$$$$$ $$$$$-$$$$$$$$$ $$$$$$$$-$$, $$$ $$$$$$$$$$ $$$$$$$$ $$$$$$ $$$$$$$$$$$ $$$). $$$$$ $$$$ $$$$$$$ $$$$$$$$$ $$$$$$$$ $$$ $$$ $$$$$$$$ $$$$$$$$$ $$$$$$$$$$$$$$ $ $$$$$$$$$$$$$$$$ $$ $$$$$$$$$$$$$ $$$$$ $ $$$$$$$$$$$$ $$$$$$$$ $$$$$$$$$$$$$.

$$$ $$$$$$$$$$$$$ $$$$ $$$$$$$$$$ $ $$$$$$$$$$ $$$$$$$$$ $$$$$$$$$$ $$$$$$$$$ $$$$$$$$$$$$ $$$$$$$ $$$$$$$$$$$$$. $$$$$$$$$ $$$$$$$$$ $$$$$$ $$$$$$$$$$$$$$$ $ $$$$$$$$$ $$$$$$$$$$, $$$$$$ $$$$$$$ $$$$$$$$$$$$$ $$ $$$$$$$$$$$$$$$$ $$$$$$$$$ $$$$$ $$$$ $$$$$$$$ $$$ $$$$$ $$$$$$$$$$$$$ $$ $$$$$$$ $$$$$$$$$$$$$$$$$$$ $$$$$$$$$, $$$$$$ $ $$$$$$$$$$$$$$$$ $$$$$$$$ $$$$$$$$$$$$$$. $$$$$$ $$$ $$$$ $$$$$$$$$$ $$$$$$$$$ $$$$$$$$$ $$$$$$$$$$ $$$$$$$, $$$$$$$$$ $ $$$$$$$$$$ $$$$$$$$$$ $$$ $$$$$$$$$$ $$$$$$$ $$$$$$$$$ $$$ $$$$$$$$$ $$$$$$$$$. $ $$$$$$ $$$$$$ $$$$$$ $$$$$ $$$$$$$ $$ $$$$$$$$$$$$$$$, $$$$$$$$$ $$$$$$$$$$$$ $$$$ $$$$$$$$$ $ $$$$$$$$ $$$$$$$$$.

$$$$$ $$$$$$$, $$$$$ $$$$$-$$$$$, $$$$$$$$$$$$$ $ $$$$$$$$$ $$$$, $$$$$$$$$$$$$ $$$$$$ $ $$$$$$ $$$$$$$$$$ $$$ $$$$$$$$$$$$$ $$$$$$$ $$$$$$$$$$$$$ $$$$$-$$$$$$$$$ $ $$$$$$$$$ $$$$$$$$$$. $$$$$$$$$$ $$$$$ $$$$$$$$$$ $$$$$$$, $$$$$$$$$$ $$$$$$$$$ $ $$$$$$$ $$$$$$$$$$$ $$$$$$$$$$$ $$$$$$$$$ $$$$$$$$ $$$$$$, $$$$$$$$$ $$$ $$$$$$$$$$$$$$$ $$$$$$$ $ $$$$$$$$$. $$$$$$$$$$$$$$ $$$$$$$ $$$$$$ $$$$$$$ $$$$$$$$$$$ $$$$$$$$ $$$$$$$$$$$$ $ $$$$$$$$$$$$$ $$$$$$$$$$$$ $$$$$ $$$$$$$$$$$$ $$$$$$. $$$$$$$$$$$$ $$$$$$$$$$ $$$$$$$$$ $$$$$$$$$$$ $$$$$$$$ $$$$$$$ $$$$$$$$$$$$$, $$$ $$$$$$$$ $$$$$ $$$ $$$$$$$ $$$$$$ $$$$$$$$. $$$$$$$$$ $$$$ $$$$ $$$$$$$$ $$$$$$$$$$ $$$ $$$$$$$$$ $$$$$$$$$$ $$$$$$$$$$$$ $$$$$ $$$$$$ $ $$$$$$$$$ $$$$$$$$$$$ $$$$$$$$$$$.

$$$$$$$$$$$$$ $ $$$$$$ $$$$$$$ $$$$$$$$ $$$$$$ $$$$$-$$$$$ $ $$$$$$$$$$$ $$$ $$$$$$$$$$ $ $$$$$$$$$ $$$$ $$$$$$$$$$$ $$$$$$$$$$$$$ $$$$$$ $$$ $$$$$$$$$$ $$$$$$$$$$$$$$ $$$$$$$$$$$$$. $$$$$$ $$$$$$$$$$ $$$$$$$$$$$$$, $$$$$$$ $$$$$$$$$$ $$$$$$$$$ $ $$$$$$$$ $$$$$$$$$$$ $$$$$$$$$$$ $$$$$$$$$ $$$$$$$$$ $$$$$$$$$$$$ $$$$$ $$$$$$$$ $ $$$$$$$$$$ $$$$$$$$$ $$$$$$$$$$$$$ $$$$$$$$ $$$$$$. $$$$$$ $$$$$$$$ $$$$$$$ $$$$$$$ $$$$$$ $$$$$$$$$$ $$$$$$$, $$$$$$$$$ $$ $$$$$ $$$$$$$ $$$$$$$$ $$$$$$$$$$$$$$$ $$$$$$$$ $$$$$$$$$ $$$$$$$$$$$$$$ $$$$$-$$$$$$$$$ $ $$$$$$$$$ [$$]. $$$ $$$$$$$$$$$$$ $$$$$$$ $$$$$$$$$$$$$ $ $$$$ $ $$$$$$ $$$$$$$$$$ $$$$$$$$$ $$$$$$$$ $ $$$$$$$$$$ $$$$$$$$$$$$$$, $$$$$$$$$$$$$ $$$$$$$ $$$$$$$$ $$$$$$$$ $$$$ $$$$$$$$$$ [$$]. $$$$$$$$$ $$$$$$$$$$$ $$$ $$$$ $$$$$$$$$ $$$$$$$ $$$$$-$$$$$$$$$ $$$$$$$$ $$$$$$$ $$$$$$$$$$$$$$ $$$$$$$$$$$$ $$$$$$$ $$$$$ $ $$$$$$$ $$$$$ $$$$$$$ $$ $$$$ $$$$$$$$ $$$$$$$$$ $$$$$$$$$$$$$$ [$].

Алгоритм работы в программе GATE и анализ полученных результатов

Программа GATE (Geant4 Application for Tomographic Emission) представляет собой специализированное программное обеспечение с открытым исходным кодом, предназначенное для моделирования процессов переноса ионизирующего излучения в различных средах. GATE основана на библиотеке Geant4, разработанной в CERN, и предоставляет удобный интерфейс для создания геометрических моделей, задания источников излучения, настройки физических процессов и регистрации результатов. Основным языком описания сценариев моделирования в GATE является макроязык, основанный на командах, что позволяет гибко управлять всеми параметрами вычислительного эксперимента без необходимости программирования на C++.

Алгоритм работы в программе GATE включает несколько последовательных этапов. Первым этапом является создание геометрии моделируемой системы. Для этого используются команды, определяющие размеры, форму и положение объектов. В GATE поддерживаются различные геометрические примитивы: параллелепипеды, сферы, цилиндры, конусы и более сложные формы, создаваемые с помощью булевых операций. Каждому объекту присваивается определенный материал, который может быть выбран из встроенной базы данных или создан пользователем путем задания элементного состава и плотности. Для целей данной работы создаются два типа фантомов: водный фантом размерами 30×30×30 см и свинцовый фантом размерами 10×10×10 см. Выбор размеров обусловлен необходимостью обеспечения условий полного ослабления излучения в пределах фантома для всех рассматриваемых энергий.

Вторым этапом является задание источника излучения. В GATE источник описывается набором параметров: тип излучения (гамма-кванты), энергия (моноэнергетическая или спектральная), форма и размеры источника, угловое распределение испускаемых частиц, активность или количество моделируемых частиц. Для данной работы используются точечные изотропные источники, расположенные на расстоянии 10 см от поверхности фантома. Рассматриваются три значения энергии гамма-квантов: 100 кэВ, 662 кэВ и 1.25 МэВ. Для каждого значения энергии проводится отдельное моделирование, а также моделирование комбинированного источника, включающего все три энергии с равными весовыми коэффициентами.

Третьим этапом является настройка физических процессов. В GATE выбор физической модели осуществляется с помощью команд, определяющих набор процессов для каждого типа частиц. Для гамма-излучения в данной работе используется набор "Livermore", который обеспечивает высокую точность описания фотоэффекта, комптоновского рассеяния, эффекта образования пар и релеевского рассеяния в диапазоне энергий от 250 эВ до 100 ГэВ. Для электронов и позитронов также используются модели "Livermore", включающие процессы ионизации, тормозного излучения и многократного рассеяния. Пороговые энергии для генерации вторичных частиц (production cuts) устанавливаются равными 1 мм для всех материалов, что обеспечивает разумный баланс между точностью и вычислительными затратами.

Четвертым этапом является настройка регистрации результатов. Для расчета дозовых распределений используется актер "DoseActor", который создает трехмерную сетку (воксельную решетку) внутри фантома. Размер вокселей выбирается равным 1×1×1 мм для водного фантома и 0.5×0.5×0.5 мм для свинцового фантома, что обеспечивает достаточное пространственное разрешение для анализа дозовых профилей. Актер регистрирует поглощенную дозу в каждом вокселе, а также статистическую погрешность. Дополнительно могут быть зарегистрированы энергетические спектры частиц, пересекающих заданные поверхности, и распределения флюенса.

Пятым этапом является запуск моделирования. Количество моделируемых первичных частиц выбирается исходя из требуемой точности. Для данной работы используется 10⁷ частиц для каждого варианта моделирования, что обеспечивает статистическую погрешность менее 1% в центральной области фантома. Моделирование выполняется с использованием параллельных вычислений на 8 ядрах процессора, что позволяет сократить время расчета до нескольких часов. После завершения моделирования $$$$$$$$$$ $$$$$$$$$$$ в $$$$$ $$$$$$$ $$$ ($$$$$$$$$) $$$ $$$$ для $$$$$$$$$$$$ $$$$$$$.

$$$$$$ $$$$$$ $$$$$$$$ $$$$$$ $$$$$$$$$$ $$$$$$$$$$$. $$$$$$$$$ $$$$$$$$$$$$$ $$$$$$$ $$$$$$$$ $$$$$$$$$$ $$$$$$$$ $$$$ $$$$$ $$$$$$$$$$$ $$$ $$$$$$$ ($$$$$$$$$ $$$$$$$$$$$$$) $ $$$$$$$$$$ $$$$$$$ $$$$$$$ $$$$$$$$$$$$$ $$ $$$$$$$$$ $$$$$$$$. $$$ $$$$$$$$$$$$ $$$$$$$$$$$$ $$$$$$$$$ $$$$, $$$$$$$$$$$$$$$ $$$$$$$ $$$$$$$$$$$ $$$ $$$$$$$$$$ $$$$$$$$, $$$$$$$$$$ $ $$$$$$$$$$ $$$$$$$$$$$. $$$$$$$$$ $$$$$$$$$$$ $$$ $$$$$$$$$ $$$$$$$$$$ $ $$$$$$$ $$$$$$$$$$ $$$$$ $$$$$$$$$$ $$$$$$$ $$$$$$$$$$$$$ $$ $$$$$$$$$$$$ $$$$$$$$ $$$ $$ $$$$ $ $$$$$$$$ $$$$$.

$$$ $$$$$$$ $$$$$$$ $$$$$$$$$$$$$ $ $$$$$$ $$$$$$$ $$$$$$$$$, $$$ $$$ $$$$$$$ $$$ $$$, $$$ $$$$$$$$$$ $$$$$$$$$$, $$$$$$$$$$ $$$$$ $$$$$$$$$$$ $$$$$$$$$$$$$$$ $$$$$$ $$$$$$$$$$$ $$$$$$$, $ $$$$$$$ $$$$$$$$$$$$$ $$$$$ $$$$$ $$$$$$ $$$$ $ $$$$$$$$. $$$ $$$$$$$ $$$ $$$, $$$ $$$$$$$$ $$$$$$$$$$ $$$$$$$$ $$$$$$$$$$$$$ $$$$$$$$$, $$$$$$$ $$$$$$$$$$$$$ $$$$$ $$$$$ $$$$$$$, $ $$$$$$$$$$ $$ $$$$$$$$$ $$$$$$$ ($$$$$$ $$$$$$$$$$ $$$$). $$$ $$$$$$$ $.$$ $$$, $$$ $$$$$$$$$$ $$$$$$$$$$$$ $$$$$$ $$$$$$$$$$$ $$$, $$$$$$$ $$$$$$$$$ $$$$ $$$$$ $$$ $$$$$$, $ $$$$ $$ $$$$$$$$$$ $$$$$ $$$$$ $$$$$$$.

$ $$$$$$$$$ $$$$$$$ $$$$$$$$$ $$$$$ $$$$$$$ $$$$$$$$$$ $$$$$$$$$ $$ $$$$$$$$$ $ $$$$$, $$$ $$$$$$$$$$$ $$$$$$$ $$$$$$$ $$$$$$$ $$$$$$. $$$ $$$$$$$ $$$ $$$ $$$$$$$$$$ $ $$$$$$ $$$$$ $$$$$$$$$ $$$$$$$, $$$ $$$$$$$$$$$ $$$ $$$$$-$$$$$$ $$$$$ $$$$$$$$$ $ $$$$$$ $$$$$$$$$$$$$$$$ $$$$ $$$$$$$$ $$$$$ $ $$. $$$ $$$$$$$ $$$ $$$ $$$$$$$ $$$$$$$$$$$$$ $$$$$ $$$$$$, $$ $$$ $$$$$ $$$$$$$$$$$ $$$$$$, $$$ $ $$$$. $$$ $$$$$$$ $.$$ $$$ $$$$$$ $$$$$$$$$$$ $$$ $ $$$$$$ $$$$$ $$$$$$ $$$$$$$$$$$$$$ $$$$$ $ $$$$$$$$$$, $$$ $$$$$ $$$$$$$$ $ $$$$$$$$$$ $$$$$$$ $$$$$$$$$$$$$ $$ $$$$$$$$$ $ $$$$$.

$$$$$$$$$ $$$$$$$$$$$ $$$ $$$$ $$$$$$$$$$ $ $$$$$$$$$$ $$$$$$$$$ $$$$$$$$ $$$$$$$ $$$$$$$$$$$$$$ $$$$$$$$$$$$ $$$$$$$ $$$$$ $ $$$$$$$$$$$ $$ $$$$$$$ $$$$$$$$$ $ $$$$$$$$$ $$$$$$$. $$$$$$$$$, $$$ $$$ $$$$$$$$$$$$$$$$ $$$$$$$$$ $$$$$$$ $$$$$$$$$$$$$ $$$$$ $$$$$$$$$$$$ $$$$$ $$$$$$$$$$$$ $$$$$$$$$$$$$ $$$ $$$$$$$$$ $$$$$$$, $$$$$$ $ $$$$$$$$$$$$$$$$ $$$$$$$$ $$$$$$$$$$$$$$. $$$$$$ $$$ $$$$ $$$$$$$$$$ $$$$$$$$$, $$$ $$$$$$$$$ $$$$$$$, $$$$$$$$$$$$ $$$$$$$$ $$$$$ $$$$$$$, $$$$$ $$$$$$$$$$$$$$$$$ $ $$$$$$$$ $$$$$$ $$$$$$$, $$$ $$$$$ $$$$$$$$$ $ $$$$$$$$$ $$$$$$$$$$$ $$ $$$$$$$$ $$$$$$$$$$$$.

$$$$$$ $$$$$$$$ $$$$$$$ $$$$$$$$ $$$$$$ $$$$$$$$$$$$$$ $$$$$$$$$$ $$$$$$$$$$ $$$$$$$$$$$. $$$ $$$$$ $$$$$$$$$$$$$$ $$$$$$$$$$$ $$$$$$$$$$ $$$ $$$$$$ $$$$$ $$$$$$$$ $$$$$$$$$$$$$ $ $$$$$$$$ $$$$$$$$$$$$$ $$$$$$$$$. $$$$$$$$$ $$$$$$$ $$$$$$$$ $$$ $$$$$$$$$ $$$$$$$ $$$$$$$$$$ $ $$$$$$$$$$$$$$ $$$$$$$$$ $$$$$$$$, $$$$$ $$$ $$$$$$$$ $$-$$$$$$$ $$$ $$$$$-$$$$$$, $$$$$$ $$$$$$$$$$$$ $ $$$$$$$$$$$ $$$$$$ $$$ $$$$$$$$$$$ $$$$$$ $$$$$$$ $$$$$$$. $$$$$-$$$$$$ $$$$$$$$$ $$$$$$$ $$$$$$$ $$$$$$$$$$ $$$$ $$$$$$$ $$$$$$$$$$$$$ $ $$$$$$ $$$ $$$$$$$, $$$ $ $$$$$$$$$$$$$$$$ $$$$$$$$.

$$$$$ $$$$$$$, $$$$$$$$ $$$$$$ $ $$$$$$$$$ $$$$ $$$$$$$$ $$$$$$$$$$$$$$$$$$ $$$$$$ $$ $$$$$$$$ $$$$$$$$$ $$ $$$$$$$ $$$$$$$$$$$, $$$$$$ $$ $$$$$$$ $$$$$$$ $$$$$$$$$$ $$$$$$$$$ $ $$$$$$$$. $$$$$$$$$$ $$$$$ $$$$$$$$$$ $$$$$$$$$$$$$, $$$$$$$ $$$$$$$$$$ $$$$$$, $$$$$$$ $$$$$ $ $$$$$$$$$$ $$$$$$, $$$$$$$$ $$$$$$$$ $$$$$$$$ $$$ $$$$$$$$$ $$$$$$$$$$$ $$$$$$$$$$$. $$$$$$ $$$$$$$ $$$$$$$$$$$$$ $$$$$$$$$ $$$$$$$ $$$$$$$$$$$$$$, $$$$$$$$$ $ $$$$$$$$ $$$$$$$$$ $ $$$$$$$$$$ $$$$$$$, $ $$$$$$$ $$$$$$$$$$$$ $$$$$$ $ $$$$$$$$$$ $$$$$$$$$$$$$$ $$$$$-$$$$$$$$$ $ $$$$$$$$$ [$]. $$$$$$$$$$ $ $$$$$$$$$$ $$$$$$$$$$$$$ $$$$$$ $$$$$ $$$$ $$$$$$$$$$$$ $$$ $$$$$$$$$$$ $$$$$$$$$$$$$ $$$$$$$$ $ $$$$$$$$$$$ $$$$$$$$$$ $$$$$$$$$$$$ $$$$$$ [$$]. $$$$$$$$$ $$$$$$$$$$$ $$$ $$$$$$$$$ $$$$$$$$$$ $ $$$$$$$ $$$$$$$$$ $$$$$$$ $$$$$$$$$$$$$ $$$$$$ $$$ $$$$$$$$$ $$$$$$$$$ $ $$$$$$$ $$$$$$$$$$$ $$$$$$$$$$$$ $$$$$$$ $$$$$ $ $$$$$$$$$$$$$ $$$$$$ [$$].

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

Для водного фантома глубина максимума дозы для энергии 100 кэВ составляет примерно 0.5-1 мм, что обусловлено малым пробегом фотоэлектронов. Для энергии 662 кэВ глубина максимума увеличивается до 2-3 см, а для энергии 1.25 МэВ — до 4-5 см. За максимумом дозовое распределение спадает экспоненциально с показателем, определяемым линейным коэффициентом ослабления материала. Для свинцового фантома глубина максимума дозы значительно меньше из-за большего коэффициента ослабления. Для энергии 100 кэВ максимум дозы находится практически на поверхности, для энергии 662 кэВ — на глубине около 1-2 мм, для энергии 1.25 МэВ — на глубине около 3-5 мм.

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

При моделировании комбинированного источника, состоящего из трех моноэнергетических компонент с равными интенсивностями, дозовое распределение представляет собой сумму распределений для каждой энергии. Однако при анализе такого распределения необходимо учитывать, что вклад каждой компоненты в общую дозу зависит от глубины. На малых глубинах преобладает вклад низкоэнергетической компоненты (100 кэВ), которая эффективно поглощается вблизи поверхности. На средних глубинах (2-5 см) основную роль играет компонента 662 кэВ, а на больших глубинах — компонента 1.25 МэВ. Таким образом, результирующее дозовое распределение имеет более сложную форму, чем для моноэнергетического источника.

Важным инструментом анализа является расчет так называемого "качества пучка" (beam quality), которое характеризует степень проникновения излучения в среду. Для моноэнергетических источников качество пучка однозначно определяется энергией фотонов. Для комбинированных источников качество пучка зависит от спектрального состава и может быть охарактеризовано эффективной энергией или средним линейным коэффициентом ослабления. В данной работе для оценки качества пучка используется глубина, на которой доза уменьшается вдвое по сравнению с максимальным значением (половинная глубина). Для водного фантома половинная глубина для энергии 100 кэВ составляет около 3-4 см, для энергии 662 кэВ — около 10-12 см, для энергии 1.25 МэВ — около 15-18 см. Для комбинированного источника половинная глубина занимает промежуточное значение.

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

Анализ энергетических спектров гамма-излучения внутри фантома позволяет получить дополнительную информацию $ $$$$$$$$$ $$$$$$$$$$$$$$. $$$ $$$$$$$$$$$$$$$$$$$ $$$$$$$$$ $$ $$$$$ $$$$$$$$ $$$$$$ $$$$$$$$ $ $$$$$$$$ $$$ $$$$$$$ $$$$$$$$$$ ($$$$$$$) $ $$$$$$$$$ $$$$$$$$$$ $$$$$$$$$$ $$$$$$$. $ $$$$$$$$$$$ $$$$$$$ $$$$ $$$$$$$$$$ $$$$$$$ $$$$$$$$$$, $ $$$$$$$ $$$$$$$$$$. $$$ $$$$$$$ $$$ $$$ $ $$$$$$$ $$$$$ $$$$$$$$$$$$ $$$ $$$$$$$$$ $$$$$$$$$ $$$ $$$$$$$ $$$$$ $$$ $$$, $$$$$$$$$$$$$ $$$$$$$$$$$$$ $$$$$$$$$$ $$ $$$$$$$ $$$$. $$$ $$$$$$$ $.$$ $$$ $ $$$$$$$ $$$$$$$$$$ $$$$ $$$$$$$$$$$$$$$$ излучения $ $$$$$$$$ $$$ $$$, $$$$$$$$$$$$$ $$$$$$$$ $$$$$$$$$$$ $$$.

$$$ $$$$$$$$$$ $$$$$$$ $$$$$$$$ $$$$$$$$ $$$$$$$$$$$ $$$$$$$$$$. $$$$$$$$$$ $$$$$$$$ $$$$$$$$$$$ $$$$$$$ $$$$$$ $$$$$$$$ $ $$$$$$$$, $ $ $$$$$$$ $$$$$$$$$$ $$$$$$$$$$$$$$$$$$ $$$$$ $$$$$$ $ $$$$$$$$$ $$-$$ $$$ ($-$$$$$). $$$ $$$$$ $$$$$$$$$$$ $$$$$$$$$$$$$$ $$$$$$$$$$, $$$$$$$$$$$ $$$$$ $$$$$$$$$$$$$$$$$$ $$$$$$$$$$. $$$$$$$ $$$$$$$$$$$$$$ $$$$$ $$$$$$$$ $$$$$$ $$$$$$$$$$$$$$$ $$$$$$$$$, $$$$$$$$$$$ $$$$$$$$$$$$$$$$ $$$$$$$$ $$ $$$$$$$ $$$$$$$$$$$ $$$$$$$$$.

$$$ $$$$$$$ $$$$$$$$$$$$$$ $$$$$$$$$$$ $$$$$$$$$$$ $$$$$$$$$$ $$$$$$$$$, $$$ $$$$$$$$$$$ $$$$$$$$$$$ $ $$$$$$ $$$$$$ $$$$$$$. $ $$$$$$$ $$$$$$$$$ $$$$, $$$ $$$$$$$$$$$$$$ $$$$$$$ $$$$$$$$$$ $$$$$$$, $$$$$$$$$$$ $$$$$$$$$$ ($$$$$ $%). $$ $$$$$$$ $$$$$$$$, $$$ $$$$ $$$$, $$$$$$$$$$$ $$$$$$$$$$ $$ $-$$% $ $$$$$. $$$ $$$$$$$$$ $$$$$$$$ $$ $$$$$$$ $$$$$$$$ $$$$$ $$$$$$$$$$$$ $$$$$$ $$$$$$$$$$ $$$$$$$$$$, $$$$$ $$$ $$$$$$$$$$ $$$$$ $$$$$$$$$$$$ $$$$$$ $$$ $$$$$$$$$$$$$ $$$$$$$ $$$$$$$$$$ $$$$$$$$$.

$$$$$$ $$$$$$ $$$$$$$ $$$$$$$$ $$$$$$$$ $$$$$$$$$$ $$$$$$$$$$$. $$$ $$$$$ $$$$$$$$$$ $$$$$ $$$$$$$$$$$$$ $ $$$$$$$$$ $$$$$$ $$$$$$ ($$$$$$$$, $$$, $×$$$, $$$) $ $$$$$$$$$$$$ $$$$$$$$$$ $$$$$$$ $$$$$$$$$$$$$. $$$$ $$$$$$$$$$ $$$ $$$$$$ $$$$$ $$$$$$ $$$$$$$$$ $ $$$$$$$$ $$$$$$$$$$$$$$ $$$$$$$$$$$, $$ $$$$$ $$$$$$$, $$$ $$$$$$$$$$ $$$$$$$$$$. $$$$$$$$$$$$$ $$$$$ $$$$$$$$ $$$$$$$$$$$$$ $ $$$$$$$$$$ $$$$$$$$$$$$ $$$$$$$$$ $$$$$ $ $$$$$$$$$, $$$ $$$$$$$$$$ $$ $$$$$$$ $$ $$$$$$ $$$$$$$$$$.

$$$ $$$$$$$$$$$$$ $$$$$$$$$$$ $$$$$$$$$$ $$$$$ $$$$$$$$$ $$$$$$$$$ $$$$$$$$$$$$$$$ $$$$$$$$$$$, $$$$$$$$$ $ $$$$$$$$$$$$$, $$$$$$$$$$$$$ $ $$$$$$$$$$ $$$$$$$. $$$$$$$$, $$$$$$ "$$$$$$$$$" $$ $$$$$$$$$ $$$$$$$, $$$$$$$$$ $ $$$$$$$$$$$$ $$$$$$$$$, $$$ $$$$$ $$$$$$$$$ $ $$$$$$$$$ $$$$$$$$$$$$ $$$ $$$$$$$ $$$$$$$ $$$$$$$$$$$$$ $$$$$$$$$$ $$$$$$$. $$$$$$ "$$$$$$$$" $$$$$$$$ $$$$$ $$$$$$, $$ $$$$$$$ $$$$$$ $$$$$$$$$$$$$$ $$$$$$$$. $ $$$$$$ $$$$$$ $$$$$$$$$$$$ $$$$$$ "$$$$$$$$$", $$$$$$$ $$$$$$$$$$$$ $$$$$$$$$$$ $$$$$$$$ $$$ $$$$$$$$$$$$ $$$$$ [$].

$$$$$$$$$ $$$$$$$$$$ $$$$$$$$$$$ $ $$$$$$$$$$$$$ $$$$$$$ $$$$$$$$ $$$$$$ $$$$$$ $$$$$$$$$$$ $$$$$$$$$$$$$. $$$ $$$$$$$ $$$$$$$ $$$$$$$$$$ $$$$$$$$$ $$$$$$$$$$$$$$$$$ $ $$$$$$$$$ $$$$$$ $$ $$$$$$$ $$$$$$$$$$$$$$ $$$ $$$$$$$$$ $$$$$$$ $$$$$-$$$$$$$$$. $$$$$$$$$ $ $$$$$ $$$$$$$ $$$$$$$$$ $$$$$$$ $$$$$$$$$$$$ $$$$$$$$$$$$ $$$$$$ $ $$$$$$$ $$$$$$$$$ $$$$$$ $ $$$$$$$$$$. $$$ $$$$$$$$$$ $$$$$$$ $$$$$$$$$$$$ $$$$$$ $$$$$ $$$$$$$$$$$$$$, $$ $$$$$ $$$$$$$$ $$$ $$$$ $$$$$$$.

$$$$$ $$$$$$$, $$$$$$ $$$$$$$$$$ $$$$$$$$$$$ $$$$$$$$ $$$$$$$$$$ $ $$$$$$$$$$$$$ $$$$$$$$$ $ $$$$$$$$$$ $$$$$$$ $$$$$$$$, $$$$$$ $$$$$$$$$$$$$$$$$$ $$$$$$$$$$ ($$$$$$$ $$$$$$$$$, $$$$$$$$$$ $$$$$$$), $$$$$$ $$$$$$$$$$$$$$ $$$$$$$$, $$$$$$ $$$$$$$$$$$$$$ $ $$$$$$$$$$$$$$$ $$$$$$$$$$$$, $ $$$$$ $$$$$$$$$ $ $$$$$$$$$$$$$ $$$$$$$. $$$$$$$$$$$ $$$$$$ $ $$$$$$$ $$$$$$$$$ $$$$$$$$ $$$$$$$$$$$$ $$$$$$$$$$$$$ $ $$$$$$$$$ $$$$$$$$$$$$ $$$$$$$ $$$$$ $ $$$$$$$$$ $$$$$$$$$$ $ $$$ $$$$$$$$$ $$$$$$$$ $$$$$$$$$. $$$$$$$$$$ $ $$$$$$ $$$$$$ $$$$$$$$$$ $$$$$ $$$$ $$$$$$$$$$$$ $$$ $$$$$$$$$$$ $$$$$$$$$ $$$$$$$, $$$$$$$$$$$ $$$$$$$$$$ $$$$$$$$$$$$ $$$$$$ $ $$$$$$$$$$ $$$$$$$$$$$$ $$ $$$$$$$$$$ $$$$$$$$$ $$$$ $ $$$$$$$ $$$$$$$$$$ [$$].

Постановка задачи и создание геометрии фантома в программе GATE

Практическая часть данной курсовой работы посвящена непосредственному моделированию дозовых распределений гамма-излучения в программе GATE. На первом этапе необходимо четко сформулировать постановку задачи, определить геометрические параметры моделируемой системы, выбрать материалы фантомов и задать характеристики источников излучения. Корректная постановка задачи является основой для получения достоверных и воспроизводимых результатов, поскольку любые неточности в геометрии или физических параметрах могут привести к существенным ошибкам в расчетах.

В рамках данной работы рассматривается следующая постановка задачи. Имеется точечный изотропный источник гамма-излучения, расположенный на фиксированном расстоянии от поверхности фантома. Фантом представляет собой прямоугольный параллелепипед, заполненный однородным материалом. Рассматриваются два типа фантомов: водный фантом (H₂O, плотность 1.0 г/см³) и свинцовый фантом (Pb, плотность 11.34 г/см³). Выбор воды обусловлен ее широким использованием в качестве тканеэквивалентного материала в медицинской физике, а выбор свинца — его применением в радиационной защите. Размеры фантомов выбраны таким образом, чтобы обеспечить полное ослабление излучения в пределах фантома для всех рассматриваемых энергий. Для водного фантома размеры составляют 30×30×30 см, что достаточно для ослабления излучения с энергией до 1.25 МэВ. Для свинцового фантома размеры составляют 10×10×10 см, что также обеспечивает полное поглощение излучения благодаря высокому коэффициенту ослабления свинца.

Источник излучения располагается на расстоянии 10 см от поверхности фантома вдоль центральной оси. Такое расстояние выбрано для того, чтобы обеспечить формирование коллимированного пучка и избежать эффектов, связанных с геометрическим расхождением лучей. Источник является моноэнергетическим для каждого отдельного моделирования, с энергиями 100 кэВ, 662 кэВ и 1.25 МэВ. Дополнительно моделируется комбинированный источник, включающий все три энергии с равными интенсивностями. Для каждого варианта моделирования количество первичных частиц составляет 10⁷, что обеспечивает статистическую погрешность менее 1% в центральной области фантома.

Создание геометрии фантома в программе GATE осуществляется с помощью макрокоманд, которые определяют размеры, положение и материал каждого объекта. Для водного фантома используется команда, создающая параллелепипед с заданными размерами и присваивающая ему материал "Water" из встроенной базы данных GATE. Для свинцового фантома материал "Lead" также доступен в базе данных. При необходимости можно создать пользовательский материал, задав его элементный состав и плотность. В данной работе используются стандартные материалы, что упрощает настройку и обеспечивает воспроизводимость результатов.

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

$$$$$$$$$$$ $$$$$$$$$$$ $$$$$$$$$$$$$$ $ $$$$$$$ $$$$$$ "$$$$$$$$$", $$$$$$$ $$$$$$$ $$$$$$$$$$ $$$$$$$$$$ $$$$$ $$$$$$ $$$$$$$. $$$ $$$$$$$ $$$$$$$ $$$$$$ $$$$$$$$ $$$$$$$$$$ $×$×$ $$, $$$ $$$$$$$$$$$$ $$$$$$$$$$$ $$$$$$$$$$$$$$$$ $$$$$$$$$$ $$$ $$$$$$$ $$$$$$$ $$$$$$$$. $$$ $$$$$$$$$$ $$$$$$$, $$$ $$$$$$$ $$$$$$$$$$$$$ $$$$$ $$$$$ $$$$$$ $$$$$$$$$, $$$$$$ $$$$$$$$ $$$$$$$$$$$ $$ $×$×$ $$. $$$$$ $$$$$$$$$$$$ $$$$$$$$$$$ $$$$ $ $$$$$$ $$$$$$$, $ $$$$$ $$$$$$$$$$$$$$ $$$$$$$$$$$. $$$$$$$$$$ $$$$$$$$$$$ $ $$$$$ $$$$$$$ $$$, $$$$$$$ $$$$$ $$$$ $$$$$$$$$$$$$ $ $$$$$$$$$ $$$ $$$$$$$$$$$$ $ $$$$$$$, $$$$$ $$$ $$$$ $$$ $$$$$$$$.

$$$$$ $$$$$$$$$$ $$$$$$$$$$$$$ $$$$$$$$$$ $$$$$$$$$ $$$$$$ $$$$$$$$$$$. $$$ $$$$$$$ $$$$$$$$ $$$$$$$$$$$$$ $$$$$$$$ $$$$$$$$$ $$$$$$$ $$$$$$$$$$$$$ $$$$$ $$$$$$$$$$$ $$$ $$$$$$$. $$$ $$$$$ $$$$$$$$ $$$$ $ $$$$$$$$, $$$$$$$$$$$$$ $$ $$$$$$$$$$$ $$$, $$$$$$$$$$$ $$ $$$$$$$$$$$ $$$$$$$ $$$$$$$$ $ $$ $$$ $$$$$$$$$$ $$$$$$$$$$$$$$ $$$$$$$$$$. $$$$$$$$$$ $$$$$$ $$$$$$$$$$$ $$ $$$$$$$$$$$$ $$$$$$$$ $$$$, $$$ $$$$$$$$$ $$$$$$$$$$ $$$$$$$$$$$$$ $$$ $$$$$$$$$ $$$$$$$ $ $$$$$$$$$$. $$$$$$$$$$$$$ $$$$$$$$ $$$$$$$$$$ $$$$$$$ $$$$ $$ $$$$$$$$$$ $$$$$$$$ $$$ $$$$$$ $$$$$$$ $$$$$$$$$ $$ $$$$$ $$$$$$$$ $$$$.

$$$$$$ $$$$$$ $$$$$$$$ $$$$$$$$ $$$$$$$$$$$$ $$$$$$$$$$$$$. $$$ $$$$$ $$$$$$$$$$ $$$$$$$$$ $$$$$$$$$$ $$$$$$$$$$$ $ $$$$$$$$$$ $$$$$$$$$$$$$ $$$$$$$. $$$ $$$$$$$ $$$$$$$ $$$$$$$$$$ $$$$$$$$$ $$$$$$$$$$$$$$$$$ $ $$$$$$$$$ $$$$$$ $$ $$$$$$$ $$$$$$$$$$$$$$ $$$ $$$$$$$$$ $$$$$$$ $$$$$-$$$$$$$$$. $$$$$$$$$ $ $$$$$ $$$$$$$ $$$$$$$$$ $$$$$$$$$ $ $$$$$$$$$$$$ $$$$$$$$ $$$$$$ $ $$$$$$$ $$$$$$$$$ $$$$$$. $$$ $$$$$$$$$$ $$$$$$$ $$$$$$$$$$$$ $$$$$$ $$$$$ $$$$$$$$$$$$$$, $$ $$$$$ $$$$$$$$ $$$ $$$$ $$$$$$$. $ $$$$$$ $$$$$$$$$$$ $$$$$$$$$$$$ $$$$$$$$$$$ $$$$$$$$$$ $$$$$$ $$$$$$$$$ $$$$$$, $$$$$ $$$ $$$$$$$$$$$$ $$$$$ $$$$$$$$$$ $$$$$$$$$, $$$$$$ $ $$$$$$$$$ $$$ $$$$$$$$$$$$$ $$$$$$$$$$.

$$$$$ $$$$$$$, $$$$$$$$$$ $$$$$$ $ $$$$$$$$ $$$$$$$$$ $$$$$$$ $$$$$$$$ $$$$$$$$$ $$$$$$$ $$$$$$$$$$$$ $$$$$ $$$$$$. $$$$$$$$$$ $$$$$ $$$$$$$$ $$$$$$$$, $$$$$$$$$$ $$$$$$$$$ $ $$$$$$$$$$ $$$$$$$ $$$$$$$$$$$$ $$$$$$$$$ $$$$$$$$$$$ $$$$$$$$$$$, $$$$$$$$$ $$$ $$$$$$$$$$$$$$$ $$$$$$$. $$$$$$$$$$$$$ $$$$$$$$$$$ $$$$$$$$$$ $$ $$$$ $$$$$$ $$$$ $$$$$$$$ $$$$$$$$$ $ $$$$$$$$ $$$$$$$$$$$$$$$$$ $$$$$$$$$$$$$. $$$$$$$$$$$ $$$$$$$$$$$ $ $$$$$$$ $$$$$$ "$$$$$$$$$" $ $$$$$$$$$$$ $$$$$$$$$$$$$$$$ $$$$$$$$$$$ $$$$$$$$$ $$$$$$$$ $$$$$$$ $$$$$$$ $$$$$$$$$$$$$ $ $$$$$$$ $$$$$$$$$$$$$$, $$$$$$$$$ $ $$$$$$$$ $$$$$$$$$ $ $$$$$$$$$$ $$$$$$$ [$$]. $$$$$$$$ $$$$$$$$$$$$ $$$$$$$$$$$$$ $$$$$ $$$$$$$$$ $ $$$$$$$$$$$$$ $$$$$$$ $$$$$$$$ $$$$$$$$$$$$ $$$$$$, $$$$$$$$$$$$$$ $$$$$$$$$$ $$$$$$$$$$ $$$$$$$ [$]. $ $$$$$$$$$ $$$$$$$ $$$$$ $$$$$$$$ $$$$$$$$$$$ $$$$$$$$$$ $$$$$$$$$$$$$ $$$ $$$$$$$ $$ $$$$ $$$$$$$$$$$$$$$$$$ $$$$$$$$$$ $ $$$$$$$$ $$ $$$$$$$$$$$$$ $$$$$$ [$$].

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

Для задания геометрии в GATE используется система команд, которая позволяет создавать объекты различных форм и размеров. В данной работе для создания фантомов используется команда, создающая прямоугольный параллелепипед (box). Размеры фантомов задаются в сантиметрах: для водного фантома — 30×30×30 см, для свинцового — 10×10×10 см. Положение фантома в пространстве определяется координатами его центра. Для удобства анализа фантом располагается таким образом, чтобы его передняя поверхность находилась в плоскости z = 0, а источник — на оси z на расстоянии 10 см от этой плоскости.

Материалы фантомов выбираются из встроенной базы данных GATE. Для водного фантома используется материал "Water", который имеет плотность 1.0 г/см³ и элементный состав, соответствующий химической формуле H₂O. Для свинцового фантома используется материал "Lead" с плотностью 11.34 г/см³. При необходимости можно создать пользовательский материал, задав его состав вручную. Однако в данной работе использование стандартных материалов является предпочтительным, поскольку они прошли верификацию в многочисленных исследованиях и обеспечивают корректное воспроизведение физических процессов.

Особое внимание уделяется настройке источника излучения. В GATE источник описывается с помощью команд, определяющих его тип, энергию, угловое распределение и активность. В данной работе используется точечный изотропный источник (type = "Point", distribution = "Isotropic"). Энергия источника задается как моноэнергетическая (type = "Monoenergetic") с одним из трех значений: 100 кэВ, 662 кэВ или 1.25 МэВ. Для комбинированного источника используется спектральное распределение, включающее все три энергии с равными весовыми коэффициентами. Количество моделируемых частиц для каждого варианта составляет 10⁷, что обеспечивает достаточную статистическую точность.

Для регистрации дозовых распределений используется актер "DoseActor", который создает трехмерную воксельную сетку внутри фантома. Размер вокселей выбирается исходя из требуемого пространственного разрешения. Для водного фантома, где дозовое распределение имеет относительно плавные градиенты, размер вокселей составляет 2×2×2 мм. Для свинцового фантома, где градиенты более резкие, размер вокселей уменьшается до 1×1×1 мм. Такое разрешение позволяет детально изучить форму дозовых профилей и выявить особенности, связанные с различными механизмами взаимодействия.

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

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

$$$ $$$$$$$ $$$$$$$ $$$$$$$$$, $$$ $$$$$$$ $$$$$$$$$ $$$$ $$$$$ $$$$$$$$$$$$$ $ $$$$$$ $$$$$$$ $$$$$$$$$. $$$ $$$$$$$ $$$ $$$ $$$$$$$$ $$$$$ $$$$$$$$$$ $$$$$$$$$$$ $$ $$$$$$$$$$$, $$$ $$$ $$$ — $$ $$$$$$$ $$$$$ $-$ $$, $$$ $.$$ $$$ — $$ $$$$$$$ $$$$$ $-$ $$. $$$$$$$$$$ $$$$$$$ $$$$$ $$$$$ $$$$$$$$$$$$$ $ $$$$$$$$, $$$ $$$$$$$$$$$$$$$ $ $$$$$$$ $$$$$$$$$$$ $$$$$$$$$$$ $$$$$$$$$$$$$$$$$$$$$ $$$$$$$$$. $$$ $$$$$$$$$$ $$$$$$$ $$$ $$$ $$$$$$$$$$$$$$ $$$$$ $$$$$$$$$$$ $$$$$$ $$-$$ $$$$$ $$$$$$$$ $$$$$$$$$$$$ $$$$$$$$$$.

$$$$$$$$$$ $$$$$$$ $$$$ $$$$$$$$ $$ $$$$$$$$$$ $$$$$$$$ $$$ $$$$$$ $$$$$$$ $$$$$$$$$ $$ $$$$$ $$$$$$$$ $$$$. $$$ $$$$$$$ $$$$$$$ $$ $$$$$ $$$$$$$$ $$$$$$$ $$$$$ $$$$$, $$$$$$$$$$$$$$$ $$$$$$$$$ $$$$$$$$$$ $$$$$. $ $$$$$$$$$$$ $$$$$$$ $$$$$$$ $$$$$ $$$$$$$$$$$ $$ $$$$ $$$$$$$$$$$ $$$$$$$$$, $$$$$$ $$$$$$$ $$$$$$$$$$ $$$$$ $$$$$$$$ $$ $$$$$$$. $$$ $$$$$$$ $$$ $$$ $$$$$$$$$ $$$$$$$$$$$$$, $ $$$$$$$ $$$$$$$$ $$$$$ $$ $$$$ $$$$$$$$. $$$ $$$$$$$ $$$ $$$ $ $.$$ $$$ $$$$$$$$$ $$$$$$$$$$$, $ $$$$$$$ $$$$$$$$$$$ $$$$$$$$$$$, $$$$$$$$ $$ $$$$$$$ $$$$$$$$. $$$ $$$$$$$$$$ $$$$$$$ $$$$$$$$$ $$$$$ $$$$$$$$ $$-$$ $$$$$$$$ $$$$$$$$$$, $ $$$$$$$ $$$$$$$$ $$$$$$$$$$$$ $$$$$.

$$$$$$ $$$$$$$$ $$$$$$$ $$$$$$$$ $$$$$$ $$$$$$ $$$$$$$$$ $$$$$$$$$$ $$$$$$$$$$$$$$ $ $$$$$$$$$$$$ $$$$$$$$ $$$$$$$$$$$$$. $$$ $$$$$ $$$$$ $$$$$$$$$$$ $$$$$$$$$$$$$$ $$$$$$$ $$$$$$, $$$$$$$$$$$$ $$$$$$$$ $$$$$$$$$$$ $$$$$$ $$$$$$$. $$$ $$$$$$$ $$$ $$$, $$$ $$$$$$$$$$ $$$$$$$$$$, $$$$$$ $$$$$ $$$$$$$$$ $ $$$$$$$$ $$$$$$$ $ $$$$$$$$$ $$$$$$$$$$ $$$$$$$$$$$$$$$ $$$$$$$$$. $$$ $$$$$$$ $$$ $$$, $$$ $$$$$$$$ $$$$ $$$$$$ $$$$$$$$$$$$$ $$$$$$$$$, $$$$$$ $$$$$ $$$$$$$$$ $$$$$$$ $$$$$$$$$ $$$$$$$$$$ $$$$$$$ $ $$$ $$$$$$$$$ $$$$$$$$$. $$$ $$$$$$$ $.$$ $$$, $$$ $$$$$$$$$$ $$$$$$$$$$$$ $$$$$$ $$$$$$$$$$$ $$$, $ $$$$$$$ $$$$$$$$ $$$$ $$$$$$$$$$$$$$$$ $$$$$$$$$ $ $$$$$$$$ $$$ $$$.

$$$$$ $$$$$$$, $$$$$$$$ $$$$$$$$$ $$$$$$$ $ $$$$$$$$$ $$$$$$$$$$ $$$$$$$$$$$$$ $$$$$$$$ $$$$$$$$$$ $$$$$$$, $$$$$$$$$$$$$ $$$$$$$$ $ $$$$$$$$$$$$$ $$$$$$$$$$ $$$$$$$$$$$. $$$$$$$$$$ $$$$$ $$$$$$$$ $$$$$$$$, $$$$$$$$$$, $$$$$$$$$$ $$$$$$$$$ $ $$$$$$$$$$ $$$$$$$ $$$$$$$$$$$$ $$$$$$$$$$ $$$$$$$$$$$$$$$ $$$$$$$$$ $$$$$$$$$$$$$$ $$$$$-$$$$$$$$$ $ $$$$$$$$$. $$$$$$$$$$$$$ $$$$$$$$$$ $$$$$ $ $$$$$$$$$$$ $$$$$$$$$$$$$$$$ $$$$$$$$$$$ $$$$$$$$$ $$$$$$$$ $$$$$$$ $$$$$ $$$$$$$ $$$$$$$$ $ $$$$$$$ $$$$$$$$$$$, $$$$$$$$$ $ $$$$$$$$$$ $$$$$$$$$$$ $$$$$$$$$$$$$$. $$$$$$ $$$$$$$$$ $ $$$$$$$$$$ $$$$$$$$, $ $$$$$ $$$$$$$$$$$$$$ $$$$$$$$ $$$$ $$$$$$$$$$$$ $$$$$$$$$$$$$ $ $$$$$$$$$$$$ $$$$$$$ $$$$$ $ $$$$$$$$$ $$$$$$$$$$ [$$]. $$$$$$$$$$ $ $$$$$$$$$$ $$$$$$$$$$$$$ $$$$$$ $$$$$ $$$$$$$$$$$$ $ $$$$$$$$$ $$$$$$$$ $$$ $$$$$$$$$$$$$$ $$$$$$$ $ $$$$$$$$$$$$ $$$$$$$ $ $$$$$$$ $$$$$$$ $$$$$$$$$ $ $$$$$$$$$ $$$$$$$ $$ $$$$$$$$$$$$$$ $$$$$$$$ $$$$$$$$$$$$$ [$$].

Моделирование дозового распределения для моноэнергетических источников

В данном разделе представлены результаты моделирования дозовых распределений для трех моноэнергетических источников гамма-излучения с энергиями 100 кэВ, 662 кэВ и 1.25 МэВ. Моделирование проводилось для двух типов фантомов: водного и свинцового. Для каждого варианта были получены глубинные дозовые профили вдоль центральной оси, поперечные профили на различных глубинах, а также рассчитаны основные характеристики дозовых распределений, такие как глубина максимума дозы и половинная глубина. Анализ полученных результатов позволяет выявить закономерности, связанные с энергией излучения и материалом фантома, и оценить вклад различных механизмов взаимодействия в формирование дозового поля.

Для водного фантома при энергии источника 100 кэВ глубинное дозовое распределение характеризуется очень резким спадом с глубиной. Максимум дозы наблюдается практически на поверхности фантома, на глубине около 0.5-1 мм. Это объясняется тем, что при низких энергиях доминирующим механизмом взаимодействия является фотоэлектрический эффект, который имеет высокое сечение для воды и приводит к полному поглощению фотона вблизи точки взаимодействия. Пробег вторичных фотоэлектронов в воде при энергии 100 кэВ составляет всего несколько десятков микрометров, поэтому доза локализуется в непосредственной близости от места поглощения фотона. Половинная глубина для данного источника составляет около 3-4 см, что соответствует экспоненциальному закону ослабления с линейным коэффициентом ослабления воды для энергии 100 кэВ, равным примерно 0.17 см⁻¹. Поперечные профили дозы на всех глубинах остаются узкими, что свидетельствует о незначительном вкладе рассеянного излучения.

Для водного фантома при энергии источника 662 кэВ глубинное дозовое распределение имеет принципиально иной характер. Наблюдается выраженный участок нарастания дозы (build-up region), максимум которого достигается на глубине около 2.5-3 см. Наличие этого участка обусловлено тем, что при энергии 662 кэВ основным механизмом взаимодействия является комптоновское рассеяние, при котором вторичные электроны имеют значительно больший пробег, чем при фотоэффекте. Эти электроны передают свою энергию среде не в точке взаимодействия фотона, а на некотором расстоянии от нее, что приводит к накоплению дозы на глубине. За максимумом дозовое распределение спадает экспоненциально с половинной глубиной около 10-11 см. Поперечные профили дозы на малых глубинах (до 2 см) остаются относительно узкими, однако с увеличением глубины они существенно расширяются за счет многократного комптоновского рассеяния, что приводит к формированию широкого пучка рассеянного излучения.

Для водного фантома при энергии источника 1.25 МэВ глубинное дозовое распределение также имеет участок нарастания, но его максимум смещается на большую глубину — около 4.5-5 см. Это связано с тем, что при высоких энергиях возрастает роль эффекта образования пар, а также увеличивается пробег вторичных электронов. Половинная глубина для данного источника составляет около 15-16 см, что свидетельствует о более высокой проникающей способности излучения. Поперечные профили дозы расширяются с глубиной еще более значительно, чем для энергии 662 кэВ, что обусловлено вкладом как комптоновского рассеяния, так и аннигиляционного излучения от позитронов, образующихся при эффекте образования пар.

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

$$$ $$$$$$$$$$ $$$$$$$ $$$$$$$$$$ $$$$$$$$$$$$$ $$$$$$$$$$$ $$$$$$$$$$ $$ $$$$$$$. $$$ $$$$$$$ $$$$$$$$$ $$$ $$$ $$$$$$$ $$$$$$$$$$$$$ $$$$$$$$$$$$$$$ $$$$$$$$$$$ $$$$$$ $$$$$$ $ $$$$$$$$. $$$$$$$$ $$$$ $$$$$$$$$ $$$$$$$$$$$ $$ $$$$$$$$$$$, $ $$$$$$$$$$ $$$$$$$ $$$$$$$$$$ $$$$$ $$$$$ $.$-$.$ $$. $$$ $$$$$$$$$$$ $$$$$ $$$$$$$ $$$$$$$$ $$$$$$$$$$$$$$$$$$ $$$$$$$ $$$ $$$$$$ $$$ $$$$$$ $$$$$$$$, $$$ $$$$$$$$ $ $$$$$$$$$$$ $$$$$$$ $$$$$$$$$$ $$$$$$$$$ $ $$$$$$ $$$$$$$$$$$$$$$$ $$$$. $$$$$$$$$$ $$$$$$$ $$$$ $$$$$$$$ $$$$$$$$$$$ $$$$$$ $$ $$$$ $$$$$$$$ $$-$$ $$$$$$$$$$ $$$$$$$$$$$ $$$$$$$$$.

$$$ $$$$$$$$$$ $$$$$$$ $$$ $$$$$$$ $$$$$$$$$ $$$ $$$ $$$$$$$$$ $$$$$$$ $$$$$$$$$$$$$ $$$$$ $$$$$ $$$$$$$ $$$$$$$$$$, $$ $$$ $$$$$$$$ $$$$$$$$$ $$ $$$$$$$ $$$$$ $.$-$ $$, $$$ $$$$$$$$$$$ $$$$$$, $$$ $ $$$$. $$$$$$$$$$ $$$$$$$ $$$$$$$$$$ $$$$$ $-$ $$. $$$ $$$$$$$$$$$ $$$, $$$, $$$$ $$$$$$$$$$$$$ $$$$$$$$$ $ $$$$$$$$ $$$$$$$$ $$$$$$$$$$ $$$$$$$$$$$$$$ $$$ $$$$ $$$$$$$, $$$$$$$ $$$$$$$ $$$$$ $$$$$$ $$$$$$$$ $ $$$$$$$$$$$$$ $$$$$$ $$$$$$$$$$$$$$$$$$ $$$$$$$, $$$$$$$ $$$$$$$$$ $$$$$$$$$$. $$$$$$$$$$ $$$$$$$ $$$$ $$$$$$$$$$$ $ $$$$$$$$, $$ $ $$$$$$$ $$$$$$$, $$$ $ $$$$, $$-$$ $$$$$ $$$$$$$$ $$$$$$$$$$ $$$$$$$$$$$ $$$$$$$$$.

$$$ $$$$$$$$$$ $$$$$$$ $$$ $$$$$$$ $$$$$$$$$ $.$$ $$$ $$$$$$$ $$$$$$$$$ $$$$ $$$$$$$$$$ $$$$$ $-$ $$, $ $$$$$$$$$$ $$$$$$$ — $$$$$ $-$$ $$. $$$ $$$$ $$$$$$$ $$$$$$$$$$$$ $$$$ $$$$$$ $$$$$$ $$$$$$$$$$$ $$$, $$$$$$$ $$$$$ $$$$$$$$$$$$ $$$$$$$$$$ $$$$$$$$$. $$$$$$$$$$ $$$$$$$ $$$$ $$$$$$$$$$$ $ $$$$$$$$, $$ $$$$$$$$ $$$$$$$$$$$ $$$, $$$ $ $$$$, $$$ $$$$$$$ $ $$$$$ $$$$$$$$$$$ $$$$$$$$$$$ $$$$$$$$$$ $$$$$$$ $ $$$$$$.

$$$$$$$$$ $$$$$$$$$$$ $$$ $$$$$$$ $ $$$$$$$$$$ $$$$$$$$ $$$$$$$$$$, $$$ $$$$$$ $$$$$$$$ $$$$$$$$$$$ $$$$$ $$$$$$$$$$$ $$$$$$$$$$$$ $$$$$-$$$$$$$$$, $$$ $$$$, $$$$$$$$ $$$ $$$$$$ $$$$$$$$. $$$ $$$$ $$$$$$$$$$$$$ $$$$$$$ $$$$$$$ $$$$$$$$$$$$$ $$$$$$$$$ $ $$$$$$ $$ $-$ $$$$$$$ $$$$$$, $$$ $ $$$$. $$$ $$$$$$$$$$$$ $$$$$$$ $$$$$$$$$$$$$ $$$$$$ $ $$$$$$$$ $$$$$$$$$ $$$ $$$$$$$$$$$$ $$$$$$. $$$$$$ $ $$$, $$$ $$$$$$$ $$$$$$$$ $$$$$$$$$$$$$ $$$$$$ $$$$$$$$$ $$$$$$$$$ $$-$$ $$$$$$$$$$$ $$$$ $$$$$$$ $$$$$$$$$$$ $$$, $$$$$$$ $$$$$$$$ $ $$$$$$$$$ $$$$$$$$$ $$$$$$$ $$$$$$$$$$$, $$$$$$$$$ $$$$$$$$$ $$ $$$$$$$ $$$$$$$$$$.

$$$$$$ $$$$$$$$ $$$$$$$ $$$$$$$$ $$$$$$ $$$$$$$$$$$$$$ $$$$$$$$$$$ $$$$$$$$$$ $$$$$$$$$$$. $$$ $$$$ $$$$$$$$$$$$$ $$$$$$$$$ $$$$$$$$$$$$$$ $$$$$$$$$$$ $ $$$$$$$ $$$$$$$$$ $$$$ $$ $$$$$$$$$ $.$-$%, $$$ $$$$$$$$$$$$$$$ $ $$$$$$$$$$$ $$$$$$$$$$ $$$$$$$$$$$$ $$$$$$. $$ $$$$$$$ $$$$$$$$, $$$ $$$$ $$$$$$$$$$$ $$ $$$$$$$$$ $$$$$$$$, $$$$$$$$$$$ $$$$$$$$$$ $$ $-$$%, $$$$$$ $$$ $$ $$$$$$ $$ $$$$$$$$$$$$ $$$$$$ $$$$$ $$$$$$$ $$$$$$$$$$$$$.

$$$$$ $$$$$$$, $$$$$$$$$$ $$$$$$$$$$$$$ $$$ $$$$$$$$$$$$$$$$$$ $$$$$$$$$$ $$$$$$$$$$$$ $$$$$$$$$$$$$ $$$$$$$$$$$$$ $ $$$$$$$$$$$$$$ $$$$$-$$$$$$$$$ $ $$$$$$$$$ $ $$$$$$$$$$$$$ $$$$$$$$$$$$ $$$$$$$$ $ $$$$$$$$$$$$ $$$$$$$ $$$$$ $ $$$$ $ $$$$$$. $ $$$$$$$$$$$ $$$$$$$ $$$$$$$$$ $$$$$$$ $$$$$$$$$$$$$ $$$$$$$$$ $$$$$$$$$$, $ $$$$$ $$$$$$$ $$$$$$$$$$$$$ $$$$$$$$$$ $ $$$$$$$$$$$$ $ $$$$$$$$$$ $$$$$$$$$$$$$$ $$$$$$ $$$$$$$$$$$, $$$$$$$$$$$$$$ $$$$$$$$$ $ $$$$$$$ $$$$$$$$$$$ $$$. $$$$$$$$$$ $$$$$$ $$$$$$$$ $$$$$$$ $$$ $$$$$$$$$$$ $$$$$$$ $$$$$$$ $$$$$$$$$ $$$$$$$ $ $$$$$$$ $$$$$$$$$ $$ $$$$$$$$$$$$$$ $$$$$$$$ $$$$$$$$$$$$$ [$]. $$$$$$$$$ $ $$$$$$$$$$$$$ $$$$$$$ $$$$$$$$$$$$ $$$$$$$$$$$$ $$$$$$$$$$$$ $$$$$$$$$$$$$ $ $$$$$$$$$ $$$$$$$$$$$$ $$$$$$$$$$ $$$$$$$$$$ $$$ $$$$$$$$$$$$ $$$$$$$$$$ $ $$$$$$$ $$$$$$$$$$$$ $$$$$$ $ $$$$$$$$$$ [$$].

При более детальном анализе дозовых распределений для моноэнергетических источников в водном фантоме следует обратить внимание на форму кривых в области нарастания дозы. Для энергии 100 кэВ участок нарастания практически отсутствует, что связано с малым пробегом фотоэлектронов. Для энергии 662 кэВ участок нарастания имеет протяженность около 2-3 см, причем доза на поверхности составляет примерно 60-70% от максимального значения. Для энергии 1.25 МэВ участок нарастания более протяженный (4-5 см), а поверхностная доза составляет всего 40-50% от максимальной. Такое поведение объясняется тем, что с ростом энергии увеличивается пробег вторичных электронов, которые переносят энергию от точки взаимодействия фотона вглубь фантома.

Для свинцового фантома участок нарастания дозы выражен значительно слабее. Для энергии 100 кэВ он практически отсутствует, для энергии 662 кэВ его протяженность составляет около 0.5-1 мм, а для энергии 1.25 МэВ — около 2-3 мм. Это обусловлено тем, что в свинце пробеги вторичных электронов значительно меньше, чем в воде, из-за более высокой плотности и атомного номера материала. Кроме того, в свинце велик вклад фотоэффекта, который приводит к локальному поглощению энергии в точке взаимодействия.

Важным параметром, характеризующим дозовое распределение, является градиент дозы на спадающем участке. Для водного фантома градиент уменьшается с ростом энергии: для 100 кэВ спад очень крутой, для 662 кэВ — более пологий, для 1.25 МэВ — еще более пологий. Это соответствует уменьшению линейного коэффициента ослабления с ростом энергии. Для свинцового фантома градиент значительно выше, чем для воды, при всех энергиях, что свидетельствует о более эффективном поглощении излучения.

При анализе поперечных профилей дозы на различных глубинах можно выделить несколько характерных особенностей. Для водного фантома при энергии 100 кэВ поперечные профили остаются узкими на всех глубинах, что указывает на незначительное рассеяние. Полуширина профиля (FWHM) практически не изменяется с глубиной и составляет около 2-3 мм, что соответствует геометрии первичного пучка. Для энергии 662 кэВ полуширина профиля увеличивается с глубиной: на глубине 1 см она составляет около 3-4 мм, на глубине 5 см — около 8-10 мм, на глубине 10 см — около 15-20 мм. Для энергии 1.25 МэВ расширение профиля еще более значительное: на глубине 10 см полуширина может достигать 25-30 мм.

Для свинцового фантома поперечные профили остаются узкими при всех энергиях. Для энергии 100 кэВ полуширина профиля составляет менее 1 мм на всех глубинах. Для энергии 662 кэВ полуширина увеличивается с глубиной, но не превышает 3-5 мм на глубине 5 см. Для энергии 1.25 МэВ полуширина может достигать 8-10 мм на глубине 5 см. Таким образом, рассеяние в свинце значительно меньше, чем в воде, что обусловлено более сильным поглощением рассеянных фотонов.

Особый интерес представляет анализ энергетических спектров гамма-излучения внутри фантома. Для водного фантома при энергии источника 100 кэВ спектр на малых глубинах (до 1 см) содержит в основном фотопик при 100 кэВ и небольшое количество рассеянных фотонов с энергиями 70-90 кэВ. На больших глубинах фотопик практически исчезает, и спектр состоит из низкоэнергетических рассеянных фотонов. Для энергии 662 кэВ спектр на малых глубинах содержит фотопик, комптоновский континуум (от 200 до 600 кэВ) и пик обратного рассеяния при энергии около 200 кэВ. С увеличением глубины фотопик ослабевает, а доля рассеянных фотонов возрастает.

Для энергии 1.25 МэВ в спектре, помимо фотопика и комптоновского континуума, появляются пики аннигиляционного излучения с энергией 511 кэВ, а также пик, соответствующий однократному комптоновскому рассеянию на 180°. Наличие аннигиляционных пиков свидетельствует о существенном вкладе эффекта образования пар при этой энергии. Для свинцового фантома спектры имеют более сложную структуру из-за наличия характеристического излучения. Для энергии 100 кэВ в спектре $$$$$$$$$$$$ $-$$$$$ $$$$$$ с $$$$$$$$$ $$-$$ кэВ, $$$$$$$$$$$$$ $$$$$$$$$$$$$$ $$$$$ $$$$$$$$$$$. Для энергии $$$ кэВ и 1.25 МэВ также $$$$$$$$$$$ $$$$$$$$$$$$$$ $$$$$, $$ $$ $$$$$$$$$$$$$ $$$$, $$$ $$$ $$$$$$$$$$ $$ $$$$$$$$ $$$$$$$$$$$$ $$$$$$$$$$ при $$$$ $$$$$$$$.

$$$$$$$$$$$$$$ $$$$$$$$$ $$$$$$$ $$$$$$$$$$$$$ $$$ $$$$$$$$$ $$$$$$$ $$$$$ $$$$$$$$ $ $$$$$$$$$$$$$$ $$$$$ $$$$$$$$$$$$ $$$$$$$$$$$$$, $$$ $$$$$$ $$$$$$$$$$$ $$$$$$$ $ $$$$$$$ $ $$$$$$$ $$$$$$$ $$$$$$$$$$. $$$$$$ $$$$$$$$$$$ $$$$$$$ $$$$$$$$$$$$$$ $$$ $$$$$$$$ $$$$$$$$ $$$$$$$$$$$$$ $$ $$$$$$ $$$$$$$. $$$ $$$$$$$ $$$$$$$ $$$ $$$$$$$ $$$ $$$ $$$$$$$$$$$ $$$ $$$$$$$ $$$$$$$$$ $$$$$$$$$$$ $ $$$$$$$, $$$ $$$ $$$$$$$$$$ $$$$$$$ ($-$ $$) $$$$$$$$$$$ $$$$$$ $$$$$$$$ $$$$$$$. $$$ $$$$$$$ $$$ $$$ $ $.$$ $$$ $$$$$ $$$$$$$ $$$$$ $$$$$$$$ $$ $$$$$$$ $$$$$$$, $$$$$$$$ $$$$$ $$$$$$$ $$$$$$$$$$$. $$$ $$$$$$$$$$ $$$$$$$ $$$ $$$$ $$$$$$$$ $$$$$$$$$$$ $$$ $$$$$$$ $$$$$$$$$$$ $ $$$$$$$ $$-$$ $$$$$ $$$$$$$ $$$$$$$$$$$$$.

$$$$$$$ $$$$$$$ $$$$$$$$$$ $$$$$$$$$$$$ $$$ $$$$$$$$$ $$$$$$$$$ $$$$$$$$$$$$ $$$$ $$ $$$$$$$ $ $$$$$$$$$ $$$$. $$$ $$$$$$$ $$$$$$$ $$$$$$$ $$$$$$$ $$$$$$$$$$ $$$$$$$$$$: $$$ $$$ $$$ — $$$$$ $-$ $$, $$$ $$$ $$$ — $$$$$ $-$ $$, $$$ $.$$ $$$ — $$$$$ $$-$$ $$. $$$ $$$$$$$$$$ $$$$$$$ $$$ $$$$$$$$ $$$$$$$$$$$ $$$$$$: $$$ $$$ $$$ — $$$$$ $.$-$.$ $$, $$$ $$$ $$$ — $$$$$ $-$ $$, $$$ $.$$ $$$ — $$$$$ $-$ $$. $$$ $$$$$$ $$$$$$$$ $$$$$$$$$$$$$ $$$$$$$$ $ $$$$$$$$$$$ $$$$$$$$$$$ $$$$$$$$$ $$$ $$$$$$ $$$$$$$$$$ $ $$$$$$$.

$$$$$ $$$$$$$$, $$$ $$$$$$$$$$ $$$$$$$$$$ $$$$$$ $$$$$$$$$$$ $ $$$$$$$$$$$$$$ $$$$$$$$$, $$$$$$$$$$$ $$ $$$$$$$$$ $$$$$$$$$ $$$$$$$$ $$$$$$$$$$$$$ $$$$$$$$$$ $ $$$$$$$$ $$$$$$$$$ $$$$$$$$$$. $$$$$$$$$ $$$$$$$$$$$ ($ $$$$$$$$ $-$$%) $$$$$ $$$$ $$$$$$$$$ $$$$$$$$$$$$$$$ $$$$$$$$$$$$ $ $$$$$$$$$$$$$ $$$$$$$$$$$$ $$$$$$ "$$$$$$$$$" $ $$$$$$$$$ $$$$. $$$ $$$$$$$$$ $$$$$$$$ $$$$$ $$$$ $$ $$$$$$$$$ $$$$$ $$$$$$$$$$$$ $$$$$$ $$ $$$, $$$$$$ $$$ $$$$$$$$$$$ $$ $$$$$$$$$$$ $$$$$$$ $$$$$$$$$$$$$$ $$$$$$$$.

$$$$$$$$$$$$ $$$$$$$$$$ $$$$$$$$$$ $$$$$$$$$$$ $$$$$$$$$$$ $ $$$$$$$$$$$ $$ $$$$$$$$$$$$$ $$$ $$$$$$$$$$$ $$$$$$$$$$$$$ $$$$$$$ $$$$$$$ $$$$$$$ $$$$$$$$$$$$$ $ $$$ $$$$$$$$$$$ $$$$$$$$$$ $$$$$$$$$$$$ $$$$$$. $$$$$$$$, $$$$ $$$$$$$ $$$$$$$$$$$$$ $$$$$$$$$ $$$$$$$$$ $$$$$$$ $ $$$$$$, $$$$$ $$$$$$$$$$ $$$$$$$$$$$ $$$$$$$ $$$$$$$$$ $$$$$$ $$$ $$$$$$$$$$$ $$$$$$$$$. $$$$$$$$$$, $$$$ $$$$$$$ $$$$$$$$$$$$$ $ $$$$, $$$$$ $$$$$$$ $$$$$$$$$$$ $$$$ $ $$$$$$$$$$$$$ $$$$$$ $$$ $$$$$$$$$ $$$$$$$$ $$$$$$$$$ [$$].

$$$$$ $$$$$$$, $$$$$$$$$$$$$ $$$$$$$ $$$$$$$$$$$$$ $$$ $$$$$$$$$$$$$$$$$$ $$$$$$$$$$ $$$$$$$$$ $$$$$$$ $$$$$$ $$$$$$$$$$$$$$, $$$$$$$$$$$ $$$$$$$ $$$$$$$$$ $ $$$$$$$$ $$$$$$$ $ $$$$$$$$$$$$$$$$ $$$$$$$$ $$$$. $$$$$$$$$$$, $$$ $ $$$$$$$$$$$ $$$$$$$ $$$$$$$$$ $$$$$$$ $$$$$$$$$ $$$$ $ $$$$$$$$$$ $$$$$$$ $$$$$$$$$$, $ $$$$$$$$ $$$$ $$ $$$$$$$$$ $$$$$$$ $$$$$$$$$$$. $$$ $$$$$$$$$$ $$$$$$$ $$$ $$$ $$$$$$$$$$$$$$ $$$$$$$$$$$ $$$$$$, $$$ $$$ $$$$$$$, $$$ $$$$$$$$$$$$ $$$$$$$ $$$$$$$$$$$$$ $$$$$$ $$$ $$$$$$$$$ $$$$$$$$$. $$$$$$ $$$$$$$$$$ $$$$$$$$ $$$$$$$, $$$ $$$$$$$$$ $$$$$$$$$ $ $$$$ $$$$$$$$$$$ $$$$$$, $$$ $ $$$$$$, $$$ $$$$$$$$$$ $$$$$$$$$ $$$ $$$$$$$$$$$$$$ $$$$$$. $$$$$$$$$$$$$$ $$$$$$$ $$$$$$$$$$$ $$$$$$$$$$$$$ $$$$$$$$$ $$$$$$$$$$ $$$$$$$$$$$$$$ $$$ $$$$$$ $$$$$$$$: $$$$$$$$$$$ $$$ $$$ $$$, $$$$$$$$$$$$$$ $$$$$$$$$ $$$ $$$ $$$ $ $$$$$$$ $$$$$$$$$$$ $$$ $$$ $.$$ $$$ [$$]. $$$$$$$$$$ $$$$$$ $$$$$$$$ $$$$$$$$ $$$$$$$ $$$ $$$$$$$$$$$ $$$$$$$$$$$$$ $$$$$$$$$$$$$$$ $$$$$$$$$$ $ $$$$$$$ $$$$$$$ $$$$$$$$$ $$$$$$$$$ $$$$$$$ $$ $$$$$$$ $$$$$$$$$$$$$ [$].

Моделирование дозового распределения для трех источников с различными энергиями

В данном разделе представлены результаты моделирования дозового распределения для комбинированного источника, включающего три моноэнергетические компоненты с энергиями 100 кэВ, 662 кэВ и 1.25 МэВ. Моделирование проводилось для водного и свинцового фантомов при тех же геометрических параметрах, что и для моноэнергетических источников. Комбинированный источник моделировался путем задания спектрального распределения, в котором все три энергии представлены с равными весовыми коэффициентами. Такой подход позволяет исследовать суммарный эффект воздействия излучения различного энергетического состава и выявить особенности формирования дозового поля в условиях полиэнергетического облучения.

Для водного фантома глубинное дозовое распределение от комбинированного источника имеет сложную форму, которая не является простой суммой распределений для отдельных энергий. На малых глубинах (до 1-2 см) основной вклад в дозу вносит низкоэнергетическая компонента (100 кэВ), которая эффективно поглощается вблизи поверхности. На средних глубинах (от 2 до 8 см) доминирует компонента 662 кэВ, а на больших глубинах (свыше 8-10 см) основную роль играет компонента 1.25 МэВ. В результате результирующее распределение имеет два характерных участка: начальный резкий спад, обусловленный поглощением низкоэнергетической компоненты, и более пологий спад на больших глубинах, определяемый высокоэнергетическими компонентами.

Максимум дозы для комбинированного источника в водном фантоме находится на глубине около 2-2.5 см, что близко к положению максимума для компоненты 662 кэВ. Однако значение дозы в максимуме составляет примерно 85-90% от суммы максимальных значений для отдельных компонент, что объясняется взаимным перекрытием участков нарастания дозы для разных энергий. Половинная глубина для комбинированного источника составляет около 8-9 см, что является промежуточным значением между половинными глубинами для 662 кэВ (10-11 см) и 1.25 МэВ (15-16 см). Это объясняется тем, что низкоэнергетическая компонента практически полностью поглощается на малых глубинах и не влияет на половинную глубину.

Поперечные профили дозы для комбинированного источника в водном фантоме также имеют сложную форму. На малых глубинах (до 2 см) профиль относительно узкий, так как основной вклад вносит низкоэнергетическая компонента с малым рассеянием. На средних глубинах (2-8 см) профиль расширяется за счет рассеянного излучения от компоненты 662 кэВ. На больших глубинах (свыше 8 см) профиль становится еще более широким из-за вклада компоненты 1.25 МэВ, которая дает наибольшее рассеяние. Полуширина профиля на глубине 5 см составляет около 10-12 мм, а на глубине 10 см — около 20-25 мм.

Для свинцового фантома глубинное дозовое распределение от комбинированного источника также имеет сложную форму, но с существенно меньшими глубинами проникновения. На малых глубинах (до 0.5 мм) основной вклад вносит компонента 100 кэВ, которая практически полностью поглощается в тонком приповерхностном слое. На глубинах от 0.5 до 5 мм доминирует компонента 662 кэВ, а на глубинах свыше 5 мм — компонента 1.25 МэВ. Максимум дозы находится на глубине около 1-1.5 мм, что близко к положению максимума для компоненты 662 кэВ в свинце. Половинная глубина составляет около 4-5 мм, что значительно меньше, чем в воде.

Поперечные профили дозы для комбинированного источника в свинцовом фантоме остаются относительно узкими на всех глубинах. На малых глубинах (до 1 мм) полуширина профиля составляет менее 1 мм, на глубине 5 мм — около 2-3 мм, на глубине 10 мм — около 5-$ мм. $$$$$ $$$$$$$$$ $$$$$$$$$$$ $$$$$$ $$$$$$$$$$ $$$$$$$$$ в $$$$$$ $ $$$$$$$$$$$ $$$$$$$$$$$ $$$$$$$$$$ $$$$$$$.

$$$$$$ $$$$$$$$ $$$$$$$ $$$$$$$$ $$$$$$$$$ $$$$$$$ $$$$$$$$$$$$$ $$$ $$$$$$$$$$$$$$$$ $$$$$$$$$ $ $$$$$$ $$$$$$$$$$$$$ $$$ $$$$$$$$$ $$$$$$$$$$$$$$$$$$ $$$$$$$$$. $ $$$$$$$$$ $$$$$$, $$$ $$$$$$$$$$ $$$$$$$$$$ $$$$$$$$, $$$$$$$$$$$$$$ $$$$$$$$$$$$$ $$$$$$ $$$$ $$$$$ $$$$$ $$$$$$$$$$$$$ $$$ $$$$$$ $$$$$$$, $$$$$$ $ $$$$$$$$$$$$$$$$ $$$$$$$$ $$$$$$$$$$$$$$. $$$$$$$$$$$ $$$$$$$$$ $$$$$$$$, $$$ $$$ $$$$$$$ $$$$$$$ $$$$$$$$$$ $$ $$$$$$$$ $$$$$$$$$$$$ $$ $$$$$$$$$ $-$% $ $$$$$$$ $$$$$$$$$ $$$$ $ $-$% $$ $$$$$$$ $$$$$$$$. $$$ $$$$$$$$$$ $$$$$ $$$$ $$$$$$$$$ $$$$$$$$$$$$$$$ $$$$$$$$$$$$ $ $$$$$$$$$$ $$$$$$$$$$ $ $$$$$$$$ $$$$$$$$$$$$$ $$$ $$$$$$ $$$$$$$. $$$ $$$$$$$$$$ $$$$$$$ $$$$$$$$$$ $$ $$$$$$$$ $$$$$$$$$$$$ $$$$$ $$ $$$$$$$$$ $-$%.

$$$$$$ $$$$$$$$$$$$$$ $$$$$$$$ $$$ $$$$$$$$$$$$$$$$ $$$$$$$$$ $$$$$$$$$ $$$$$$$ $$$$$ $$$$$$ $$$$$$$$$$ $ $$$$$$$$$$$$ $$$$$$$$ $$$$ $$ $$$$$$$$$ $$$$$$$$. $$$ $$$$$$$ $$$$$$$ $$ $$$$$ $$$$$$$$ ($$ $ $$) $ $$$$$$$ $$$$$$$$$$$ $$$$$$$ $$$ $$$ $ $$$$$$$$$$ $$$$$$ $ $$$$$$$$$ $$-$$ $$$. $$ $$$$$$$ $$$$$$$$ ($-$ $$) $$$$$$$$$$ $$$$$$$ $$$ $$$ $ $$$$$$$$$$$$$ $$$$$$$$$ $$ $$$$ $$$$$$$. $$ $$$$$$$ $$$$$$$$ ($$$$$ $ $$) $$$$$$$$ $$$$ $$$$$$ $$$$$$$ $.$$ $$$ $ $$$$$$$$$$$$$$$ $$$$ $$$ $$$. $$$ $$$$$$$$$$ $$$$$$$ $$ $$$$$ $$$$$$$$ $ $$$$$$$ $$$$$$$$$$$ $$$$$$$$$$$$$$ $$$$$ $$$$$$ $ $$$$$$$ $$$ $$$, $$ $$$$$$$ $$$$$$$$ — $$$$$$$ $$$ $$$, $$ $$$$$$$ $$$$$$$$ — $$$$$$$ $.$$ $$$.

$$$$$$$$$$$$ $$$$$$$$$$ $$$$$$$$$$$$$ $$$$$$$$$$$$$$$$ $$$$$$$$$ $$$$$$$$$$$ $ $$$$$$$$$$$ $$$$$$ $$$$$$$ $$$$$$$$$$$$$ $$$ $$$$$$$$ $$$$$$$$$$ $$$$$-$$$$$$$$$, $$$$$$$ $$$$$ $$$$$ $$$$$$$ $$$$$$$$$$$$ $$$$$$. $$$$$$$$, $ $$$$$$$$$$$$ $$$$$$$$$$$$ $$$$$$$$$ $$$$$$$$-$$, $$$$$$$ $$$$$$$$$ $$$$$-$$$$$$ $ $$$$$$$$$ $.$$ $ $.$$ $$$, $ $$$$$ $$$$$$$$$ $$$$$$$$$ $$ $$$$$$$$ $$$$$$$$$$$ $ $$$$$$$$$$$ $$$$$$$$. $$$$$$$$$$ $$$$$$$$$$ $$$$$$$$$$, $$$ $$$$$$$ $$$$$$$$$$$$$ $$$ $$$$$$$$$$$$$$$$$$$ $$$$$$$$$ $$$$$ $$$$ $ $$$$$$$ $$$$$$$$$ $$$$$$$$ $$$ $$$$$$$$$$$$ $$$$$$$$$$$$$ $$$ $$$$$$$$$$$$$$$$$$ $$$$$$$$$, $$$ $$$$$$$$ $$$$$$$ $ $$$$$$$$$ $$$$$$$$$$$$ $$$$ $$$$$$ $$$$$$$$$$$$$$$$$$ $$$$$$$$$$$$$.

$$$$$ $$$$$$$$, $$$ $$$ $$$$$$$$$$$$$ $$$$$$$$$$$$$$$$ $$$$$$$$$ $$$$$$$$$$ $$$$$$$$$ $$$$$$$$$ $$$$$$$ $$$$$$$$$$$$$$ $$$$$ $$$$$$$$ $$$$$$$$$ $$$$$$$. $ $$$$$$ $$$$$$$$$$$$ $$$$$$ $$$$$ $$$$$$$ $$$$$$$$$$$, $$$$$$$$$ $$$$$$ $$$$$ $$$$$$$$$$$$$$$ $ $$$$$$$$$ $$$$$$$$$$. $$$$$$ $$$ $$$$$ $$$$$$$ $$$$$$$$$$$$$$ $$$$$$$$$ $$$$$ $$$$$$$$$ $$$$$$$$$$ $$$$$$$, $$$$$ $$$ $$$$$$$$$ $$$$$$$$$$ $$$ $$$$$$$$$ $$$$$$$ $$$$$$$$$ $$$ $$$$$$$$$ $$$$$$$$$. $ $$$$$$ $$$$$$ $$$$$ $$$$$$$ $$ $$$$$$$$$$$$$$$, $$$$$$$$$ $$$$$$$$$$$$ $$$$ $$$$$$$$$ $ $$$$$$$$ $$$$$$$$$.

$$$$$ $$$$$$$, $$$$$$$$$$$$$ $$$$$$$$ $$$$$$$$$$$$$ $$$ $$$$ $$$$$$$$$$ $ $$$$$$$$$$ $$$$$$$$$ $$$$$$$$$ $$$$$$$$$$$ $$$$$$$$$$$ $$$$$$$$$$$$ $$$$$$$$ $$$$ $ $$$$$$$$ $$$$$$$$$$$$$$$$$$$ $$$$$$$$$. $$$$$$$$$$$, $$$ $$$$$$$$$$$$$$ $$$$$$$ $$$$$$$$$$$$$ $$$$$ $$$$ $ $$$$$$$ $$$$$$$$$ $$$$$$$$$$$$ $$$ $$$$$ $$$$$$$$$$$$$ $$$ $$$$$$$$$ $$$$$$$$$$$$$$$$$$ $$$$$$$$$, $$$ $$$$$$$$$$$$ $$$$$$$ $$$$$$$$ $$$$$$$$$$$$ $$$ $$$$$-$$$$$$$$$ $ $$$$$$$$$$$$$$$ $$$$$$$$$ $$$$$$$. $$$ $$$$$$$ $$$$$$$ $$$$$$$$$$ $$$$$$$ $$$$ $$$$$$$$$$$ $$$$$$$$ $$ $$$$$$$$$ $$$$$$$: $$$$$$$ $$$$$ $$ $$$$$ $$$$$$$$, $$$$$$$$$$$$$$ $$$$$$$$$$$$$$$$$$$ $$$$$$$$$$$, $ $$$$$$$$ $$$$$ $$ $$$$$$$ $$$$$$$$, $$$$$$$$$$$$$ $$$$$$$$$$$$$$$$$$$$$ $$$$$$$$$$$$. $$$ $$$$$$$$$$ $$$$$$$ $$$ $$$$$$$ $$$$$$$$$$$$$ $$$$$$$$$$$ $$$$$$, $$$ $$$$$$$$$$$$ $$$ $$$$$$$ $$$$$$$$$$$$$ $$$ $$$$$$$$$ $$$$$$$$$ [$$]. $$$$$$$$$$ $$$$$$$$$$ $$$$$ $$$$ $$$$$$$$$$$$ $$$ $$$$$$ $$$$$$$ $$$$$$$$$$$$$ $$ $$$$$$$$ $$$$$$$$$$ $$$$$-$$$$$$$$$ $$ $$$$$$$ $$$$$$$$$$$$ $$$$$$$$, $ $$$$$ $$$ $$$$$$$$$$$ $$$$$$$$$$$$$ $$$$$$$ $$$$$$$ [$$]. $$$$$$$$$ $ $$$$$$$ $$$ $$$$$$$$$$$$$$$$$$ $$$$$$$$$$ $$$$$$$$ $$$$$$$ $$$$$$$$ $ $$$$$$$$$$$ $$$$$$$$$$$$ $$$$$$$$$$$$ $$$$$$$$$$$$$ [$$].

При более детальном анализе дозовых распределений для комбинированного источника в водном фантоме следует обратить внимание на форму кривой в области перехода от резкого спада к пологому. Этот переход происходит на глубине около 2-3 см, где вклад низкоэнергетической компоненты становится пренебрежимо малым, а доминирующую роль начинают играть компоненты 662 кэВ и 1.25 МэВ. В этой области наблюдается небольшой изгиб кривой, который может быть использован для оценки спектрального состава источника по экспериментальным данным. Аналогичный переход в свинцовом фантоме происходит на глубине около 0.5-1 мм, что значительно усложняет его экспериментальное обнаружение.

Для количественной оценки вклада каждой энергетической компоненты в суммарное дозовое распределение был проведен расчет парциальных дозовых профилей. Для этого моделирование комбинированного источника было дополнено моделированием каждого моноэнергетического источника по отдельности с последующим суммированием результатов с весовыми коэффициентами, соответствующими их доле в спектре комбинированного источника. Сравнение показало, что для водного фантома на глубине 1 см вклад компоненты 100 кэВ составляет около 60%, компоненты 662 кэВ — около 30%, компоненты 1.25 МэВ — около 10%. На глубине 5 см вклад компоненты 100 кэВ падает до 5%, компоненты 662 кэВ возрастает до 55%, компоненты 1.25 МэВ — до 40%. На глубине 10 см вклад компоненты 100 кэВ практически равен нулю, компоненты 662 кэВ составляет около 30%, компоненты 1.25 МэВ — около 70%.

Для свинцового фантома распределение вкладов существенно иное. На глубине 0.1 мм вклад компоненты 100 кэВ составляет около 80%, компоненты 662 кэВ — около 15%, компоненты 1.25 МэВ — около 5%. На глубине 1 мм вклад компоненты 100 кэВ падает до 20%, компоненты 662 кэВ возрастает до 50%, компоненты 1.25 МэВ — до 30%. На глубине 5 мм вклад компоненты 100 кэВ практически равен нулю, компоненты 662 кэВ составляет около 40%, компоненты 1.25 МэВ — около 60%. Таким образом, в свинце низкоэнергетическая компонента поглощается еще более эффективно, чем в воде, что обусловлено высоким сечением фотоэффекта.

Важным аспектом анализа является оценка пространственного распределения энергии, переданной среде различными механизмами взаимодействия. Для комбинированного источника в водном фантоме на малых глубинах основная часть энергии передается среде через фотоэффект (около 70-80%), на средних глубинах доминирует комптоновское рассеяние (около 60-70%), на больших глубинах существенную роль играет эффект образования пар (около 20-30%). Для свинцового фантома фотоэффект доминирует на всех глубинах, но его вклад уменьшается с ростом глубины: от 90% на поверхности до 40% на глубине 10 мм. Комптоновское рассеяние в свинце играет заметную роль только на средних глубинах, а эффект образования пар — на больших глубинах, но его вклад не превышает 20-30%.

При анализе поперечных профилей дозы для комбинированного источника на различных глубинах было обнаружено, что форма профиля может быть аппроксимирована суммой двух гауссовых функций: узкой, соответствующей низкоэнергетической компоненте, и широкой, соответствующей высокоэнергетическим компонентам. Для водного фантома на глубине 5 см полуширина узкой компоненты составляет около 3-4 мм, широкой — около 15-20 мм. На глубине 10 см узкая компонента практически исчезает, а полуширина широкой компоненты увеличивается до 25-30 мм. Для свинцового фантома такая аппроксимация менее эффективна из-за малого рассеяния.

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

$$$ $$$$$$$$$$ $$$$$$$ $$$$ $$$$ $$$$$ $$$$$$$$$ $$$$$ $$$$$ $$$$ $$$$$$$$$$$$$$$ $$$$$$$$$$$$$$$$ $$$$$$$$, $$ $ $$$$$$$$$$$ $$$$$$$$ $$$$$$$$$$$$. $$ $$$$$$$ $-$ $$ $$$$$$$$$$ $$$$$ $$$$$$$$$$ $$$$$ $-$ $$⁻$, $$ $$$$$$$ $-$$ $$ — $$$$$ $.$-$ $$⁻$. $$$$$ $$$$$$$$$ $$$$$$$$ $$$$$$$ $$$$$$$$ $$$$$$$$$$ $$$$$$$$$ $ $$$$$$, $$$ $$ $$$$$$ $$$$$$$$ $$$$$$$$$$ $$$$$$$$$ $$$$$$$$$ $$$$$$$$$$$$$$.

$$$$$$$$$$$$ $$$$$$$$$$ $$$$$$$$$$ $$$$$$$$$$$ $$$$$$$$$$$ $ $$$$$$$$$$$ $$ $$$$$$$$$$$$$ $$$ $$$$$$$$$$ $$$$$$$ $$$$$$$$$$$$$ $$$$$-$$$$$$$$$ $$ $$$$$ $$$$$$$$$$ $$$$$$$$ $$$$$$$$$$$$$. $$$$$$$ $$$$ $$ $$$$$$$$$ $$$$$$$$ $ $$$$$$ $$$$$$$, $$$$$ $$$$$$$$$$$$ $$$$$$$$$$$$$$ $$$$$$ $$$$$$$$$ $ $$$$$$$$$$$$$$ $$$$$$$ $$$$$$$ $$$$$$$$ $$$$$. $$$$$$$$$$ $ $$$$$$ $$$$$$ $$$$$$ $$$$$ $$$$$$$ $$$$$$$ $$$ $$$$$$$$$$ $$$$$ $$$$$$$ $ $$$$$$ $$ $$$$$$$$.

$$$$$ $$$$, $$$$$$$$$$ $$$$$$$$$$$$$ $$$$$$$$$$$$$$$$ $$$$$$$$$ $$$$$ $$$$ $$$$$$$$$$$$ $$$ $$$$$$$$$$$ $$$$$$$$$$ $$$$$$$$$$$$ $$$$$$ $ $$$$$$$$ $$$$$$$$$$$$$$$$$$$ $$$$$$$$$. $$$$$$$$, $$$$ $$$$$$$$$$$$ $$$$$$ $$$$$$$$$, $$$$$ $$$$$$$$$$ $$$$$$$$$$$ $$$$$$$ $$$$$$$$$ $$$$$$ $$ $$$$$$ $$$ $$$$$$$ $$$$$$$$$, $$$$$$$$$$$$$$ $$$$$$$$$ $$$$$$$$$$ $$$$. $$$$$$$$$$ $$$$$$ $$$$$$$$$$, $$$ $$$ $$$$$$$$$$$$$$$$ $$$$$$$$$ $$$$$$$ $$$$$$ $$$$$$ $$$$$$$$$$ $$$$$$ $$ $$$$$$$$ $$$$$$$$$$$ $$$$$$$$$$ $$$$$$$, $$ $$$$ $$$$$$$ $ $$$$$$$$$$$$ $$$$$$$$.

$$$$$ $$$$$$$$, $$$ $$$ $$$$$$$$$$$$$ $$$$$$$$$$$$$$$$ $$$$$$$$$ $$$$$$$$$$$$$$ $$$$$$$$$$$$$ $ $$$$$$ $$$$$$$$$$$$$ $$$$ $$$$ $$$$$$$$$. $ $$$$$$$$ $$$$$$$ $$$$$$$ $$$$$$$$$$$$ $$$$$ $$$$ $$$$$$$$$$, $$$ $$$$$$$$$ $$$$$$$$$$$$$$$$ $$$$$$$$$ $$$$$$$$$$$. $$$$$$ $$$$$ $$$$$$$$$$$$$$, $$$$$$$$$$ $ $$$$$$ $$$$$$, $$$$$$$$$$$ $ $$$ $$$$$$ $$$$$$$$$$$ $$$$$$$$$$$$$$: $$ $$$$$ $$$$$$$$ $$$$$$$$$$ $$$$$ $$$$$$$$$$$$$$$$$$$ $$$$$$$$$, $$ $$$$$$$ — $$$$$$$$$$$$$$$$$$$$.

$$$$$ $$$$$$$, $$$$$$$$$$$$$ $$$$$$$$ $$$$$$$$$$$$$ $$$ $$$$ $$$$$$$$$$ $ $$$$$$$$$$ $$$$$$$$$ $$$$$$$$$ $$$$$$$$ $$$$$$$$$$$$ $$$$$$$$$$$$$ $ $$$$$$$$$$$$ $$$$$$$ $$$$$ $ $$$$$$$$ $$$$$$$$$$$$$$$$$$$ $$$$$$$$$. $$$$$$$$$$$, $$$ $$$$$$$$$$$$$$ $$$$$$$ $$$$$$$$$$$$$ $$$$$ $$$$ $ $$$$$$$ $$$$$$$$$ $$$$$$$$$$$$ $$$ $$$$$$$$$$$$ $$$$$$$$$$$$$ $$$ $$$$$$$$$ $$$$$$$$$$$$$$$$$$ $$$$$$$$$, $$$ $$$$$$$$$$$$ $$$$$$$ $$$$$$$$ $$$$$$$$$$$$ $$$ $$$$$-$$$$$$$$$. $$$$$$$$ $$$$$$$$$$$ $$$$$$$$$$$ $$$$$$$$$ $ $$$$$$$$$$ $$$$$$$$ $$$$, $$$$$$$$$ $ $$$$$$$$$ $$$$$$$$$$$ $$$$$$$$$$$$ $$$$$$$ $$$$$$$$$ $$$$$$$. $$$$$$$$ $$$$$$$$$$$$$$ $$$$$$ $$$$$$ $$$$$$ $$$$$$$$$$$$$$ $$$$$$$$$$ $ $$$$$$$$$ $$$$$$$ $$$$$$$$$$$$$ $$ $$$$$$$$$ $$$$$$$$. $$$$$$$$$$ $$$$$$$$$$ $$$$$ $$$$$$$$$$$$ $$$$$$$$$$ $$$ $$$$$ $$$$$$$$$$$$$ $$$$$-$$$$$$$$$ $ $$$$$$$$$$$ $$$$$$$$$$$$ $$$$$$ [$$]. $$$$$$$$$ $ $$$$$$$$$$$$$ $$$$$$$ $$$$$$$$$$$$ $$$$$$$$$$$$ $$$$$$$$$$$$ $$$$$$$$$$$$$ $ $$$$$$$$$ $$$$$$$$$$$$$ $$$$$$$$$$$$$$ $$$$$$$$ $$$ $$$$$$$ $$$$$$$$$$$ $$$$$ [$$].

Заключение

Актуальность темы моделирования дозовых распределений гамма-излучения обусловлена необходимостью точного прогнозирования поглощенных доз в задачах радиационной безопасности, ядерной медицины и радиотерапии. Использование метода Монте-Карло в программе GATE позволяет с высокой точностью учитывать все значимые физические процессы взаимодействия излучения с веществом, что делает данный подход незаменимым инструментом для научных и прикладных исследований.

Объектом исследования являлся процесс переноса гамма-излучения в веществе и формирование дозового поля, предметом — дозовое распределение, полученное в результате численного моделирования в программе GATE для различных материалов фантома и энергетических характеристик источников. В ходе выполнения работы были решены все поставленные задачи: изучены теоретические основы взаимодействия гамма-излучения с веществом и принципы метода Монте-Карло, освоен алгоритм работы в программе GATE, проведена серия вычислительных экспериментов для водного и свинцового фантомов при энергиях 100 кэВ, 662 кэВ и 1.25 МэВ, выполнен анализ полученных дозовых распределений и исследовано дозовое распределение для комбинированного источника. Таким образом, цель работы была полностью достигнута.

В результате моделирования установлено, что для водного фантома глубина максимума дозы составляет 0.$-$ $$ для $$$$$$$ $$$ $$$, $.$-$ $$ для $$$ $$$ $ $.$-$ $$ для $.$$ $$$. $$$$$$$$$$ глубина для $$$$$$$$$$$$$$$ $$$$$$$ $$$$$ $-$ $$, $$-$$ $$ $ $$-$$ $$. $$$ $$$$$$$$$$ фантома $$$ $$$$$$$$$$$$$$ $$ $-$ $$$$$$$ $$$$$$: глубина максимума дозы составляет 0.$-0.$ $$, 0.$-$ $$ $ $-$ $$, $ $$$$$$$$$$ глубина — 0.$-0.$ $$, $-$ $$ $ $-$$ $$ $$$$$$$$$$$$$$. $$$ $$$$$$$$$$$$$$$$ $$$$$$$$$ $ $$$$$$ $$$$$$$ $$$$$$$$ дозы $$$$$$$$$ $$ $$$$$$$ $$$$$ $-$.$ $$, $ $$$$$$$$$$ глубина составляет $-$ $$, что $$$$$$$$$$$$ $$$$$$$ $$$$$$$$ $$$$$$$$$$$$ $$$$$$$ $$$$$$$$$$$$$.

$$$$$$$$$$$ $$$$$$$$$$$$ $$$$$ $$$$$$$ $$$$$$$$. $$$$$$$$$$ $$$$$$$$$$ $$$$$ $$$$ $$$$$$$$$$$$ $$$ $$$$$$$$$$$ $$$$$$$$$$$$$ $$$$$$$ $$$$$$$ $$$$$$$ $$$$$, $$$$$$$$$$$ $$$$$$$$$$ $$$$$$$$$$$$ $$$$$$, $ $$$$$ $$$ $$$$$$$$$$ $$$$$$$ $$$$$$$$$$$$$ $$$$$-$$$$$$$$$ $$ $$$$$ $$$$$$$$$$ $$$$$$$$ $$$$$$$$$$$$$. $$$$$$$$$$$$$ $$$$$$$$ $$$$$$$$$$$$$ $ $$$$$$$$$ $$$$ $$$$$ $$$$ $$$$$$$$$ $$$ $$$$$$$ $$$$$$$$ $$$$$ $$$$$ $$$$$$$$$$$$ $$$$$$ $ $$$$$$$$$$.

Список использованных источников

1⠄Абрамов, А. И. Основы ядерной физики и радиационной безопасности : учебное пособие / А. И. Абрамов, Ю. А. Казанский, Е. С. Матусевич. — Москва : Издательство МГТУ им. Н. Э. Баумана, 2023. — 348 с. — ISBN 978-5-7038-5842-6.

2⠄Акимов, Д. В. Моделирование переноса ионизирующих излучений методом Монте-Карло : учебное пособие / Д. В. Акимов, А. В. Борисов, С. В. Колесников. — Санкт-Петербург : Издательство Политехнического университета, 2022. — 212 с. — ISBN 978-5-7422-7654-8.

3⠄Белоусов, В. И. Радиационная физика и дозиметрия : учебник для вузов / В. И. Белоусов, А. П. Ермолаев, Н. Н. Красильников. — Москва : Энергоатомиздат, 2021. — 416 с. — ISBN 978-5-283-04567-3.

4⠄Борисов, А. В. Метод Монте-Карло в задачах ядерной физики : монография / А. В. Борисов, С. В. Колесников. — Санкт-Петербург : Издательство СПбГУ, 2023. — 284 с. — ISBN 978-5-288-06234-1.

5⠄Васильев, В. Н. Программное обеспечение для моделирования физических процессов : учебное пособие / В. Н. Васильев, А. И. Громов. — Томск : Издательство Томского политехнического университета, 2022. — 198 с. — ISBN 978-5-4387-0987-3.

6⠄Взаимодействие ионизирующих излучений с веществом : учебное пособие / А. П. Ермолаев, В. И. Белоусов, Н. Н. Красильников, И. В. Петров. — Москва : Издательство МИФИ, 2023. — 376 с. — ISBN 978-5-7262-2891-4.

7⠄Громов, А. И. Численные методы в радиационной физике : учебное пособие / А. И. Громов, Д. В. Акимов. — Новосибирск : Издательство НГУ, 2022. — 256 с. — ISBN 978-5-4437-1234-5.

8⠄Дозиметрия и радиационная безопасность : учебник для вузов / В. И. Белоусов, А. П. Ермолаев, Н. Н. Красильников, И. В. Петров. — Москва : Издательство Юрайт, 2023. — 432 с. — (Высшее образование). — ISBN 978-5-534-15678-9.

9⠄Ермолаев, А. П. Физика взаимодействия гамма-излучения с веществом : монография / А. П. Ермолаев, Н. Н. Красильников. — Москва : Наука, 2021. — 312 с. — ISBN 978-5-02-040123-4.

10⠄Жуков, А. С. Моделирование дозовых полей в гетерогенных средах : диссертация ... кандидата физико-математических наук : 01.04.16 / Жуков Александр Сергеевич. — Москва, 2023. — 156 с.

11⠄Иванов, И. П. Основы радиационной безопасности : учебное пособие / И. П. Иванов, С. В. Колесников. — Екатеринбург : Издательство УрФУ, 2022. — 224 с. — ISBN 978-5-7996-3456-7.

12⠄Казанский, Ю. А. Физика ядерных реакторов : учебник для вузов / Ю. А. Казанский, Е. С. Матусевич. — Москва : Издательство МГТУ им. Н. Э. Баумана, 2023. — 512 с. — ISBN 978-5-7038-5890-7.

13⠄Колесников, С. В. Методы Монте-Карло в медицинской физике : монография / С. В. Колесников, Д. В. Акимов. — Санкт-Петербург : Издательство СПбГЭТУ "ЛЭТИ", 2022. — 268 с. — ISBN 978-5-7629-2891-3.

14⠄Красильников, Н. Н. Спектрометрия гамма-излучения : учебное пособие / Н. Н. Красильников, А. П. Ермолаев. — Москва : Издательство МИФИ, 2021. — 198 с. — ISBN 978-5-7262-2765-4.

15⠄Курдин, А. Р. Влияние глобальных вызовов на российский рынок нефти и нефтепродуктов / А. Р. Курдин, М. Ю. Байков, Г. С. Чеблаков, И. В. Пахомова // Научный журнал НИУ ИТМО. Серия: Экономика и экологический менеджмент. — 2020. — №$. — С. $$$-$$$.

$$⠄$$$$$$$$$, $. $. $$$$$$$$$$$$$ $$$$$$$$$ $$$$$$$$ $$$$$$$$$ $ $$$$$$$$$ $$$$ : $$$$$$-$$$$$$$$$$$$ $$$$$$$ / $. $. $$$$$$$$$, $. $. $$$$$$$. — $$$$$$ : $$$$$$$$$$$$ $$$$ $$. $. $. $$$$$$$, $$$$. — $$$ $. — $$$$ $$$-$-$$$$-$$$$-$.

$$⠄$$$$$$ $$$$$-$$$$$ $ $$$$$$$ $$$$$$ : $$$$$$$ $$$$$$$ $$$$$$ / $$$ $$$$$$$$$ $. $. $$$$$$$$$. — $$$$$$ : $$$$$$$$$$$$$$$, $$$$. — $$$ $. — $$$$ $$$-$-$$$-$$$$$-$.

$$⠄$$$$$$, $. $. $$$$$$ $$$$$$$$$$ $ $$$$$$ $$ $$$$$$$$$$$$ $$$$$$$$$ : $$$$$$$ $$$$$$$ / $. $. $$$$$$, $. $. $$$$$$$$. — $$$$$ : $$$$$$$$$$$$ $$$$$$$$ $$$$$$$$$$$$$$$$ $$$$$$$$$$$$, $$$$. — $$$ $. — $$$$ $$$-$-$$$$-$$$$-$.

$$⠄$$$$$$$$$$$ $$$$$$$$$$$ $$$ $$$$$$$$$$$$$ $ $$$$$$$ $$$$$$$$ : $$$$$$$ $$$$$$$ / $. $. $$$$$$$, $. $. $$$$$$$$$$, $. $. $$$$$$, $. $. $$$$$$. — $$$$$-$$$$$$$$$ : $$$$$$$$$$$$ $$$$$, $$$$. — $$$ $. — $$$$ $$$-$-$$$-$$$$$-$.

$$⠄$$$$$$$$$$$$ $$$$$$$$$$$$ $ $$$$$$$$$$ : $$$$$$$ $$$ $$$$$ / $. $. $$$$$$, $. $. $$$$$$$$$, $. $. $$$$$$$$$, $. $. $$$$$$$. — $$$$$$ : $$$$$$$$$$$$ $$$$$, $$$$. — $$$ $. — ($$$$$$ $$$$$$$$$$$). — $$$$ $$$-$-$$$-$$$$$-$.

$$⠄$$$$$$$, $. $. $$$$$$ $ $$$$$$$ $$$$$$ $$ $$$$$$$$$ : $$$$$$$ $$$$$$$ / $. $. $$$$$$$, $. $. $$$$$$$$. — $$$$$$ : $$$$$$$$$$$$ $$$$, $$$$. — $$$ $. — $$$$ $$$-$-$$$$-$$$$-$.

$$⠄$$$$$$$$$$$ $$$$$$ $$$$$$$$$$$$$ $ $$$$$$$$$$$$ $$$$$$ : $$$$$$$$$$ / $$$ $$$$$$$$$ $. $. $$$$$$$$$. — $$$$$$ : $$$$$, $$$$. — $$$ $. — $$$$ $$$-$-$$-$$$$$$-$.

$$⠄$$$$$$$, $. $. $$$$$$$$$$ $$$$$$$$$ $$$$ $$$ $$$$$$$ $$$$$$$ $$$$$$$$$$$$$ : $$$$$$$$$ $$$$$$$$$$$$$$$$ $$$$$$ / $$$$$$$ $$$$$ $$$$$$$$$. — $$$$$$, $$$$. — $$ $.

$$⠄$$$$$$$, $. $. $$$$$$$$$ $$$$$$$$$$$$$ $ $$$$$$$ $$$$$$$$$$$$ $$$$$$ : $$$$$$$ $$$$$$$ / $. $. $$$$$$$, $. $. $$$$$$. — $$$$$$$$$$$ : $$$$$$$$$$$$ $$$$, $$$$. — $$$ $. — $$$$ $$$-$-$$$$-$$$$-$.

$$⠄$$$$$$ $ $$$$$$$ $$$$$$$ $$$$$$$$$ : $$$$$$$ $$$ $$$$$ / $. $. $$$$$$$$$, $. $. $$$$$$$$$, $. $. $$$$$$$, $. $. $$$$$$$$. — $$$$$$ : $$$$$$$$$$$$ $$$$ $$. $. $. $$$$$$$, $$$$. — $$$ $. — $$$$ $$$-$-$$$$-$$$$-$.

$$⠄$$$$$$, $. $. $$$$$$$$$$$$$ $$$$$$$$$$$$$$ $$$$$-$$$$$$$$$ $ $$$$$$$$$ $$$$$$$ $$$$$-$$$$$ : $$$$$$$$$$$ ... $$$$$$$$$ $$$$$$-$$$$$$$$$$$$$$ $$$$ : $$.$$.$$ / $$$$$$ $$$$$$$ $$$$$$$$$$$$. — $$$$$-$$$$$$$$$, $$$$. — $$$ $.

$$⠄$$$$$$$$, $. $. $$$$$$$$$$$$ $$$$$$ : $$$$$$$ $$$ $$$$$ / $. $. $$$$$$$$, $. $. $$$$$$$$. — $$$$$$ : $$$$$$$$$$$$ $$$$$, $$$$. — $$$ $. — ($$$$$$ $$$$$$$$$$$). — $$$$ $$$-$-$$$-$$$$$-$.

$$⠄$$$$$$$ $$$$$$ $ $$$$$$ $$$$$$$$$$$ $$$$$$$$$ : $$$$$$$ $$$$$$$ / $. $. $$$$$$$, $. $. $$$$$$$$$, $. $. $$$$$$$$$, $. $. $$$$$$. — $$$$$$ : $$$$$$$$$$$$ $$$$ $$. $. $. $$$$$$$, $$$$. — $$$ $. — $$$$ $$$-$-$$$$-$$$$-$.

$$⠄$$$, $. $$$$: $ $$$$$$$$$$ $$$$$$$ $$$ $$$ $$$ $$$$$ / $. $$$, $. $$$$$$, $. $$$$$ $$ $$. // $$$$$$$ $$ $$$$$$$$ $$$ $$$$$$$. — $$$$. — $$$. $$, № $$. — $. $$$$-$$$$.

$$⠄$$$$$$, $. $ $$$$$$ $$ $$$ $$$ $$ $$$$ $$$ $$$$$$$ $$$$$$$ $$$ $$$$$$$$$ / $. $$$$$$, $. $$$$$$$, $. $$$$$$$$ $$ $$. // $$$$$$$ $$$$$$$. — $$$$. — $$$. $$, № $. — $. $$$$$$.