2025-11-13T23:34:11.432936

Reinforcement learning-based statistical search strategy for an axion model from flavor

Nishimura, Miyao, Otsuka
We propose a reinforcement learning-based search strategy to explore new physics beyond the Standard Model. The reinforcement learning, which is one of machine learning methods, is a powerful approach to find model parameters with phenomenological constraints. As a concrete example, we focus on a minimal axion model with a global $U(1)$ flavor symmetry. Agents of the learning succeed in finding $U(1)$ charge assignments of quarks and leptons solving the flavor and cosmological puzzles in the Standard Model, and find more than 150 realistic solutions for the quark sector taking renormalization effects into account. For the solutions found by the reinforcement learning-based analysis, we discuss the sensitivity of future experiments for the detection of an axion which is a Nambu-Goldstone boson of the spontaneously broken $U(1)$. We also examine how fast the reinforcement learning-based searching method finds the best discrete parameters in comparison with conventional optimization methods. In conclusion, the efficient parameter search based on the reinforcement learning-based strategy enables us to perform a statistical analysis of the vast parameter space associated with the axion model from flavor.
academic

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

Основная информация

  • ID статьи: 2409.10023
  • Название: Reinforcement learning-based statistical search strategy for an axion model from flavor
  • Авторы: Satsuki Nishimura, Coh Miyao, Hajime Otsuka (Kyushu University)
  • Классификация: hep-ph (Феноменология физики высоких энергий), cs.LG (Машинное обучение), hep-th (Теория физики высоких энергий)
  • Дата публикации: arXiv:2409.10023v2 hep-ph 11 октября 2025
  • Ссылка на статью: https://arxiv.org/abs/2409.10023

Аннотация

В данной работе предложена стратегия поиска на основе обучения с подкреплением для исследования новой физики за пределами Стандартной модели. Обучение с подкреплением как один из методов машинного обучения является мощным инструментом для поиска параметров модели, удовлетворяющих феноменологическим ограничениям. В качестве конкретного примера авторы сосредоточились на минимальной аксионной модели с глобальной симметрией U(1) флейвора. Обучаемый агент успешно нашел распределения U(1)-зарядов кварков и лептонов, решающие проблемы флейвора и космологии в Стандартной модели, обнаружив более 150 реалистичных решений для кваркового сектора с учетом эффектов перенормировки. Для решений, найденных методом обучения с подкреплением, авторы обсудили чувствительность будущих экспериментов по обнаружению аксионов — бозонов Намбу-Голдстоуна спонтанно нарушенной симметрии U(1). Авторы также проверили скорость метода поиска на основе обучения с подкреплением по сравнению с традиционными методами оптимизации при поиске оптимальных дискретных параметров.

Исследовательский контекст и мотивация

Определение проблемы

  1. Нерешённые проблемы Стандартной модели: Стандартная модель содержит нерешённые проблемы иерархии флейвора (иерархия масс кварков и лептонов и углы смешивания), сильную CP-проблему, происхождение тёмной материи и механизм инфляции
  2. Вызовы поиска в пространстве параметров: Теории за пределами Стандартной модели обычно включают огромное пространство параметров, где традиционные методы оптимизации неэффективны при поиске дискретных параметров
  3. Сложность аксионной модели: Минимальная аксионная модель объединяет механизм Фроггатта-Нильсена (FN) и механизм Печеи-Куинна (PQ), требуя одновременного удовлетворения ограничений физики флейвора и космологии

Исследовательская мотивация

  • Традиционные методы градиентного спуска не могут эффективно решать задачи оптимизации дискретных параметров
  • Требуется статистический метод для систематического исследования обширного пространства параметров
  • Обучение с подкреплением при ограниченных данных способно самостоятельно обнаруживать решения, что подходит для данного класса задач

Основные вклады

  1. Первое применение обучения с подкреплением к поиску параметров аксионной модели: Разработана стратегия поиска на основе Deep Q-Network (DQN)
  2. Обнаружение большого количества реалистичных решений: Найдено 156 финальных решений, удовлетворяющих ограничениям на перенормированные массы и смешивание
  3. Повышение эффективности: Время вычисления сокращено с минимум 55 дней до 6 дней по сравнению с традиционными методами
  4. Статистический анализ: Проведён систематический статистический анализ распределения решений при различных энергетических масштабах
  5. Феноменологические предсказания: Вычислены связи аксион-фотон, обеспечивающие предсказания для будущих экспериментов, таких как DMRadio-m3

Подробное описание методов

Определение задачи

Входные данные: Вектор распределения U(1)-зарядов Qa={q(Qi),q(ui),q(di),q(ϕ)}Q_a = \{q(Q_i), q(u_i), q(d_i), q(\phi)\} (i=1,2,3) Выходные данные: Распределения U(1)-зарядов, удовлетворяющие ограничениям физики флейвора и космологии Ограничения:

  • Воспроизведение перенормированных масс (Eα<1.75E_\alpha < 1.75)
  • Воспроизведение матрицы смешивания CKM (EijC<0.2E_{ij}^C < 0.2)
  • Положительные целые индексы FN nijn_{ij}

