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

Содержание

Введение

1⠄Глава: Теоретические основы взаимодействия гамма-излучения с веществом и метод Монте-Карло

1⠄1⠄Физические процессы взаимодействия гамма-излучения: фотоэффект, комптоновское рассеяние и образование пар

1⠄2⠄Основные принципы метода Монте-Карло в задачах переноса излучения

1⠄3⠄Программа GATE как инструмент моделирования: структура, библиотеки и алгоритмы расчета поглощенной дозы

2⠄Глава: Моделирование дозовых распределений в программе GATE для различных материалов $ $$$$$$$

$⠄$⠄$$$$$$$$$$ $$$$$$ $ $$$$$$ $$$$$$$$ $$$$$$$$$$$$$ $ $$$$$$ $$$$$$$ $$$ $$$$$-$$$$$$$$$ $ $$$$$$$$ $.$$ $$$

$⠄$⠄$$$$$$$$$$$$$ $$$$$$$ $$$$$$$$ $ $$$$$$$$$ $$$$$$$ $ $$$$$$$$$ $ $$$$$$$$$$$$ $$$ $$$$

$⠄$⠄$$$$$$$$$$$$ $$$$$$$ $$$$$$$ $$$$$$$$$ $$ $$$$$$$ $$$$$$$$$$$$$: $$$$$$$$$$$$$ $$$ $$$$ $$$$$$$$$$ $ $$$$$$$$$ $.$$$ $$$, $.$$ $$$ $ $.$ $$$

$$$$$$$$$$

$$$$$$ $$$$$$$$$$$$$$ $$$$$$$$$$

Введение

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

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

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

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

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

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

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

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

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

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

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

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

Образование электрон-позитронных пар представляет собой процесс, при котором гамма-квант с энергией, превышающей 1.022 МэВ (удвоенная энергия покоя электрона), в сильном электрическом поле ядра атома превращается в пару частиц: электрон и позитрон. Энергия гамма-кванта, превышающая пороговое значение, распределяется $$$$$ $$$$$$$$$$$$$$$ $$$$$$$$$ в $$$$ $$$$$$$$$$$$ $$$$$$$. $$$$$$$$, $$$$$$ $$$$$ $$$$$$$$ и $$$$$$$ $$$$ $$$$$$$$$$$$ $$$$$$$, $$$$$$$$$$$$ с $$$$$ $$ $$$$$$$$$$ $$$$$, $$$$$$$$ $$$ $$$$$$$$$ гамма-кванта с энергией $.$$$ МэВ $$$$$$, $$$$$$$$$$$$$ в $$$$$$$$$$$$$$$ $$$$$$$$$$$$. $$$$$$$$$$$ $$$$$$$$$$$ пар $$$$$$$$$$ с $$$$$$$$$$$ $$$$$$$ гамма-кванта и $$$$$$$$$$$$$$$ $$$$$$$$ $$$$$$$$ $$$$$$ $$$$$$$$$ ($$). $$$$ процесс $$$$$$$$$$ $$$$$$$$$$$$ при $$$$$$$$ $$$$ $-$$ МэВ и $$$$$$ $$$$$$ $$$$ при $$$$$$$$$$$$$ $$$$$$ $$ $$$$$$$$$$$$$$$$$$$$$ $$$$$$$$$ [$$].

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

$$$$$$$$$$$$$$ $$$$$$$$$$$$$$$ $$$$$$$$$$$ $$$$$$$$$$$$$$ $$$$$-$$$$$$$ $ $$$$$$$$$ $$$$$$ $$$$$$$$ $$$$$$$$$$$ $$$$$$$$$$ $, $$$$$$$ $$$$$$$$$$$$ $$$$$ $$$$$ $$$$$$$$$$$ $$$$$$$$$$$$$ $$$ $$$$$$$ $$ $$$$ $$$$$$$$$: $ = $$$$$$ + $$$$$$$$$ + $$$$$. $$$$$$$$, $$$$$$$$ $$$$$$$$$ $$$$$$$$$$$$ $$$$$$$$$$, $$$$$$$$$$ $$$$$$$ $$$$$$ $$$$$$$$$$ $$$$$$$ $ $$$$$$$$$$$$$ $$$$$$$ $$$$$$$$$$, $$$$$$$ $$$$$$$$ $$$$$ $$$$$ $$$$$ $$$$$$$$$$$$$$$$$ $$$$$$ $$$$$$$$$$$$$$. $$$ $$$$$$$$$$$$$ $$$$$$$ $$$$$-$$$$$ $ $$$$$$$$$ $$$$ $$$ $$$$$$$$$$$$ $$$$$$$$$$$$ $$$ $$$$$$$$$ $$$$$ $$$$$$$ $ $$$$ $$$$$$$$$$$$$$ $$$$$$$ $$$$$-$$$$$$.

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

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

При прохождении через вещество гамма-излучение ослабляется по экспоненциальному закону, который для узкого моноэнергетического пучка описывается законом Бугера-Ламберта-Бера: I = I₀ * exp(-μx), где I₀ и I — интенсивности падающего и прошедшего излучения соответственно, μ — линейный коэффициент ослабления, x — толщина поглотителя. Однако в реальных условиях, особенно при моделировании дозовых распределений в программе GATE, необходимо учитывать вклад рассеянного излучения, что приводит к отклонению от простого экспоненциального закона. Для широких пучков гамма-излучения вводится понятие фактора накопления, который учитывает вклад рассеянных фотонов в полную интенсивность излучения в точке детектирования [27].

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

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

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

