Основные результаты получены по следующим четырем тесно связанным направлениям проекта:
1) теоретические оригинальные идеи
=========================
1. Внедрение N-нуклонных конфигураций (флуктонов) и верификация МК-генераторов
Проведен анализ выполнения законов сохранения (энергии, электрического, барионного и лептонного зарядов) в генераторах EPOS4, UrQMD, PHQMD и PYTHIA для энергий NICA. Выявлены существенные нарушения законов сохранения в EPOS4 и PHQMD (зависящие от прицельного параметра), что ограничило их применение для изучения подпороговых эффектов.
Модифицирован исходный код генератора UrQMD: часть пар нуклонов заменена на компактные двухнуклонные (флуктоны F) с барионным числом B=2. Введено три типа флуктонов (NN, PN, PP). Реализован алгоритм распада флуктона. Введен новый класс столкновений F+N→N+N+N.
Проведено моделирование подпорогового и околопорогового рождения антипротонов и кумулятивных процессов в ядро-ядерных столкновений.
Произведена начальная настройка параметров, связанных с флуктонов. Сделан обзор литературы по подпороговому рождению и кумулятиву
(с точки зрения подходящих экспериментальных данных).
Ведется работа по внедрению флуктонов в генератор SMASH.
2. Развитие модели мультипомеронного обмена (ММПО) и механизмов рождения странности
В рамках модели мультипомеронного обмена вычислена величина эффективного натяжения струны, проведено успешное сравнение с феноменологическими подходами, основанными на анализе экспериментальных данных.
Модифицирован подход для учета флуктуаций эффективного натяжения струны вокруг его среднего значения (с учетом требований КХД-конфайнмента). Показано, что учет этих флуктуаций оказывает существенное влияние, особенно для мультистранных частиц, и позволяет улучшить описание выходов в ядро-ядерных столкновениях в широком диапазоне энергий.
Произведена модификация модели ММПО, включающая подход Galoyan-Ribon-Uzhinsky для расчета вероятностей рождения странных частиц в ограниченном фазовом объеме для странных и мультистранных частиц. Проведена верификация на экспериментальных данных pp и AA столкновений в диапазоне энергий NICA/SPS.
Проанализированы отношения средней плотности поперечной энергии для мезонов (φ, K, π) в диапазоне √sNN = 39 ГэВ – 5.44 ТэВ. Подтверждена гипотеза об универсальности механизмов рождения частиц с различным кварковым составом.
3. Гидродинамическое моделирование столкновений ядер с голографическим уравнением состояния (УрС)
Разработан метод расширения голографической термодинамики за счет учета химических потенциалов странности и электрического заряда. Свободные параметры модели откалиброваны с помощью методов машинного обучения (ML) на основе данных решеточной КХД (LQCD). Разработанное голографическое УрС интегрировано в гидродинамические коды MUSIC и vHLLE. Выполнено численное моделирование полного цикла эволюции (фреймворки iEBE-MUSIC и SMASH-vHLLE).
Рассчитаны спектры по быстроте, поперечному импульсу и поперечной массе для каонов, пионов, протонов и лямбда-барионов при энергиях SPS и RHIC, проведено сравнение с экспериментальными данными.
4. Моделирование эффектов нарушения изотопической симметрии.
Проанализированы данные, опубликованные экспериментом NA61/SHINE и другими коллаборациями по отношению выходов каонов (заряженных к нейтральным) R_K и выяленным статистически значимым отклонением от единицы (значимость 4.7σ – 5.3σ) в диапазоне энергий NICA-SPS, что соответствует нарушению изотопический симметрии в сильных взаимодействиях. Проанализирован и сделан обзор имеющихся в литературе различных механизмов и подходов, которые могут давать вклад в этот эффект (массы кварков, электромагнитные взаимодействия, пороговые явления при распадах, дисбаланс рождения u и d кварков). В рамках модифицированной модели UrQMD, основанной на параметризованном нарушении симметрии при разрыве струн (P(u):P(d)=3:1), проведено моделирование нарушения изотопической симметрии в столкновении Ar+Sc, Xe+Xe, Xe+W, Bi+Bi, C+C при энергиях соответствующих российским ускорителям NICA, У70, У1000.
Произведена оценка требований к статистике и разрешению детектора MPD для исследования явления в будущих экспериментальных данных.
5. Моделирование локального нарушения P-четности в сильных взаимодействиях.
Выполнено исследование условий для локального нарушения пространственной четности, генерации топологического заряда и кирального дисбаланса, оценка наблюдаемых сигнатур.
Обосновано применение методов парциально-волнового анализа (ПВА) для выделения сигнала расщепления поляризаций в дилептонных каналах распада легких векторных мезонов (ρ,ω→e+e−,μ+μ−). Определены оптимальные кинематические переменные (инвариантная масса, pT, углы вылета в системе Готфрида-Джексона) для эффективного разделения вкладов различных поляризационных состояний на фоне низкого уровня комбинаторного адронного фона.
Сформулировано предложение по поиску эффектов расщепления поляризаций в адронных распадах с использованием многочастичных корреляционных методов, однако показано, что они осложнены необходимостью количественной оценки искажений угловых распределений вследствие адронных перерассеяний на поздних стадиях эволюции среды. Продемонстрирована необходимость привлечения микроскопических транспортных моделей для корректного учета фоновых процессов.
Произведено обоснование оптимального энергетического диапазона для NICA для данного типа исследования. Проведен критический анализ энергетической зависимости эффектов кирального дисбаланса. Обосновано, что сверхвысокие температуры (LHC) приводят к тепловому и магнитному экранированию CP-нечетных доменов. В то же время, среда с высокой барионной плотностью стабилизирует локальные домены и максимизирует плотность топологических флуктуаций. Сделан фундаментальный вывод о приоритетности поиска данных эффектов на коллайдерах средних энергий (NICA, SPS, FAIR, RHIC BES, У70, У1000), где достигается оптимальный баланс для регистрации эффектов локального несохранения четности в сильных взаимодействиях. Полученные теоретические результаты и оценки формируют методическую базу для будущих анализов данных детекторов MPD и SPD (NICA).
6. Исследование сильно-интенсивных переменных и минимизация объемных флуктуаций
Разработан систематический подход к конструированию сильно-интенсивных переменных, устраняющих зависимость от флуктуаций числа источников в рамках модели независимых источников. На основе анализа совместных кумулянтов наблюдаемых множественностей и свойств эмиссии единичного источника аналитически вычислены независимые комбинации второго (Σ, Δ), третьего (Γ) и четвертого (Υ ) порядков. Особое внимание уделено наблюдаемой 4-го порядка Υ (Upsilon) и строгому определению условий, при которых флуктуации числа источников полностью сокращаются.
Построена сильно-интенсивная переменная для анализа вперед-назад быстротных корреляций множественностей. Проведен анализ влияния конфигураций цветовых струн на эти корреляции, что закладывает теоретическую основу для получения предсказаний переменных Γ и Upsilon
для энергий LHC (эксперимент ALICE).
Обоснована критическая проблема искажения кумулянтов высших порядков (негауссовых флуктуаций сохраняющихся величин: электрического заряда, странности, барионного числа) тривиальными объемными флуктуациями. Показано, что использование множественности заряженных частиц для отбора классов центральности существенно искажает сигналы критической точки. Предложен и теоретически обоснован метод минимизации этих эффектов путем перехода от множественности заряженных частиц к числу нуклонов-участников (Npart) при определении начальных условий столкновения.
7. Разработка модели d-d взаимодействия в рамках классического глауберовского подхода
Дейтрон-дейтронное взаимодействие исследовалось в рамках глауберовского подхода с учётом отличия от случая столкновений тяжелых ядер (в частности, симметричное расположение нуклонов относительно центра масс дейтрона). Получены оценки влияния свойств взаимодействия на нуклонных масштабах на флуктуации числа нуклонов-участников и числа нуклонных столкновений. Разработан алгоритм монте-карловского (МК) моделирования процесса рассеяния дейтронов. Были получены выражения для среднего значения и дисперсии числа раненых нуклонов, дисперсии общего числа нуклонов-участников, среднего значения и дисперсии числа нуклонных столкновений для dd взаимодействия при фиксированном значении прицельного параметра. Также были оценены вероятности двух последовательных столкновений нуклона из одного дейтрона с двумя нуклонами другого дейтрона и вероятности взаимодействия двух нуклонов одного дейтрона с двумя нуклонами другого при
фиксированном значении прицельного параметра.
8. Развитие грид-инфраструктуры СПбГУ и интеграция с ОИЯИ
Обеспечена стабильная интеграция в GRID ОИЯИ в рамках эксперимента SPD уровня продакшн на кластере СПбГУ (что подтверждено в коллаборации). Ведется подготовка стораджа СПбГУ для работы в распределенной СХД ОИЯИ (EOS). Обеспечено функционирование внутренней сети кластера на полной скорости (10 Гбит). Проведено масштабное сетевое тестирование внешнего 10-гигабитного канала (до ОИЯИ) и подтверждены его скоростные характеристики в многозадачных нагрузках. Реализована авторизация и аутентификация удаленных пользователей на основе не только X.509 сертификатов, но и токенов в трех экспериментах ОИЯИ. Кластер успешно подключен к грид-инфраструктуре MPD+BM@N (основанной на middleware DIRAC), обеспечена стабильная одновременная работа с экспериментом SPD и локальными пользователями.
2) цифровые технологии (включающие в себя МК, МО и ИИ),
=======================================
В отчётном периоде выполнены работы по разработке и верификации
нейросетевого подхода к событийной оценке прицельного параметра (b) в
столкновениях тяжёлых ионов. Результаты, достигнутые в соответствии с
поставленными задачами этапа НИР, представлены ниже.
Задача 1. Подготовка модельных данных и моделирование детекторной системы.
Сформированы модельные данные о столкновениях Au+Au при энергии √s_NN =
11 ГэВ, полученные с помощью двух независимых генераторов событий — QGSM
и EPOS — по 200 000 событий на каждый генератор. Диапазон прицельного
параметра 0 < b < 16 фм; распределение событий по b выбрано
реалистичным, с учётом геометрической зависимости вероятности
столкновения. Из событий извлекалась информация о вторичных заряженных
частицах и их импульсах; при моделировании ограниченного акцептанса
учитывалось попадание частиц в геометрию детектора. В качестве модели
детекторной системы рассматривалась конфигурация на основе
микроканальных пластин (МКП), рассматриваемая для применения в
экспериментах на коллайдере NICA: два кольца по разные стороны от точки
столкновения (внешний радиус 25 см, внутренний радиус 2,5 см, расстояние
от точки столкновения 4 м, временное разрешение 50 пс).
Задача 2. Оценка прицельного параметра по полной (4π) информации о событии.
Исследовано качество восстановления прицельного параметра при
доступности полной информации о событии в геометрии 4π, что позволяет
оценить предельное качество восстановления и сравнить различные
представления данных. Рассмотрены три представления события: (i)
двумерное изображение 40×40 пикселей — проекция углового распределения
частиц на плоскую сетку; (ii) массив характеристик отдельных частиц
3×1600 (время пролета, косинус угла вылета, знак продольной компоненты
импульса); (iii) набор из 10–20 физических наблюдаемых (множественность,
угловые и времяпролетные распределения, передне-задняя асимметрия). В
качестве моделей применялись сверточные нейронные сети (1–2 сверточных
слоя + полносвязные слои), архитектура ResNet, полносвязные сети и
архитектура «deep reconstruction neural network» (DRNN). Наилучший
результат (RMSE 0,49 фм для QGSM и 0,58 фм для EPOS) получен сверточной
сетью для представления 3×1600. Для оценки качества применялась метрика
RMSE с определением статистической неопределённости bootstrap-методом
(10 000 выборок; неопределённость RMSE порядка 0,01–0,02 фм)
Задача 3. Исследование устойчивости алгоритмов к выбору модели генерации
событий.
Показано, что при обучении сети на данных одного генератора и
тестировании на другом качество восстановления существенно ухудшается:
для представления 3×1600 ошибка возрастает с 0,50 до 1,68 фм (обучение
на QGSM, тест на EPOS) и составляет 1,43 фм в обратном направлении;
аналогичная зависимость наблюдается для других представлений. Это
демонстрирует, что высокая точность внутри одного модельного набора не
гарантирует переносимости алгоритма на экспериментальные данные.
Предложено и проверено решение: смешивание событий двух генераторов в
обучающем наборе (соотношение 1:1) в сочетании с архитектурой DRNN
позволяет существенно снизить чувствительность к выбору генератора —
RMSE составляет 0,66 фм (тест на QGSM) и 0,62 фм (тест на EPOS), т.е.
качество восстановления на независимых модельных наборах становится близким.
Задача 4. Построение физически интерпретируемого промежуточного
представления события.
Показано, что набор физических наблюдаемых содержит значительную часть
информации, необходимой для восстановления прицельного параметра: при
включении множественности RMSE составляет 0,50 фм (QGSM) и 0,61 фм
(EPOS), при её исключении — 0,66 и 0,83 фм соответственно. Проведён
SHAP-анализ значимости признаков: множественность заряженных частиц
является доминирующим признаком (относительный вклад 0,68 для QGSM и
0,64 для EPOS), остальные наблюдаемые (времяпролетные характеристики,
угловые распределения, асимметрии) обеспечивают дополнительную
информацию. Важно, что общая структура значимости наблюдаемых
сохраняется для двух независимых моделей генерации, что подтверждает
физическую обоснованность предложенного представления.
Задача 5. Проверка двухэтапной схемы оценки прицельного параметра по
данным детектора с ограниченным акцептансом.
Разработана и проверена двухэтапная процедура: на первом этапе нейронная
сеть восстанавливает по данным детектора (вход — массив 3×150: время
пролета, номер ячейки, номер кольца) десять физических наблюдаемых
(полная множественность заряженных частиц, числа частиц в передней и
задней полусферах, абсолютная передне-задняя асимметрия, средний косинус
угла вылета, среднее и среднеквадратичное время пролета, 10-й, 50-й и
90-й процентили распределения времени пролета); на втором этапе
прицельный параметр оценивается по реконструированным (а не истинным)
наблюдаемым. Качество первого этапа: RMSE полной множественности 67
частиц (QGSM) и 174 частицы (EPOS), среднего косинуса угла 0,066 и 0,047
соответственно. Показано, что введение промежуточного физически
интерпретируемого слоя не приводит к существенной потере точности по
сравнению с прямым восстановлением: 0,69 фм против 0,68 фм (QGSM) и 1,60
фм против 1,53 фм (EPOS) при обучении/тестировании на отдельных наборах;
для смешанных данных — 0,73/1,60 фм (двухэтапная схема) против 0,73/1,61
фм (прямой метод).
Задача 6. Исследование N-нуклонных кластеров ядерной материи6.1 Внедрение N-нуклонных конфигураций в процесс генерации ядер в Монте-Карлогенератор событий UrQMD4. Настройка соотношения вероятностей образования 2N конфигураций по сравнению с имеющимися данным по кумулятивному рождению пи-мезонов и протонов
В генератор событий UrQMD4внедрена модель кумулятивного рождения частиц на основе флуктонного механизма.Для этого добавлены новые частицы — флуктоны F,а также каналы взаимодействия F+N и F+F, с сохранением барионного числа, заряда иэнергии-импульса в каждом узле. Состояния корреляций формируются с измеренным изоспиновым составом (np:pp:nn≈90:5:5), что согласуется с экспериментальными данными по короткодействующим корреляциям.
Модель настроена по данным для p+C при 12 ГэВ/с. В столкновениях Au+Auпри √sNN = 4.5 и √sNN = 8 ГэВ в области x>1.7 модифицированная модель показывает увеличенное количество кумулятивных протонов по сравнению со стандартной версией UrQMD. Этот избыток растёт с x и достигаетобласти x = 2.5–3.2. Это — ожидаемая сигнатура многонуклонных корреляций,отсутствующая в стандартной версии UrQMD.
Мы продемонстрировали, чтомодификация UrQMD с введением в состав ядер флуктонов дает возможность получения кумулятивныхчастиц в области высоких значений кумулятивных переменных.
6.2 Моделирование столкновения ядер Xe и W при энергиях NICA с использованием PHQMD, UrQMD4 и UrQMD4-mod
Предложенная намим одификация UrQMD4-mod, учитывающая наличие в ядрах флуктонов, увеличивает выход антибарионов (ˉp, ˉn, ˉΛ, ˉΣ) в подпороговой области в 2–3 раза.
В генераторе PHQMDобнаружено нарушение сохранения энергии, величина ошибки зависит от прицельного параметра.
Выход ˉΛ в UrQMD4,UrQMD4-mod и в PHQMD примерно на порядок меньше выхода ˉp и ˉn.
Различия в изменениисуммарной множественности протонов и нейтронов в UrQMD4, UrQMD4-mod и PHQMD обусловлены разной динамикой рождениягиперонов, в основном Λ и Σ.
Результаты моделированиявыхода Λ-гиперонов с помощью модели UrQMD4-mod, учитывающей наличие в ядрах флуктонов, по-видимому, лучше соответствует результатам экспериментов BM@N, где предсказания транспортных моделей (DCM-SMM, UrQMD и PHSD) систематически превышают измеренные выходы. Поскольку PHQMD дает заметно завышенный выход посравнению с UrQMD 4.0, выход Λ-гиперонов в этом генераторе, по-видимому, сильно завышен по сравнению с экспериментальными данными.
Часто используемое предположение о равенстве среднего выхода протонов и нейтронов ⟨p⟩=⟨n⟩ в диапазоне энергий NICA не выполняется. При этом PHQMD, UrQMD4 и UrQMD4-mod дают близкие значения для ⟨p⟩ илидля ⟨n⟩ вобласти √sNN от 2.6 до 5ГэВ, но их предсказания начинают сильно расходиться при увеличении √sNN до 9 ГэВ.
6.3 Изучение нарушенияизотопической симметрии при столкновении ядер
Представлены результаты Монте-Карло моделирования в рамках стандартной и модифицированной версий генератора UrQMD4 для столкновений 40Ar+45Sc, 12C+12C,продемонстрировано преимущество симметричных систем для независимой проверки данного эффекта. Произведена оценка потенциала детектора MPD на коллайдере NICAдля исследования данного явления в столкновениях d+d, Xe+Xe, Bi+Bi, Xe+W. Сформулированытребования к статистике, разрешению по импульсу и идентификации частиц, а также возможные стратегии анализа данных для верификации нарушения изоспиновой симметрии в новых энергетических диапазонах.
6.4 Корректировка формул C, P, T симметрий для киральных полей
Операции зарядового сопряжения С, пространственной инверсии (четности) Р иинверсии времени Т играют важную роль в современной квантовой теории поля. Однакоформулы для операторов C, P и T были получены в период зарождения квантовойтеории поля еще до открытия киральной асимметрии слабых взаимодействий. Мы показали,что они либо неприменимы для киральных операторов, либо содержат ошибки вформулах или их интерпретации. Показано, что стандартные формулы CPT сопряжения неприменимы в кирально-несимметричных теориях. Исправлены квантовополевые формулы преобразований Лоренца и C, P, T сопряжений спиноров. Они обеспечивают правильную унитарность или антиунитарность дискретных преобразований и киральность результатов.
Научная и практическая значимость по направлению 2.
Полученные результаты показывают перспективность физически
интерпретируемого нейросетевого подхода к определению параметров
столкновений тяжёлых ионов: предложенная схема связывает ограниченную
информацию детектора с физическими характеристиками события и позволяет
независимо исследовать и сопоставлять восстановленные величины с
экспериментальными данными, что критически важно при переходе от
модельных данных к данным экспериментов (в частности, на коллайдере
NICA). Метод не требует явной реконструкции траекторий и идентификации
частиц, что расширяет возможности анализа событий с использованием
детекторов на МКП.
Исследование с применением технологий искусственного интеллекта
выполнено в области физики высоких энергий / ядерной физики
(физика столкновений тяжёлых ионов), где технологии искусственного
интеллекта выступают в качестве инструмента для решения конкретной
исследовательской задачи — событийной оценки прицельного параметра
столкновений тяжёлых ионов, который не может быть измерен напрямую и
восстанавливается по признакам события, связанным с его начальной
геометрией (множественность заряженных частиц, угловые и времяпролетные
распределения).
Для решения задачи используются уже имеющиеся (стандартные) технологии ИИ:
сверточные нейронные сети (CNN) с одним-двумя сверточными слоями и
полносвязными слоями;
архитектура ResNet для анализа двумерных представлений событий;
полносвязные нейронные сети (MLP) для работы с наборами физических
наблюдаемых;
архитектура «deep reconstruction neural network» (DRNN) для
построения представления события, устойчивого к изменению источника
данных (модели генерации);
метод SHAP для количественной оценки значимости признаков
(физических наблюдаемых) при восстановлении скрытого параметра события.
Результат исследования вносит вклад в область ядерной и физики высоких
энергий: предложен и проверен на модельных данных двух независимых
генераторов (QGSM, EPOS) физически интерпретируемый нейросетевой метод
восстановления параметров столкновений (прицельного параметра) по данным
детектора с ограниченным акцептансом на основе микроканальных пластин,
рассматриваемого для экспериментов на коллайдере NICA (Дубна).
Результаты могут быть использованы при планировании и анализе
экспериментов по изучению свойств ядерного вещества в экстремальных
условиях (определение центральности событий, геометрии столкновений) и
представляют интерес для смежных задач анализа больших массивов данных
физики высоких энергий.
3) детекторные технологии,
===================1. Разработан малошумящий детектор на МКП, предназначенный для регистрации редких процессов с вылетом частиц в кинематически запрещенную область, который предалагается использовать в качестве триггера для регистрации кумулятивных частиц в ядро-ядерных столкновениях на пучках коллайдера NICA. 2. Представлена оптимизированная конструкция, которая обеспечивает устойчивость к синфазным электромагнитным наводкам. 3. Выполнены проверка и настройка работы четырехканального блока электроники на быстрых компараторах, что позволило обеспечить высокую эффективность регистрации космических лучей (на уровне 95+-5 %.
4) практические приложения для диагностики -- протонная томография (pCT)
==================================================.
4.1. Создание общей концепции протонноготомографа
Проведен анализ имеющихсяна сегодняшний день разработок в области протонной томографии. Основной задачейпротонного томографа является получение трёхмерного изображенияпространственного распределения тормозной способности опухолевого образованиявнутри тела пациента по отношению к протонам. Путем измерений координат и углов,под которыми протон входит в тело пациента, и соответствующих углов и координатна выходе из тела, а также при знанииэнергии протона, можно идентифицировать наиболее вероятная траектория и потери энергиипротонов вдоль неё. Проделав эти измерения для множества протонов, можнореконструировать пространственное распределение тормозной способностиоблучаемого объекта. На первом этапе работ был проанализирован опыт группыуниверситета Бергена по созданию протонного томографа на основе монолитныхактивных пиксельных сенсоров (МАПС - детекторы). В итоге, в СПбГУ былапредложена более совершенная концепция, включающая в себя передние трековыеплоскости (перед пациентом), задние трековые плоскости (после пациента) ицифровой трековый калориметр (ЦТК). Все трековые плоскости будут содержатьдетекторы МАПС, работа с которыми интенсивно ведется в СПбГУ. Такая концепцияпозволить с большой точностью проследить треки отдельных протонов и определятьих остаточную энергию, после прохождения тканей и органов пациента. Такимобразом, будет получена информация об относительных энергетических потеряхпротонов в теле пациента, что будет способствовать точному определению энергиипротонного пучка ускорителя в методах протонной терапии. Планируемая ошибка вопределении максимальной дозы на опухоли (пик Брэгга) составит порядка 0.5 мм(современные диагностические методы компьютерной томографии дают ошибки в 6-20мм).
4.2. Создание модели прототипа протонноготомографа, состоящего из трековых плоскостей и цифрового трекового калориметра.Моделирование трекового калориметра протонного томографа.
В рамках работ по созданию протонноготомографа на основе детекторов МАПС была разработана программа теоретических иэкспериментальных исследований. Была создана модель цифрового трековогокалориметра (ЦТК) для протонного томографа. ЦТК включал в себя слои Монолитных АктивныхПиксельных Сенсоров (МАПС), созданных по технологии КМОП, чередующихся со слоямипоглотителей из различных материалов (выбор материала и его толщины определялсяпутем моделирования). Высокая степень гранулярности пиксельных сенсоровпозволяет регистрировать и разделять большое количество протонных треков вкаждом цикле считывания данных, а их рекордное быстродействие удовлетворяеттребованию протонной томографии с высокой частотой сканирования ( >10 МГц),что сводит к минимуму время облучения пациента. Для предложенной концепции ЦТК, былопроведено моделирование транспорта протонов через его слои. На основе анализатранспорта протонов с энергиями 150 МэВ и 200 МэВ (выбранные наиболее частоиспользуемые энергии из диапазона работы современных ускорителей, используемыхв протонной терапии онкологических заболеваний) через слои калориметра былиоптимизированы геометрические параметры предложенной модели. Это позволиловыбрать материал поглотителя цифрового трекового калориметра и оптимизироватьего толщину. В итоге оптимальным оказался никелевый поглотитель с толщиной 1.5мм. С использованием данного поглотителя была показана экономия в МАПСдетекторах, используемых при построении протонного томографа. Другимоптимальным материалом с точки зрения дешевизны и минимума дополнительнойрадиационной нагрузки на пациента, оказался алюминий толщиной 3.5 мм. Был также разработан алгоритмреконструкции треков на основе клеточного автомата. Данный алгоритм сочетает всебе высокую эффективность восстановления прямых треков и высокое быстродействиев их реконструкции. Применяя указанный алгоритм, получили высокую эффективностьреконструкции протонных треков (90%) для широкого пучка диаметром 20 мм. Прискорости считывания МАПС ~5-10 мкс-1 для достижения такого количества частиц всобытии потребуется интенсивность протонного пучка ~10^6 с-1. Для узкого«карандашного» пучка (диаметр 2 мм), а также при более высокой интенсивностипротонного пучка (100 протонов в событии) эффективность алгоритма не будетпревышать 50%, а время обработки может составлять несколько часов. Важный практическийвывод по результатам данного моделирования – использование протонных пучков сширокой апертурой.
На следующем этапе исследований, было выполнено моделирование всегопротонного томографа с учетом трековых плоскостей и ЦТК, с применением веществвходящих в состав организма человека. Это так называемые фантомы: кости, жир,мозг, вода, воздух. Диапазон энергий протонов был расширен и составил от 125МэВ до 250 МэВ (граничные энергии современных ускорителей, используемых впротонной терапии онкологических заболеваний). Был разработан алгоритм повосстановление начальной энергии протонов в ЦТК, были получены карты тормозныхспособностей протонов в веществе фантомов: кости, жир, мозг, вода, воздух. В итоге былаполучена остаточная энергия после прохождения указанных фантомов и оцененоэнергетическое разрешение этого метода. Варьирование энергии ускорителя позволило найти нижний предел поразрешению ЦТК. Также было проведено сравнение используемого метода с методомпоглощения протонов в веществе фантомов и методом отклонения треков в веществефантомов. Было показано, что методика с использованием цифрового трековогокалориметра наиболее оптимальна.
4.3. Экспериментальные и теоретическиеисследования с использованием телескопа монолитных активных пиксельных сенсорови гибридных пиксельных детекторов в контексте задач протонной томографии.
В рамках работ для проведения экспериментальных исследований, былиспользован стенд, в который входит детекторный телескоп на основе монолитныхактивных пиксельных сенсоров (МАПС). Данный экспериментальный стенд применялсядля оценки идентификационной способности цифрового трекового калориметра (ЦТК). Стенд использовался дляпроведений экспериментальных исследований с использованием телескопа МАПС напучке протонов с энергией 200 МэВ синхроциклотрона СЦ-1000 в НИЦ КИ ПИЯФ. Витоге, впервые были получены карты сработавших пикселей (хитов) сиспользованием телескопа МАПС для указанных энергий протонов. Также былпроведен кластерный анализ с оценкой средней множественности сработавшихпикселей, формирующих конкретный пиксельный кластер. Важность этих исследованийобусловлена тем, что по множественности пикселей и соответственно идентификацииобразовавшихся кластеров можно определить энергию протонов, оставленнуюнепосредственно в детекторном слое. Была выполнена калибровка по энергии взависимости от среднего числа пикселей, входящих в пиксельный кластер. Этоувеличивает точность в определении остаточной энергии протонов. Был также созданэкспериментальный модуль на основе гибридных пиксельных детекторов. Этот модульбудет использован на пучках протонов энергий от 30 МэВ до 250 МэВ, вэкспериментах в Медицинском радиологическом научном центре им. А. Ф. Цыба (г.Обнинск), в рамках решения задач протонной томографии.
Дополнительныеэкспериментальные работы
Поскольку детекторы МАПС будут смонтированы на специальныхуглекомпозитных панелях поддержки и охлаждения (разработки СПбГУ, патентыСПбГУ), то необходимо было провести исследования тепловыхсвойств таких углекомпозитной панелей, для определения эффективности охлаждениядетекторов МАПС. Для этих целей был создан специальный экспериментальный стенд.Стенд включал в себя пространство из материала AIREX, внутри которого находитсярадиатор охлаждения, нагреватель с исследуемой углекомпозитной панелью итермосенсорами. Углекомпозитная панель охлаждалась только с одного края спомощью специального чиллера. Была проведена серия измерений с подачей тепловоймощности на углекомпозитную панель и подачей соответствующего охлаждения.Тепловая мощность подбиралась в соответствии с тепловыделением отдельных системи электронных компонент детекторов МАПС. Были получены и проанализированыраспределения температуры от мощности наквадратный сантиметр. Важным результатом стало, что при мощности 48 мВт/см2(сегодняшние детекторы МАПС имеют мощность в районе 40 мВт/см2) охлаждение в предложенномварианте работает не очень эффективно, поскольку максимальная температура на углекомпозитнойпанели зафиксирована в отметке +47 С, а рабочая температура детекторов МАПСсоставляет +25 С – +30 С. В качестве рекомендаций, было предложено охлаждать углекомпозитную панель с двух сторон. Измеренияпоказали, что в данном случае максимальная температура будет в центреисследуемого образца, но она не будет превышать отметки в +30 С.
Дополнительныеэкспериментальные работы
Работы по изучениюпроцессов вторичной электронной эмиссии. Для изучения указанных процессов былаиспользована разработанная в СПбГУ многосенсорная система мониторинга пучковзаряженных частиц и тяжелых ионов (патент СПбГУ). Система предназначена длявизуализации профилей пучков протонов, дейтронов, альфа-частиц и тяжелых ионовв широком диапазоне энергий. Принципдействия системы основан на эффекте вторичной электронной эмиссии. Былипроанализированы результаты экспериментов, выполненных на Уникальной научнойустановке (УНУ) "Физико-технический институт имени А.Ф. Иоффе циклотронтипа У-120" с использованием пучковпротонов с энергией 3 МэВ и пучков ионов 40Ar+8 с энергией 53 МэВ. В ходеэкспериментов были визуализированы профили пучков, а также исследованыхарактеристики сигналов, формируемых сенсорами системы при взаимодействии счастицами пучка. Обработка экспериментальных данных позволила определитькоэффициенты вторичной электронной эмиссии золотого покрытия — материала сенсоров.Полученные результаты стали основой для проверки теоретических моделейпроцессов вторичной электронной эмиссии при взаимодействии вещества сзаряженными частицами и тяжёлыми ионами низких и средних энергий.
4.4. Сравнительный анализ результатовдиагностических исследований на имеющихся рентгеновских томографах ирекомендации для разрабатываемой концепции протонного томографа
Первым направлениемнаучных изысканий стало изучение доз облучения при исследованиях с помощью ужесуществующих компьютерных томографов (для различных областей человеческоготела). Для поиска и сбора информации были изучены как данные современнойнаучной литературы, так и информация, полученная от коллег-рентгенологов,работающих в различных медицинских учреждениях СПб. К настоящему моменту «вилка»доз облучения достаточно большая, так как находится под влияниям целого рядафакторов. Таким образом, сбор и обобщение информации продолжается.
Ответственныеза каждый пункт (курсивом выделены участники, которые не входят в основныеисполнители проекта):
1. В.И. Жеребчевский, В.В. Петров, Н.А.Мальцев
2. В.И. Жеребчевский, В.В. Петров, Н.А.Мальцев
3. В.И. Жеребчевский, В.В. Петров, Н.А.Мальцев, Е.О. Землин, Н.А. Прокофьев,А.О. Скалецкий
4. Д.А. Пую