Климанов Дозиметрическое планирование лучевой ч2 2008

δ-электронов с кинетической энергией Te > Temin для материала с плотностью электронов ne рассчитывается из формулы:

где re – радиус электрона; Temin определяет граничную энергию, выше

которой электрон рассматривается как δ-электрон. Сечение равняется нулю при Temax ≤ Temin .

Используя формулу (2.2), можно аналитически определить полное макроскопическое сечение

Полные потери энергии протоном при прохождении пути dz в веществе с плотностью ρ и атомным весом А в результате взаимодействия с электронами описываются формулой Бете – Блоха

(ФББ):

где EI – средний ионизационный потенциал атома/молекулы; Z – атомный зарядовый номер; v – скорость протона; k = 8·π·q2·e4; e – заряд электрона; q – заряд протона; aBarcas – поправка Баркаса (поправка на эффект плотности); ashell – поправка на эффект связи электрона на

оболочке.

Электромагнитное упругое взаимодействие с ядрами атомов служит основной причиной рассеяния протонов. В результате актов рассеяния протоны отклоняются от направления первоначального движения. Угловое распределение таких частиц описывается известной формулой Резерфорда для дифференциального сечения упругого рассеяния:

101

где mnuc и Znuc – масса и заряд ядра соответственно; Zp –заряд протона; η – параметр экранирования; θs – угол рассеяния протона.

Как видно из формулы (2.5) рассеяние протонов происходит в подавляющем большинстве случаев на малые углы. Рассеяние на большие углы, соответствующие малому прицельному расстоянию, является редким событием. Так как большая часть рассеяний происходит на небольшие углы, то для отклонения протона от первоначального направления на заметный угол требуется большое число взаимодействий. Общая теория, описывающая многократное рассеяние, известна как теория Мольера. Распределение Мольера включает несколько членов, первый из которых имеет форму гауссиана. На практике при расчетах часто ограничиваются учетом только первого члена. Это приводит к приближению малых углов и Гауссовскому распределению с шириной, которая согласно Росси и Грейзену после прохождения слоя dz в материале с радиационной длиной X0 равняется

где p – протонный момент; Es – постоянный параметр, не зависящий от энергии протона и состава материала. Его значение определяют подгонкой под результаты точных расчетов.

Образование вторичных электронов в процессе ионизации приводит к флуктуации энергии первичных фотонов. Это явление влияет на крутизну дозового распределения за пиком Брэгга, его необходимо учитывать при расчете дозовых распределений. Данный эффект автоматически учитывается, если образовании δ-электронов

происходит при энергиях выше Temin . Другими словами, при выборе Temin достаточно малым расчет дозового распределения методом Монте-

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

102

флуктуационная функция моделируется гауссианом со средним

где ∆z – геометрический шаг при моделировании траекторий. Похожий подход к учету флуктуации пробега протонов,

обусловленными флуктуациями в потерях энергии, применен в работе [9]. Авторы [9] в приближении теории Ферми – Эйджа решали задачу о пространственном распределении протонов от источника в виде тонкого луча (ТЛ). Они предложили аппроксимировать уменьшение флюенса первичных частиц путем умножения среднего квадрата радиального расширения ТЛ на каждой глубине из-за рассеяния на корректирующий фактор С(z):

где σ2 (R0 ) –дисперсия полной длины пробега для частиц с энергией E0 и пробегом R0;

S/ρ – массовая тормозная способность протонов.

3.2. Ядерные взаимодействия

Вероятность ядерных реакций составляет ~5 % от вероятности ионизационных взаимодействий для 50 МэВ протонов при Temin =0,1

МэВ. С увеличением энергии эта соотношение уменьшается, например для 200 МэВ протонов она меньше 1 %. Поэтому иногда ядерные взаимодействия рассматриваются как поправка к электромагнитным процессам.

Согласно публикации МКРЕ 46 [10] мягкая ткань состоит, в основном, из водорода, углерода, азота и кислорода. В публикации МКРЕ 63 микроскопические сечения взаимодействия протонов для этих элементов, нормированные на атомную массу, имеют примерно одинаковую величину. Поэтому в расчетах часто предполагают, что с точки зрения ядерных реакций протонов вода ведет себя подобно

103

мягкой ткани. Такая замена не совсем справедлива для скелета человека, содержащего 5 – 20 % кальция, так как для кальция величина нормализованных сечений ядерных реакций примерно на 25 % меньше, чем для кислорода. Однако такая разница мало сказывается на дозовых распределениях в протонной терапии, поэтому далее анализируются ядерные взаимодействия протонов только с элементами воды.

3.2.1. Ядерные упругие взаимодействия протон-протон