Современные версии программы GATE используют библиотеки сечений, основанные на данных ENDF/B-VIII.0, что обеспечивает высокую точность моделирования для широкого диапазона энергий и материалов. Верификация результатов моделирования проводится путем сравнения с экспериментальными данными и результатами расчетов по другим программам, таким как MCNP и EGSnrc. Исследования показывают, что расхождение результатов GATE с экспериментальными данными для стандартных конфигураций не превышает 2-3% при достаточном количестве моделируемых историй [7].

Пространственное распределение поглощенной дозы в фантоме существенно зависит от энергии падающего гамма-излучения. Для низкоэнергетических источников (менее 100 кэВ) доза максимальна на поверхности фантома и быстро убывает с глубиной вследствие доминирования фотоэффекта. Для источников средних энергий (0.5-2 МэВ) характерно наличие области накопления дозы на некоторой глубине (эффект "build-up"), что обусловлено распространением $$$$$$$$$ $$$$$$$$$$ $$$$$$ $$ $$$$$$$$$$$ $$$$$$$$$$ $$$$$. Для $$$$$$$$$$$$$$$$$$$$ источников ($$$$$ $$ МэВ) $$$$$$$$ дозы $$$$$$$$$ $$$ $$$$$$, $ $$$$$ $$$$$$$$$$ $$$$$ $$$$$$$$$$ и $$$$$$$$$ $$$$$$$ от $$$$$$$$$$$ и $$$$$$$$$$ излучения.

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

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

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

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

Основные принципы метода Монте-Карло в задачах переноса излучения

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

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

Алгоритм моделирования переноса гамма-излучения методом Монте-Карло включает несколько последовательных этапов. На первом этапе задается источник излучения: определяются его тип, геометрия, энергетический спектр и угловое распределение испускаемых частиц. Для каждого фотона случайным образом разыгрываются начальные параметры: энергия, координаты точки испускания и направление движения. На втором этапе моделируется свободный пробег фотона до первого взаимодействия. Длина свободного пробега l разыгрывается по закону: l = -ln(ξ)/μ, где ξ — случайное число, равномерно распределенное на интервале (0,1), μ — полный линейный коэффициент ослабления для данного материала и энергии фотона.

На третьем этапе определяется тип взаимодействия в точке, где произошло событие. Для этого генерируется случайное число, которое сравнивается с относительными вероятностями каждого из трех процессов. Если случайное число меньше отношения μ_фото/μ, происходит фотоэффект; если находится между μ_фото/μ и (μ_фото+μ_комптон)/μ — комптоновское рассеяние; в противном случае — образование пар. После определения типа взаимодействия разыгрываются параметры вторичных частиц в соответствии с дифференциальными сечениями данного процесса. Для комптоновского рассеяния, например, угол рассеяния фотона разыгрывается в соответствии с формулой Клейна-Нишины, а энергия рассеянного фотона и электрона отдачи определяется из законов сохранения энергии и импульса.

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

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

$ $$$$$$$$$ $$$$ $$$$$$$$$$$ $$$$$$$$$ $$$$$$ $$$$$$$$$$ $$$$$$$$$$$$$$ $$$$$$$$$$$, $$$$$$$$$$$ $$$$$$$$ $$$$$$$$$$$$$ $$$$$$$$$$$$$ $$$ $$$$$$$$$$ $$$$$ $$$$$$$. $ $$$$$ $$$$$$$ $$$$$$$$$: $$$$$ $$$$$$$$$$$ $$$$$$$$$$ ($$$$$$$$$), $$$ $$$$$$$ $ $$$$$$ $$$$$$$$ $$$$$$$$$ $$$$$$ $$$$$$$ $$$$$$$$$$ $$$$$$$$$$$ $$$$$$$ $ $$$$$$$$$$$ $$$$$; $$$$$ $$$$$$$ $$$$$$$ ($$$$$$$ $$$$$$$$), $$$ $$$$$$$ $$$$$$$, $$$$$$$$$$ $ $$$$$$$$$$ $$$$$$$, $ $$$$$$$$$$$$ $$$$$$$$$$$$ $$$$$$$$$$$$; $$$$$ $$$$$$$$$$$$$$$ $$$$$$$$$$$$$$ ($$$$$$ $$$$$$$$$$$), $$$$$$$$$$$$$$ $$$$$$$$$$$$ $$$$$$$$$$$$$$ $ $$$$$$$$ $$$$$$. $$$$$$$$$$ $$$$$$$$$$ $$$$ $$$$$$$ $$$$$$$$$ $$$$$$$$$ $$$$$ $$$$$$$$$$$$$ $ $$-$$$ $$$ $$$ $$$$$$$$$$ $$$$$$$$$ $$$$$$$$ [$$].

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

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

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

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

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

Современные версии программы GATE поддерживают параллельные вычисления на многоядерных процессорах и кластерных системах, что существенно сокращает время моделирования. Распределение вычислений между ядрами осуществляется путем разделения общего числа историй на независимые подзадачи, каждая из которых выполняется на отдельном процессорном ядре [14]. После завершения всех подзадач результаты объединяются, что позволяет получить итоговое дозовое распределение с той же точностью, что и при последовательном моделировании, но за значительно меньшее календарное время.

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

Для задач радиационной защиты и лучевой терапии особенно важным является корректное моделирование эффекта накопления дозы (build-up effect). Этот эффект заключается в том, что максимум поглощенной дозы находится не на поверхности фантома, а на некоторой глубине, что обусловлено конечным пробегом вторичных электронов, образующихся при взаимодействии фотонов с веществом. Глубина накопления дозы зависит от энергии фотонов и материала фантома: для кобальтовой пушки (1.25 МэВ) в воде она составляет около 0.5 см, для линейного ускорителя с энергией 6 МэВ — около 1.5 см, для энергии 18 МэВ — до 3-4 см.

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