Архитектура модели

Проектирование окружения

  • Пространство состояний: 10-мерный целочисленный вектор, каждая компонента в диапазоне -9, 9
  • Пространство действий: 20 действий (изменения каждого заряда на ±1)
  • Условия ограничений:
    • Ограничение на массу топ-кварка: q(Qˉ3Hcu3)=0q(\bar{Q}_3H^cu_3) = 0
    • Требование положительных целых индексов FN

Архитектура нейронной сети

Входной слой(Z^10) → Скрытый слой 1(R^64) → Скрытый слой 2(R^64) → Скрытый слой 3(R^64) → Выходной слой(R^20)
  • Функция активации: SELU (скрытые слои), Softmax (выходной слой)
  • Оптимизатор: ADAM
  • Функция потерь: потеря Хубера

Проектирование механизма вознаграждения

R(Q,a) = {
    V(Q') - V(Q)  если V(Q') - V(Q) > 0
    -10           если V(Q') - V(Q) ≤ 0
    +100          если Q' — терминальное состояние
}

Внутренняя функция стоимости

V(Q)=minη[Mquark+C]V(Q) = -\min_\eta [M_{quark} + C]

где:

  • Mquark=α=u,dEαM_{quark} = \sum_{\alpha=u,d} E_\alpha (ошибка массы)
  • C=i,jEijCC = \sum_{i,j} E_{ij}^C (ошибка угла смешивания)
  • Eα=log10(mα/mα,RG)E_\alpha = |\log_{10}(|m_\alpha|/|m_{\alpha,RG}|)|

Технические инновации

  1. ε-жадная стратегия: Балансировка исследования и использования, значение ε убывает от 1 до 0.01
  2. Анализ при различных энергетических масштабах: Рассмотрение эффектов перенормировки при четырёх энергетических масштабах M=101417M = 10^{14-17} ГэВ
  3. Двухэтапная оптимизация: Сначала использование RL для поиска дискретных параметров, затем оптимизация непрерывных констант Юкавы методом Монте-Карло
  4. Регулировка статистических пороговых значений: Корректировка порога ошибки массы с 1.0 до 1.75 для адаптации к перенормированным массам

Экспериментальная установка

Набор данных

  • Перенормированные массы: Основаны на массах кварков при различных энергетических масштабах из литературы 29
  • Элементы матрицы CKM: Экспериментально измеренные значения и погрешности
  • Космологические параметры: Ограничения из результатов Planck 2018

Метрики оценки

  • Точность воспроизведения масс: 1.78×102rmass56.21.78 \times 10^{-2} \leq r_{mass} \leq 56.2
  • Точность углов смешивания: 0.63rmixings1.580.63 \leq r_{mixings} \leq 1.58
  • Условия терминального состояния: V(Q)<10.0|V(Q)| < 10.0, Eα<1.75E_\alpha < 1.75, EijC<0.2E_{ij}^C < 0.2

Детали реализации

  • Параметры обучения: 20 агентов, 10510^5 эпизодов, 32 шага/эпизод
  • Скорость обучения: α=2.5×104\alpha = 2.5 \times 10^{-4}
  • Размер пакета: 32
  • Диапазон поиска VEV: 0.01η0.30.01 \leq |\eta| \leq 0.3

Результаты экспериментов

Основные результаты

Статистика обнаружения терминальных состояний

Энергетический масштаб (ГэВ)Кол-во терминальных состоянийПоложительные nijn_{ij}Извлечено при V1=1.0V_1=1.0
101410^{14}71043444
101510^{15}55532352
101610^{16}37423624
101710^{17}54632336
Итого2,1851,316156

Сравнение эффективности

  • Метод RL: 6 дней (одиночный процессор)
  • Оценка традиционного метода: >55 дней
  • Коэффициент ускорения: >9 раз

Результаты статистического анализа

  1. Зависимость от энергетического масштаба: При M=1014,1015M = 10^{14}, 10^{15} ГэВ обнаружено больше высококачественных решений
  2. Распределение числа доменных стенок: Большинство решений находятся при NDW30N_{DW} \approx 30, минимальное значение 20
  3. Распределение внутренней стоимости: Медианные значения при различных энергетических масштабах близки, что указывает на стабильность обучения

Космологические ограничения

  • Тёмная материя: Определяется через связь угла неправильной ориентации θi\theta_i и масштаба PQ faf_a
  • Возмущения кривизны: Накладывают верхние ограничения на масштаб инфляции HinfH_{inf}
  • Ограничения инфляции: Модели с M1015M \gtrsim 10^{15} ГэВ более предпочтительны с космологической точки зрения

Результаты лептонного сектора

  • M=1015M = 10^{15} ГэВ: Найдено 23 модели, удовлетворяющие ограничениям углов смешивания на уровне 3σ
  • M=1016M = 10^{16} ГэВ: Найдено 7 моделей, удовлетворяющих ограничениям
  • Массы нейтрино: Удовлетворяют ограничению mν<85\sum m_\nu < 85 мэВ

Связанные работы

Применение машинного обучения в физике частиц

  • Harvey & Lukas (2021): Применение RL к моделям кварковых масс
  • Предыдущая работа авторов 8: RL в исследовании структуры флейвора лептонного сектора

Исследования аксионных моделей

  • Ema и др. (2017): Предложение модели flaxion
  • Calibbi и др. (2017): Построение минимальной аксионной модели
  • Традиционные методы в основном полагаются на аналитические оценки и ограниченную выборку

Сравнение методов оптимизации

  • Традиционные методы: библиотека scipy, самый быстрый SLSQP требует 0.274 мс/итерация
  • Метод Монте-Карло: низкая эффективность случайного поиска
  • Генетические алгоритмы: систематическое сравнение не проводилось в данной работе

Выводы и обсуждение

Основные выводы

  1. Эффективность метода: RL успешно нашёл большое количество реалистичных распределений U(1)-зарядов, доказав эффективность метода
  2. Вычислительная эффективность: Значительное повышение эффективности поиска по сравнению с традиционными методами
  3. Статистические закономерности: Выявлены характеристики распределения решений при различных энергетических масштабах и космологические предпочтения
  4. Экспериментальные предсказания: Предоставлены конкретные параметрические предсказания для будущих экспериментов, таких как DMRadio-m3

Ограничения

  1. Ограничение на один флейвон: Рассмотрен только один флейвон-поле, что не позволяет обрабатывать нарушение CP
  2. Фокус на кварковый сектор: Основной анализ сосредоточен на кварковом секторе, анализ лептонного сектора относительно ограничен
  3. Зависимость от гиперпараметров: Проектирование вознаграждения и выбор пороговых значений содержат элементы субъективности
  4. Вычислительные ресурсы: Требуются значительные вычислительные ресурсы для обучения

Направления будущих исследований

  1. Расширение на несколько флейвонов: Введение двух флейвон-полей для обработки нарушения CP
  2. Комплексные константы Юкавы: Рассмотрение случаев с комплексными константами связи
  3. Оптимизация архитектуры сети: Разработка более крупных и универсальных нейронных сетей
  4. Другие аксионные модели: Обобщение метода на другие аксионные модели и модели флейвора

Глубокая оценка

Преимущества

  1. Высокая инновационность: Первое систематическое применение обучения с подкреплением к поиску параметров аксионной модели
  2. Высокая практическая ценность: Значительное повышение вычислительной эффективности, делающее возможным поиск в больших пространствах параметров
  3. Богатые результаты: Не только обнаружено большое количество решений, но и проведён детальный статистический анализ
  4. Ориентация на эксперимент: Предоставлены конкретные экспериментальные предсказания, связывающие теорию и эксперимент

Недостатки

  1. Универсальность метода: Проектирование функции вознаграждения сильно зависит от конкретной задачи, что ограничивает обобщаемость
  2. Теоретическая глубина: Отсутствует глубокий теоретический анализ того, почему RL эффективен для данной задачи
  3. Недостаточное сравнение: Сравнение с другими современными методами оптимизации (например, байесовской оптимизацией) недостаточно
  4. Ограниченная верификация: Основная верификация проведена на одной конкретной модели, требуется верификация на большем количестве моделей для подтверждения универсальности

Влияние

  1. Междисциплинарное взаимодействие: Способствует взаимопроникновению машинного обучения и теоретической физики высоких энергий
  2. Методологический вклад: Предоставляет новый подход для решения аналогичных задач оптимизации дискретных параметров
  3. Экспериментальное руководство: Обеспечивает конкретное руководство по пространству параметров для экспериментов по поиску аксионов
  4. Вычислительная физика: Способствует применению методов вычислительной физики в теоретической физике

Применимые сценарии

  1. Модели физики флейвора: Другие модели флейвора с дискретными симметриями
  2. Суперсимметричные модели: Суперсимметричные модели с большим количеством дискретных параметров
  3. Модели с дополнительными измерениями: Теории с дополнительными измерениями, требующие поиска в больших пространствах параметров
  4. Феноменологический анализ: Любые феноменологические исследования, требующие поиска дискретных параметров при наличии ограничений

Библиография

Основные цитируемые работы включают:

  • 1,2 Ema и др., Calibbi и др.: Первоначальное предложение модели flaxion
  • 8 Nishimura и др.: Предыдущее применение RL в физике флейвора авторами
  • 25 Sutton & Barto: Фундаментальная теория обучения с подкреплением
  • 29 Huang & Zhou: Точные вычисления перенормированных масс
  • 9 Сотрудничество DMRadio: Будущие эксперименты по поиску аксионов

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