Макроскопическое поперечное сечение Σpp(Tp) протон-протонного столкновения в воде, нормализованное на плотность ρ, определяется из микроскопического сечения σpp(Tp) по формуле:

где NA – число Авогадро; wH – весовая доля водорода в воде; AH – атомный вес водорода. Значения σpp(Tp) приводятся в ряде работ, например в [12]. Хорошая аналитическая аппроксимация этих данных в диапазоне энергий 10 – 300 МэВ предлагается в работе [5]:

где Tp – энергия протона в МэВ.

Угловое распределение протонов после столкновения можно также оценить, используя работу [12]. Однако в системе центра

распределению энергии протонов в лабораторной системе координат.

3.2.2. Ядерные упругие взаимодействия протон-кислород

Полное макроскопическое сечение упругого столкновения протона с ядром кислорода определяется аналогично выражению (2.10):

где w0 – весовая доля кислорода в воде. Значения σel(Tp) приводятся в публикации МКРЕ 63 [11]. Аналитическая аппроксимация этих данных в диапазоне энергий 50 – 250 МэВ предлагается в работе [5]:

104

где Tp – выражено в МэВ.

Максимальная энергия, передаваемая ядру отдачи, определяется по формуле (2.1) с заменой в ней me на mO. Средняя кинетическая энергия ядер кислорода в зависимости от энергии протона табулирована в работе [11]. Аналитическая аппроксимация, предлагаемая в работе [5] имеет вид:

TΟ = 0,65exp(−0,0013 Tp ) −0,71exp(−0,0177 Tp ). (2.14)

Распределение ядер кислорода по энергии TO после упругого столкновения с протоном имеются в работе [12]. Эти данные аппроксимируются в работе [5] следующей формулой:

Пробег ядер отдачи кислорода в воде очень мал, поэтому считается, что энергия, переданная этому ядру, поглощается локально. Направление движения протона после такого столкновения изменяется

всоответствии с переданной энергией.

3.2.3.Ядерные неупругие взаимодействия протонкислород

При неупругих столкновениях протонов с ядром кислорода в рассматриваемой области энергий возможны следующие основные реакции:

Полное макроскопическое поперечное сечение неупругого взаимодействия определяется из микроскопического сечения

105

где значения σin(Tp) можно взять из работы [12]. Аналитическая аппроксимация этих данных в диапазоне энергий 7 – 250 МэВ, предложенная в работе [5], имеет вид:

Кроме вторичных протонов и нейтронов в этих взаимодействиях могут рождаться также дейтроны, тритоны, альфа-частицы, тяжелые фрагменты, фотоны, электроны и т.д. Полный учет всех каналов реакции проводится в таких сложных программах ядерной физики, как GEANT4 [13], MCNPX [14] и др. Однако такое детальное рассмотрение вряд ли нужно в области протонной терапии. Согласно данным публикации МКРЕ 63 [12] и библиотеки сечений [15] приблизительно 50 % энергии первичного протона при неупругом столкновении передается вторичному протону. Значительная доля энергии переходит к нейтронам, дейтронам, альфа-частицам и тяжелым ядерным фрагментам. Но из этих частиц только нейтроны распространяются на большие расстояния. Остальные частицы можно считать поглощаемыми локально. На рис. 2.4 показаны вклады в дозу, которые создаются разными каналами реакций при первичном взаимодействии 150 МэВ протонов в воде [5]. Например, кривая, обозначенная как “нейтроны”, представляет вклад в дозу от всех частиц, создаваемых нейтронами, которые образовались при неупругих ядерных реакциях первичного протона с ядрами кислорода.

3.3. Массовая тормозная способность

Полная массовая тормозная способность протонов представляет собой сумму электронной тормозной способности (см. приложение табл. П.5) и ядерной тормозной способности. При кинетической энергии выше 1 МэВ ядерная тормозная способность пренебрежимо мала, т.е. протоны замедляются благодаря, в основном, неупругому взаимодействию с атомными электронами.

Для ткани человека можно выразить тормозную способность материала с плотностью ρ через функцию от массовой тормозной способности для воды Sw(Tp)/ρw:

106

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

Рис. 2.5. Отношения массовых тормозных способностей протонов с энергией 10 и 100 МэВ для разных материалов к массовой тормозной способности воды, рассчитанные из данных работы [10]

ивычисленные по формуле (2.20)

3.4.Ограниченная массовая тормозная способность

108

в котором происходит торможение протонов. К недостаткам циклотрона относится и большое потребление электроэнергии, так как он работает в непрерывном режиме. На рис. 2.6 приводится фотография циклотрона, производимого для ПЛТ компанией IBA (Бельгия).

Рис. 2.6. Циклотрон компании IBA на максимальную энергию 235 МэВ диаметром 4 м и электропотреблением 400 кВт [4]

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

110