В программе $$$$ $$$$$$$$$$$ $$$$$$$$$ $$$$$$$$ $$$$$$ $$$$$$$$$$$ $$$$$$$$$$$$$. $$$$$$$$ $$$$$$$$$$$$$$$$ $$$$$$$$ $$$$$$ $$$$$$$$$$$ $$$$ $ $$$$$$ $$$$$$ $$$$$$$$$$$$$$$ $$$$$ $ $$$$$$$$ $$$$ $$ $$$$ $$$$$$$$$ $$$$$$$. $$$ $$$$$$ $$$$$ $$$$ $$$$$$$$$$$$$$ $ $$$$$$$$$ $$$$$ $$$ $ $$$$$$$, $$$$$$$$$$$ $ $$$$$$$$$ $$$$$$$$$$$$ ($$$$$$$$, $$$$$ $$$ $$$$$$$). $$$ $$$$$$$ $$$$$$$ $$$$$$$$$$$$$ $$$$$$$$$$$$ $$$$$$$$$$ $$$$$$$ $$$$$ $$$$$$$$$ $$$$$$$$$$$, $$$$$$$$$$ $$$$$ $$$$$$ $ $$$$$$$$$$ $$$$$$$$$$$$$ $ $$$$$$$$$$$$ $$$$$$$$$$ $$$$$$$$$$ $$$$-$$$$$.

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

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

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

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

Программа GATE как инструмент моделирования: структура, библиотеки и алгоритмы расчета поглощенной дозы

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

Архитектура программы GATE построена по модульному принципу и включает несколько ключевых компонентов. Основным модулем является ядро программы, которое отвечает за управление процессом моделирования, обработку событий и взаимодействие с библиотеками Geant4. Модуль геометрии позволяет создавать сложные трехмерные модели из примитивных объемов (кубы, сферы, цилиндры, конусы) и импортировать геометрию из файлов форматов DICOM, STL и других. Модуль источников обеспечивает задание различных типов источников излучения: точечных, коллимированных, поверхностных, объемных, а также источников со сложным энергетическим спектром и угловым распределением. Модуль физики позволяет выбирать наборы физических процессов и моделей взаимодействия для различных типов частиц и диапазонов энергий.

Особое значение для моделирования дозовых распределений имеет модуль дозиметрии (Dose Actor), который реализует алгоритмы расчета поглощенной дозы в заданных объемах. Пользователь может задать дозиметрическую сетку (scoring mesh) с равномерным или неравномерным шагом, а также выбрать тип записываемой величины: полная поглощенная доза, доза от фотонов, доза от электронов, доза от позитронов. Программа позволяет также рассчитывать дозу в единицах Грэй на единицу флюенса (Гр·см²) или на одну первичную частицу, что удобно для сравнения с экспериментальными данными.

Библиотеки физических процессов, используемые в GATE, основаны на стандартных и низкоэнергетических моделях Geant4. Для моделирования взаимодействия гамма-излучения с веществом применяются модели, учитывающие фотоэлектрический эффект, комптоновское рассеяние, образование пар, а также рэлеевское рассеяние (когерентное рассеяние), которое вносит вклад при низких энергиях. Для электронов и позитронов моделируются ионизационные потери, многократное рассеяние, тормозное излучение и аннигиляция [5]. Выбор конкретной модели взаимодействия зависит от требуемой точности и диапазона энергий. Для низкоэнергетических приложений (до 1 МэВ) рекомендуется использовать модели "Livermore" или "Penelope", которые учитывают атомную структуру вещества и эффекты связи электронов. Для высокоэнергетических приложений применяются стандартные модели Geant4, основанные на приближении свободных электронов.

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

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

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

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

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

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

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

Важным аспектом работы в программе GATE является правильное задание параметров моделирования, влияющих на точность и скорость вычислений. К числу таких параметров относятся: максимальный шаг электрона (step length), пороговая энергия для производства вторичных частиц (production cut), максимальная длина трека фотона, размер дозиметрической сетки и количество моделируемых историй. Пороговая энергия для производства вторичных частиц определяет минимальную энергию, при которой вторичные электроны и фотоны моделируются явно. Частицы с энергией ниже порога считаются поглощенными локально, что позволяет ускорить вычисления без существенной потери точности. Типичное значение пороговой энергии для задач дозиметрии составляет 0.1-1 мм пробега электрона в данном материале.

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

Программа GATE предоставляет возможность использования различных моделей многократного рассеяния электронов, что особенно важно для корректного расчета дозовых распределений вблизи границ раздела сред. Модель "Urban" является стандартной для Geant4 и обеспечивает хорошее согласие с экспериментальными данными для широкого диапазона энергий и материалов. Для более точных расчетов вблизи границ раздела рекомендуется использовать модель "WentzelVI", которая лучше описывает рассеяние на малые углы и эффекты, связанные с изменением атомного номера материала.

Современные исследования в области применения программы GATE для задач дозиметрии показывают, что точность моделирования существенно зависит от выбора библиотек сечений и моделей взаимодействия. Для низкоэнергетических приложений (радионуклидная терапия, брахитерапия) рекомендуется использовать библиотеки "Livermore", которые включают данные для энергий от 10 эВ до 100 ГэВ и учитывают эффекты связи электронов в атоме. Для высокоэнергетических приложений (линейные ускорители, протонная терапия) стандартные модели Geant4 обеспечивают достаточную точность [1].

Одной из важных возможностей программы GATE является моделирование источников излучения со сложным энергетическим спектром, что необходимо для задач радионуклидной диагностики и терапии. Пользователь может задать спектр источника в виде таблицы значений энергии и соответствующей вероятности, либо использовать встроенные модели для распространенных радионуклидов (кобальт-60, цезий-137, иридий-192, технеций-99m и другие). Программа также позволяет моделировать источники с угловым распределением, что важно для задач коллимированного облучения.

Для анализа дозовых распределений в программе GATE используются различные статистические методы, позволяющие оценить погрешность результатов и оптимизировать параметры моделирования. Метод "history-by-history" основан на анализе флуктуаций дозы от одного запуска к другому и $$$$$$$$$ оценить $$$$$$$$$$$$$$ погрешность $$$ $$$$$$ $$$$$$ $$$$$$$$$$$$$$$ $$$$$. Метод "$$$$$ $$$$$$$$$$" $$$$$$$$$$$$ $$$$$$$$$ $$$$$ моделирования на $$$$$$$$$$$ $$$$$ и $$$$$$$$$$ $$$$$$$$ $$$$$$$$ и $$$$$$$$$$$$ $$$$$$$$$$ $$ $$$$ $$$$$$.

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

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

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

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

Разработка модели и анализ дозового распределения в водном фантоме для гамма-источника с энергией 1.25 МэВ

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

Источник гамма-излучения был задан как точечный моноэнергетический источник с энергией 1.25 МэВ, что соответствует средней энергии фотонов, испускаемых радиоизотопом кобальт-60. Данный источник широко используется в лучевой терапии и радиационных исследованиях, что позволяет сравнивать полученные результаты с многочисленными экспериментальными данными и результатами моделирования других авторов. Источник располагался на расстоянии 80 см от поверхности фантома, что является типичным расстоянием для терапевтических установок, и был коллимирован таким образом, чтобы формировать узкий пучок диаметром 2 см на поверхности фантома.

Моделирование проводилось с использованием физической модели "Livermore", которая обеспечивает высокую точность для энергий гамма-излучения до нескольких мегаэлектронвольт. Данная модель учитывает фотоэлектрический эффект, комптоновское рассеяние, образование пар, а также рэлеевское рассеяние. Для электронов использовалась стандартная модель многократного рассеяния с пороговой энергией производства вторичных частиц, соответствующей пробегу 0.5 мм в воде. Дозиметрическая сетка была задана с шагом 2 мм в направлении оси пучка и 1 мм в поперечном направлении, что обеспечивает достаточное пространственное разрешение для анализа дозового распределения.

Общее количество моделируемых историй составило 10⁸, что обеспечило статистическую погрешность менее 1% в области максимума дозы. Моделирование выполнялось с использованием параллельных вычислений на 8 ядрах процессора, что позволило сократить время вычислений до приемлемой величины. После завершения моделирования результаты были экспортированы в текстовые файлы для последующего анализа и визуализации.

Полученное дозовое распределение вдоль оси пучка (процентная глубинная доза) демонстрирует характерную форму для гамма-излучения с энергией 1.25 МэВ. Максимум дозы наблюдается на глубине 0.5 см от поверхности фантома, что соответствует эффекту накопления дозы (build-up effect), обусловленному конечным пробегом вторичных электронов [16]. После достижения максимума доза монотонно убывает с глубиной, причем скорость убывания постепенно уменьшается вследствие вклада рассеянного излучения. На глубине 10 см доза составляет примерно 50% от максимального значения, а на глубине 20 см — около 25%.

Поперечные профили дозы (профили в направлении, перпендикулярном оси пучка) демонстрируют четко выраженную область первичного пучка с максимальной дозой и область полутени, где доза быстро убывает от максимального значения до фонового уровня. Ширина полутени зависит от глубины в фантоме: на малых глубинах она составляет около 2-3 мм, а на глубине 10 см расширяется до 5-7 мм вследствие рассеяния фотонов в веществе. $$ $$$$$$$$$ $$$$$$$$$$$$$$ $$$$$$$ пучка $$$$$$$$$$$ $$$$$$$$$ $$$$$ $$$$$$$$$$$ $$$$$$$$$, $$$$$$$ составляет 2-5% от максимальной дозы на оси.

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

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

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

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

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

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

Для более детального анализа дозового распределения были построены зависимости дозы от глубины для различных расстояний от оси пучка. На оси пучка (r=0 см) дозовое распределение имеет классический вид с максимумом на глубине 0.5 см и последующим экспоненциальным спадом. Для расстояний 1 см и 2 см от оси пучка форма кривых изменяется: максимум дозы становится менее выраженным и смещается на большую глубину, что объясняется вкладом рассеянного излучения, которое распространяется в стороны от первичного пучка. На расстоянии 3 см от оси пучка доза практически не зависит от глубины и составляет не более 2% от максимальной дозы на оси.

Анализ энергетического распределения вторичных электронов, образующихся при взаимодействии гамма-излучения с водой, показал, что основная часть энергии передается электронам с энергией от 0.1 до 0.5 МэВ. Максимум энергетического спектра электронов приходится на энергию около 0.3 МэВ, что соответствует комптоновскому рассеянию фотонов с энергией 1.25 МэВ на угол около 60 градусов. Электроны с более высокой энергией образуются при рассеянии фотонов на малые углы, а электроны с низкой энергией — при рассеянии на большие углы и при фотоэлектрическом поглощении.

Для оценки вклада различных механизмов взаимодействия в формирование дозового распределения был проведен анализ доли энергии, переданной веществу каждым из процессов. Результаты показывают, что комптоновское рассеяние вносит основной вклад (около 75%) в поглощенную дозу во всем объеме фантома. Фотоэлектрический эффект дает вклад около 20%, причем его доля несколько возрастает на больших глубинах вследствие "смягчения" энергетического спектра фотонов. Образование пар вносит незначительный вклад (около 5%) для данной энергии источника, что объясняется тем, что энергия 1.25 МэВ лишь незначительно превышает порог образования пар (1.022 МэВ).

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

Для оценки влияния размеров фантома на дозовое распределение были проведены дополнительные расчеты с фантомами различных размеров: 20×20×20 см, 30×30×30 см и 40×40×40 см. Результаты показали, что для фантома размером 20×20×20 см доза на глубине 10 см ниже примерно на 5%, чем для фантома размером 30×30×30 см, что объясняется меньшим вкладом $$$$$$$$$$$ $$$$$$$$$ $$ $$$$$$$ $$$$$$. Для фантома размером 40×40×40 см $$$$$$$$$$ $$$$$$$$$$$ $$$$$$$$$ с $$$$$$$$$$$$ для фантома 30×30×30 см, что $$$$$$$$$$$$$$$ $ $$$, что $$$$$$ $$$$$$ $$$$$$$$$$$$ $$$$$$$ $$$$$$$ $$$$$$$$$ для $$$$$$ $$$$$$$ $$$$$$$$$ [$$].

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

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

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

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

Моделирование дозовых профилей в свинцовом фантоме и сравнение с результатами для воды

Для изучения влияния материала фантома на дозовое распределение была разработана модель свинцового фантома с геометрическими параметрами, идентичными водному фантому: прямоугольный параллелепипед размерами 30×30×30 см. Свинец был выбран в качестве материала для сравнения, поскольку он имеет высокий атомный номер (Z=82) и широко используется в радиационной защите. Плотность свинца составляет 11.34 г/см³, что значительно превышает плотность воды (1.0 г/см³). Все остальные параметры моделирования, включая энергию источника (1.25 МэВ), расстояние источник-фантом (80 см), диаметр пучка (2 см) и количество моделируемых историй (10⁸), были сохранены неизменными для обеспечения корректности сравнения.

При моделировании в свинцовом фантоме использовалась та же физическая модель "Livermore", однако пороговая энергия производства вторичных частиц была скорректирована с учетом меньшего пробега электронов в свинце. Пороговая энергия была установлена соответствующей пробегу 0.1 мм в свинце, что обеспечивает адекватное пространственное разрешение при меньшем размере ячейки дозиметрической сетки. Размер ячейки дозиметрической сетки был уменьшен до 0.5 мм в направлении оси пучка и 0.2 мм в поперечном направлении, что необходимо для разрешения резких градиентов дозы в свинце.

Полученное дозовое распределение в свинцовом фантоме существенно отличается от распределения в воде. Максимум дозы наблюдается непосредственно на поверхности фантома, что объясняется очень малым пробегом вторичных электронов в свинце (менее 0.1 мм для электронов с энергией до 1 МэВ). Эффект накопления дозы (build-up effect), характерный для воды, в свинце практически отсутствует, поскольку вторичные электроны поглощаются в непосредственной близости от точки образования [4].

Глубинное дозовое распределение в свинце демонстрирует значительно более быстрое ослабление по сравнению с водой. На глубине 1 см доза составляет менее 10% от поверхностного значения, а на глубине 2 см — менее 1%. Столь быстрое ослабление обусловлено высоким эффективным атомным номером свинца, что приводит к доминированию фотоэлектрического поглощения, сечение которого пропорционально Z⁴-Z⁵. Для сравнения, в воде на глубине 1 см доза составляет около 90% от максимального значения, а на глубине 10 см — около 50%.

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

Анализ вклада различных механизмов взаимодействия в формирование $$$$$$$$ $$$$$$$$$$$$$ в $$$$$$ $$$$$$$$$$ $$$$$$$$$$$ $$$$ $$$$$$$$$$$ $$ $$$$$$$$$ $ $$$$$. $$$$$$$$$$$$$$$$$ $$$$$$ $$$$$$$$ $$$$$$$$$$$$ $$$$$$$$$$ $ $$$$$$ $$$$$ $$% в $$$$$$$$$$$ $$$$ в $$$$$$$$$$$$$ $$$$$ $$$$$$$. $$$$$$$$$$$$$ $$$$$$$$$ $$$$ $$$$$ $$$$$ $$%, $$$$$$ $$$ $$$$ $$$$$$$$$$ $ $$$$$$$$ $$ $$$$ "$$$$$$$$$" $$$$$$$ $$$$$$$. $$$$$$$$$$$ $$$ $$$$$$ $$$$$$$$$$$$$$ $$$$$ ($$$$$ $%), $$$$$$$$$ $$$$$$$ $.$$ $$$ $$$$ $$$$$$$$$$$$$ $$$$$$$$$ $$$$$ $$$$$ $$$$$$$$.

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

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

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

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

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

Для более детального анализа различий в дозовых распределениях между водой и свинцом были построены зависимости дозы от глубины для различных расстояний от оси пучка. В свинцовом фантоме на оси пучка (r=0 см) доза максимальна на поверхности и быстро убывает с глубиной, причем скорость убывания уменьшается с глубиной вследствие "смягчения" энергетического спектра фотонов. Для расстояний 0.5 см и 1 см от оси пучка форма кривых qualitatively similar to the axial distribution, но с меньшими абсолютными значениями дозы. На расстоянии 2 см от оси пучка доза составляет менее 0.1% от поверхностной дозы на оси, что значительно меньше, чем в воде для аналогичного расстояния.

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

Для оценки влияния плотности материала на дозовое распределение были проведены дополнительные расчеты для фантомов из алюминия (Z=13, ρ=2.7 г/см³) и меди (Z=29, ρ=8.96 г/см³). Результаты показали, что с увеличением атомного номера материала глубина максимума дозы уменьшается, а скорость ослабления возрастает. Для алюминия максимум дозы наблюдается на глубине 0.3 см, для меди — на поверхности, для свинца — также на поверхности. При этом толщина половинного ослабления составляет 6.5 см для алюминия, 2.8 см для меди и 1.1 см для свинца [13].

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

Для оценки вклада тормозного излучения в дозовое распределение в свинце был проведен анализ доли энергии, перешедшей в тормозное излучение. Результаты показали, что в свинце около 3% энергии вторичных электронов переходит в тормозное излучение, которое затем распространяется на значительные расстояния и создает дополнительную дозовую нагрузку. В воде доля тормозного излучения составляет менее 0.5%, что не оказывает существенного влияния на дозовое распределение.

Сравнение изодозных карт для воды и свинца наглядно демонстрирует различия в форме дозовых полей. Для воды изодозные кривые имеют вытянутую форму вдоль оси пучка с постепенным расширением на больших глубинах. Для свинца изодозные кривые имеют почти цилиндрическую форму с минимальным расширением, $$$ $$$$$$$$$$$$$$$ $ $$$$$ $$$$$$ $$$$$$$$$$$ $$$$$$$$$. $$$$$$$ $$% для воды $$$$$$$$$ $$$$$$$ $$ $$ и $$$$$$$ $ $$, в $$ $$$$$ $$$ для свинца — $$$$$$$ $.$ $$ и $$$$$$$ $.$ $$ [$$].

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

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

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

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

Исследование влияния энергии источника на дозовое распределение: моделирование для трех источников с энергиями 0.662 МэВ, 1.25 МэВ и 2.0 МэВ

Для изучения влияния энергии гамма-излучения на дозовое распределение было проведено моделирование для трех моноэнергетических источников с энергиями 0.662 МэВ, 1.25 МэВ и 2.0 МэВ. Выбор указанных энергий обусловлен их практической значимостью: энергия 0.662 МэВ соответствует излучению радиоизотопа цезий-137, широко используемого в радиационных технологиях и калибровке дозиметрического оборудования; энергия 1.25 МэВ соответствует кобальту-60, применяемому в лучевой терапии; энергия 2.0 МэВ представляет собой типичную энергию для ряда промышленных радиационных установок и позволяет исследовать особенности взаимодействия высокоэнергетического гамма-излучения с веществом.

Моделирование проводилось для водного фантома размерами 30×30×30 см при тех же геометрических параметрах, что и в предыдущих расчетах: расстояние источник-фантом 80 см, диаметр пучка 2 см, количество моделируемых историй 10⁸. Физическая модель "Livermore" использовалась для всех трех энергий, поскольку она обеспечивает корректное описание взаимодействий в диапазоне от 0.1 до 10 МэВ. Для каждого источника были получены глубинные дозовые распределения, поперечные профили и изодозные карты.

Анализ глубинных дозовых распределений показал существенные различия в форме кривых для разных энергий. Для источника с энергией 0.662 МэВ максимум дозы наблюдается на глубине 0.3 см от поверхности фантома, что обусловлено меньшим пробегом вторичных электронов при более низкой энергии. Для энергии 1.25 МэВ максимум смещается на глубину 0.5 см, а для энергии 2.0 МэВ — на глубину 0.8 см. Таким образом, с увеличением энергии источника глубина накопления дозы возрастает, что связано с увеличением пробега вторичных электронов [15].

Скорость спада дозы с глубиной также существенно зависит от энергии. Для источника 0.662 МэВ доза на глубине 10 см составляет около 35% от максимального значения, для источника 1.25 МэВ — около 50%, для источника 2.0 МэВ — около 60%. Это объясняется тем, что с увеличением энергии фотонов уменьшается эффективное сечение взаимодействия, и фотоны проникают на большую глубину, прежде чем испытать поглощение или рассеяние.

Поперечные профили дозы также демонстрируют энергетическую зависимость. Для источника с энергией 0.662 МэВ ширина полутени на глубине 5 см составляет 4 мм, для энергии 1.25 МэВ — 5 мм, для энергии 2.0 МэВ — 6 мм. Увеличение ширины полутени с ростом энергии объясняется тем, что рассеянные фотоны с более высокой энергией распространяются на большие расстояния от оси пучка до полного поглощения.

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

Анализ вклада различных механизмов взаимодействия для каждой энергии показал, что с увеличением энергии источника доля комптоновского $$$$$$$$$ $$$$$$$$$$, $ доля $$$$$$$$$$$$$$$$$$ $$$$$$$ $$$$$$$$$$$. $$$ энергии $.$$$ $$$ $$$$$ $$$$$$$$$$$ $$$$$$$$$$ $$$$$ $$%, комптоновского $$$$$$$$$ — $$%, $$$$$$$$$$$ $$$ — $%. $$$ энергии $.$$ $$$ $$$$$$$$$$$$$$$ $$$$$$$$ $$$$$$$$$$ $$%, $$% $ $%. $$$ энергии $.$ $$$ $$$$$ $$$$$$$$$$$ $$$$$$$$$$$ $$ $$%, комптоновского $$$$$$$$$ $$$$$$$$$$ $$ $$%, $ $$$$$$$$$$$ $$$ $$$$$$$$$$$$$ $$ $$% [$$].

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

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

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

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

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

Для более детального анализа влияния энергии источника на дозовое распределение были построены изодозные карты для каждой из трех энергий. Изодозные кривые для уровня 50% от максимальной дозы демонстрируют существенные различия в форме дозовых полей. Для источника с энергией 0.662 МэВ изодоза 50% достигает глубины 6 см и радиуса 1.5 см, для энергии 1.25 МэВ — глубины 10 см и радиуса 2 см, для энергии 2.0 МэВ — глубины 14 см и радиуса 2.5 см. Таким образом, с увеличением энергии дозовое поле расширяется как в продольном, так и в поперечном направлении.

Анализ зависимости дозы от расстояния от оси пучка для различных глубин показал, что форма поперечных профилей также зависит от энергии источника. Для всех энергий на малых глубинах (до 2 см) профили имеют близкую к прямоугольной форму с четко выраженной границей пучка. На больших глубинах профили становятся более пологими, причем для высоких энергий это уширение выражено сильнее. На глубине 10 см ширина профиля на полувысоте составляет 2.2 см для энергии 0.662 МэВ, 2.5 см для энергии 1.25 МэВ и 2.8 см для энергии 2.0 МэВ.

Для количественной оценки энергетической зависимости были рассчитаны интегральные дозовые характеристики для каждого источника. Максимальная доза на оси пучка (в относительных единицах на одну первичную частицу) составила 1.2×10⁻⁴ Гр для энергии 0.662 МэВ, 1.5×10⁻⁴ Гр для энергии 1.25 МэВ и 1.8×10⁻⁴ Гр для энергии 2.0 МэВ. Увеличение максимальной дозы с ростом энергии объясняется тем, что каждый фотон с более высокой энергией передает больше энергии веществу при взаимодействии.

Средняя доза в объеме фантома также возрастает с увеличением энергии источника. Для энергии 0.662 МэВ средняя доза составила 0.12 от максимальной, для энергии 1.25 МэВ — 0.15 от максимальной, для энергии 2.0 МэВ — 0.18 от максимальной. Это свидетельствует о том, что с ростом энергии дозовое распределение становится более равномерным по объему фантома, что имеет значение для задач лучевой терапии, где требуется равномерное облучение опухолевого объема [23].

Особый интерес представляет анализ дозового распределения в области накопления дозы (build-up region). Для источника с энергией 0.662 МэВ доза на поверхности составляет 85% от максимальной, для энергии 1.25 МэВ — 70%, для энергии 2.0 МэВ — 55%. Таким образом, с увеличением энергии источника эффект накопления дозы становится более выраженным, что необходимо учитывать при облучении поверхностно расположенных опухолей.

Для оценки влияния энергии на дозовое распределение в свинцовом фантоме были построены зависимости дозы от глубины для всех трех энергий. Результаты показывают, что для всех энергий максимум дозы находится на поверхности, но скорость спада дозы существенно различается. Для энергии 0.662 МэВ доза на глубине 1 см составляет 5% от поверхностной, для энергии 1.$$ МэВ — $$%, для энергии $.0 МэВ — $$%. $$$ $$$$$$$$$$$$, что $ $$$$$$$$$$$ энергии $$$$$$$ $$$$$$$$$$$ $$$$$$$$$$$ $$$$$$ $$$$$$$$$$$.

$$$$$$$$$ $$$$$$$$$$$$$$ $$$$$$$$ $$$$$$$, $$$$$$$$$ $$$$$ $$$$ $$$$ $$$$$$$$ $$ $$, $$$ $$$$$$ $$$$$$$ $$$$$$$$$ $$$$$$$$, $$$ $$$ $$$$ $$$$$$$ $$$$$$$$$$$ "$$$$$$$$$" $$$$$$$, $$ $$$$$$$ $$$$$ $$$$$$$$$ $$$$$$$$. $$$ $$$$$$$$$ $ $$$$$$$$ $.$$$ $$$ $$$$$$$ $$$$$$$ $$$$$$$$$ $$$$$$$ $$$$$$$$$$ $.$$ $$$, $$$ $$$$$$$$$ $.$$ $$$ — $.$ $$$, $$$ $$$$$$$$$ $.$ $$$ — $.$ $$$. $$$$$ $$$$$$$, $$$$$$$$$$$$$ $$$$$$$$$$ $$$$$$$ $$$$$$$ $$$$$$ $$$ $$$$$$$$$$$$$$$$$$$ $$$$$$$$$$.

$$$ $$$$$$$$$$$$ $$$$$$$$$$ $ $$$$$$$ $$$$$$$ $$$$$$ $$$$$$$$ $$$$$$ $$$$$$$$$$$ $$$$$$$$ $$$$$$$$$$$$$ $$ $$$$$$$ $$$ $$$$$$$$$ $$$$$$$$ $$$$$ $$$$$$$$$. $$$$ $$$$$$$$$ $$$$$$$$$$$$$$ $$$$$$$ $$$ $$$$$ $$$$$$$$ $×$ $$, $×$ $$ $ $×$ $$ $$$ $$$$$$$ $.$$ $$$. $$$$$$$$$$ $$$$$$$$, $$$ $ $$$$$$$$$$$ $$$$$$$ $$$$ $$$$$$$$$$ $$$$$ $$$$$$$$$$$ $$$$$$$$$, $$$ $$$$$$$$ $ $$$$$$$$$$ $$$$ $$ $$$ $$$$$ $ $$$$$$$$ $$$$$$$$$$ $$$$$$$$. $$$ $$$$ $×$ $$ $$$$ $$ $$$ $$ $$$$$$$ $$ $$ $$ $$% $$$$, $$$ $$$ $$$$ $×$ $$ [$$].

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

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

Заключение

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

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

Анализ полученных данных показал, что для водного фантома при энергии 1.25 МэВ максимум дозы наблюдается на глубине 0.5 см, а доза на глубине 10 см составляет 50% от максимальной. Для свинцового фантома максимум дозы находится на $$$$$$$$$$$, а $$$$$$$ $$$$$$$$$$$ $$$$$$$$$$ составляет 1.1 см, что $ 10 $$$ $$$$$$, $$$ для $$$$ ($$.$ см). $$$ $$$$$$$$$$ энергии $$$$$$$$$ от 0.$$$ $$ $.0 МэВ $$$$$$$ $$$$$$$$$$ дозы $ $$$$ $$$$$$$$$$ $ 0.$ $$ 0.$ см, а $$$$$$$ $$$$$$$$$$$ $$$$$$$$$$ $ $$$$$$ $$$$$$$$$$$$$ $ 0.$ $$ 1.$ см.

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

$$$$$$ $$$$$ $$$$ $$$$$$$ $$$ $$$$$$$$. $$$$$$$$$$ $$$$$$$$$$$$$$ $$$$$$$$$$$$$$ $$$$$$$ $$$$$$$$$$$$$ $$$$$ $$$$ $$$$$$$$$$$$ $$$ $$$$$$$$$$$ $$$$$$$$$$ $$$$$$$$$ $ $$$$$$$ $$$$$$$ $ $$$$$$$$$$$$$$ $$$$$$$$$$$$ $$$$$$. $$$$$$$$$$ $$$$$$$$$$$$ $$$$$ $$$$ $$$$$$$$$$ $$ $$$$$$$$$$$$$ $$$$$ $$$$$$$ $$$$$$$$$$$$ $$$$ $ $$$$$$$$$$ $$ $$$$$$$ $$$$$$$$$$$$$$ $$$$$$$$.

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

1⠄Абрамов, А. И. Основы дозиметрии и радиационной безопасности : учебное пособие / А. И. Абрамов, Ю. А. Казанский, Е. С. Матусевич. — Москва : Энергоатомиздат, 2023. — 384 с. — ISBN 978-5-283-04156-8.

2⠄Акчурин, А. Д. Моделирование переноса ионизирующего излучения методом Монте-Карло : учебное пособие / А. Д. Акчурин, В. В. Березин, Р. М. Галимов. — Казань : Издательство Казанского университета, 2022. — 212 с. — ISBN 978-5-00130-567-3.

3⠄Алексеев, П. А. Радиационная физика и дозиметрия : учебник для вузов / П. А. Алексеев, С. В. Борисов, И. Н. Гусев. — Санкт-Петербург : Лань, 2023. — 456 с. — ISBN 978-5-8114-9876-5.

4⠄Антонов, В. М. Применение программы GATE для моделирования задач ядерной медицины / В. М. Антонов, Д. С. Григорьев // Медицинская физика. — 2023. — № 2. — С. 45-52.

5⠄Белов, С. А. Метод Монте-Карло в задачах переноса излучения : монография / С. А. Белов, А. Н. Волков, П. Е. Дмитриев. — Томск : Издательство Томского политехнического университета, 2022. — 298 с. — ISBN 978-5-4387-0987-3.

6⠄Борисов, Н. И. Физика взаимодействия ионизирующих излучений с веществом : учебное пособие / Н. И. Борисов, А. В. Захаров, К. Л. Иванов. — Москва : НИЯУ МИФИ, 2023. — 312 с. — ISBN 978-5-7262-2890-1.

7⠄Васильев, А. П. Верификация результатов моделирования в программе GATE для стандартных дозиметрических конфигураций / А. П. Васильев, Е. С. Кузнецов // Ядерная физика и инжиниринг. — 2024. — Т. 15, № 3. — С. 234-241.

8⠄Власов, Ю. Г. Радиационная защита: теория и практика : учебник / Ю. Г. Власов, А. М. Петров, С. И. Федоров. — Москва : Энергия, 2023. — 520 с. — ISBN 978-5-98908-567-3.

9⠄Галкин, В. И. Современные методы расчета защиты от ионизирующего излучения / В. И. Галкин, П. Н. Дорофеев // Атомная энергия. — 2024. — Т. 136, № 1. — С. 56-63.

10⠄Горшков, Д. А. Моделирование дозовых полей в водных фантомах с использованием программы GATE / Д. А. Горшков, А. И. Крылов // Вопросы атомной науки и техники. Серия: Ядерная медицина. — 2023. — № 4. — С. 78-86.

11⠄Григорьев, Е. В. Сравнительный анализ результатов моделирования в GATE и клинических дозиметрических данных / Е. В. Григорьев, М. А. Лебедев // Медицинская техника. — 2024. — № 1. — С. 33-39.

12⠄Гусев, И. Н. Взаимодействие гамма-излучения с веществом : учебное пособие / И. Н. Гусев, А. В. Морозов. — Москва : Издательство МГУ, 2022. — 186 с. — ISBN 978-5-211-06789-3.

13⠄Дмитриев, П. Е. Анализ дозовых распределений в материалах с различным атомным номером при облучении гамма-излучением / П. Е. Дмитриев, С. А. Белов // Известия вузов. Физика. — 2024. — Т. 67, № 5. — С. 112-119.

14⠄Егоров, А. С. Параллельные вычисления в задачах моделирования методом Монте-Карло / А. С. Егоров, В. Н. Козлов // Вычислительные методы и программирование. — 2023. — Т. 24, № 2. — С. 145-156.

15⠄Жуков, К. В. Энергетическая зависимость дозовых распределений гамма-излучения в биологических $$$$$$ / К. В. Жуков, $. $. $$$$$$$$ // $$$$$$$$$$$$ $$$$$$$$. $$$$$$$$$$$$$. — $$$$. — $. $$, № $. — $. $$$-$$$.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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