TreeLevel
Русский
ПлатформаВерсия

Уроки

Четыре маршрута — от первой диаграммы до физики за пределами Стандартной модели. Каждый занимает минут десять; приведённые числа — те, что показывает TreeLevel.

1. e+e−→μ+μ−e^+e^-\to\mu^+\mu^-: первый расчёт, затем пик Z

Две диаграммы e+e- → μ+μ-
  1. Новый документ (⌘NCtrl Nкнопка Создать документ в обозревателе). В палитре уже выбрана модель StandardModel; оставьте направление времени на время →.
  2. Инструмент Вершина: щёлкните дважды по холсту, слева и справа, чтобы поставить две вершины.
  3. Инструмент Линия, частица γ\gamma: протяните от левой вершины к правой. Это пропагатор фотона.
  4. Частица e−/e+e^-/e^+: протяните из пустоты слева внизу к левой вершине — палитра сообщит «входит e−e^-». Повторите слева сверху: стрелка фермионного потока пойдёт назад во времени, и линия будет прочитана как входящий e+e^+. (Порядок, в котором вы проводите две линии, неважен: важно направление стрелки.)
  5. Частица μ−/μ+\mu^-/\mu^+: протяните от правой вершины в пустоту справа вверху (выходящий μ−\mu^-), затем из пустоты справа внизу к вершине (стрелка возвращается к вершине: выходящий μ+\mu^+).
  6. Взгляните на вкладку Анализ: Процесс: e−e+→μ−μ+e^-e^+\to\mu^-\mu^+, диаграмма согласована. Если анализ сообщает о несуществующей вершине, одна из внешних линий читается наоборот: выделите её и отметьте Античастица либо перенесите её конец на другую сторону вершины.
  7. Разверните Правила Фейнмана и амплитуда: вершина −ieγμ-ie\gamma^\mu дважды, пропагатор −igμν/(p1+p2)2-ig_{\mu\nu}/(p_1+p_2)^2, спиноры vˉ(p2)…u(p1)\bar v(p_2)\ldots u(p_1) и uˉ(p3)…v(p4)\bar u(p_3)\ldots v(p_4).
  8. Вкладка Расчёт: s=40 ГэВ\sqrt s = 40\ \text{ГэВ}, Вычислить. σ=62,3 пб\sigma = 62{,}3\ \text{пб}. Проверьте вручную: σ=4πα2/3s\sigma = 4\pi\alpha^2/3s с α=1/127,9\alpha = 1/127{,}9 (значение модели при MZM_Z) и 1 ГэВ−2=3,894×108 пб1\ \text{ГэВ}^{-2} = 3{,}894\times10^8\ \text{пб} даёт 62,3 пб62{,}3\ \text{пб}.
  9. Вернитесь в палитру и выберите Сгенерировать › Добавить недостающие диаграммы: появляется диаграмма с Z, а также обмены h и G0G^0, пропорциональные mem_e и пренебрежимо малые (их убирает пункт Без обмена скалярами в меню генерации). Теперь расчёт идёт по сумме.
  10. График по энергии: от 40 до 200 ГэВ, 300 точек, логарифмическая шкала, Построить. Пик Z достигает 1 979 пб1\,979\ \text{пб} при s=MZ=91,19 ГэВ\sqrt s = M_Z = 91{,}19\ \text{ГэВ}, с шириной ΓZ=2,495 ГэВ\Gamma_Z = 2{,}495\ \text{ГэВ}, объявленной в модели. Кривая проходит через минимум около 70 ГэВ70\ \text{ГэВ} (25 пб25\ \text{пб}), прежде чем подняться к резонансу; интерференция γ\gamma–Z, пропорциональная малым векторным константам gVg_V, в полном сечении остаётся незаметной — она видна в асимметрии вперёд–назад.

Дальше: во вкладке Анализ запустите аналитическое ⟨|ℳ|²⟩ с опущенными нулевыми массами: один только фотонный член возвращает классическое e4(t2+u2)/s2e^4(t^2+u^2)/s^2. Затем нарисуйте рассеяние Бхабхи (e+e−→e+e−e^+e^-\to e^+e^-): две диаграммы, относительный знак и полное сечение, требующее отбора ∣cos⁡θ∣≤0,9|\cos\theta|\le0{,}9.

2. Время жизни мюона

Распад мюона
  1. Новый документ, модель StandardModel. Поставьте две вершины.
  2. Линия W+/W−W^+/W^- между вершинами: виртуальный W.
  3. Частица μ−/μ+\mu^-/\mu^+: протяните из пустоты слева к первой вершине (входящий μ−\mu^-). Частица νμ\nu_\mu: протяните от первой вершины в пустоту (выходящий νμ\nu_\mu).
  4. У второй вершины: выходящий e−e^- (от вершины в пустоту) и выходящий νˉe\bar\nu_e — тяните линию из пустоты к вершине, чтобы стрелка шла назад во времени.
  5. Проверьте во вкладке Анализ: μ−→νμe−νˉe\mu^-\to\nu_\mu e^-\bar\nu_e, согласовано. Классическая ошибка — нарисовать νˉe\bar\nu_e как частицу: запрошенной вершины ve e W~ не существует, и инспектор перечислит соседние, которые есть (ve~ e W, e~ ve W~) — иначе говоря, одну линию надо развернуть.
  6. Вкладка Расчёт: распад трёхчастичный, фазовый объём интегрируется точно (Далиц). Γ=3,007×10−19 ГэВ\Gamma = 3{,}007\times10^{-19}\ \text{ГэВ}, τ=2,19 мкс\tau = 2{,}19\ \text{мкс}, T1/2=1,52 мксT_{1/2} = 1{,}52\ \text{мкс}. Измеренное значение 2,197 мкс2{,}197\ \text{мкс}; разница — радиационные поправки, которых на древесном уровне нет.

Вариант: в модели Hadrons нарисуйте n→pe−νˉen\to pe^-\bar\nu_e через обмен W (τ=946 с\tau = 946\ \text{с} на древесном уровне при gA=1,27g_A = 1{,}27 против измеренных 879 с879\ \text{с}) или π+→μ+νμ\pi^+\to\mu^+\nu_\mu, где всего одна контактная вершина.

3. e+e−→W+W−e^+e^-\to W^+W^-: построить, сложить и увидеть калибровочное сокращение

  1. Новый документ. Поставьте две вершины и проведите внешние линии: входящие e−e^- и e+e^+ слева, выходящие W−W^- и W+W^+ справа (частица W+/W−W^+/W^-: проведённая от вершины в пустоту линия — выходящий W+W^+; из пустоты к вершине, то есть против времени, — выходящий W−W^-). Внутреннюю линию пока не рисуйте.
  2. Сгенерировать › Заменить всеми древесными диаграммами: четыре диаграммы — обмен γ\gamma, Z, h (пренебрежимо мал, ∝me\propto m_e) и νe\nu_e в tt-канале.
  3. Вкладка Расчёт, s=200 ГэВ\sqrt s = 200\ \text{ГэВ}: σ=19,2 пб\sigma = 19{,}2\ \text{пб}, как на LEP (измерено 17 пб17\ \text{пб}, с радиационными поправками).
  4. Удалите диаграмму с нейтрино (выделите её в списке, −) и пересчитайте: 18,1 пб18{,}1\ \text{пб}. Почти ничего не меняется… при 200 ГэВ200\ \text{ГэВ}.
  5. График по энергии, от 170 до 1 000 ГэВ, логарифмическая шкала: без нейтрино σ\sigma растёт (118 пб118\ \text{пб} при 500 ГэВ500\ \text{ГэВ}, 385 пб385\ \text{пб} при 1 ТэВ1\ \text{ТэВ}); со всеми четырьмя диаграммами она падает, как и должна (7,1 пб7{,}1\ \text{пб} при 500 ГэВ500\ \text{ГэВ}, 2,4 пб2{,}4\ \text{пб} при 1 ТэВ1\ \text{ТэВ}). Это и есть сокращение между калибровочными диаграммами — и смысл существования трёхбозонных вершин Стандартной модели.

4. За пределами Стандартной модели: SU(5) и MSSM

SU(5): бозон, превращающий кварки в лептоны

X → u u и X → e+ d̄
  1. Выберите модель SU5. В палитре появляются два бозона, X (заряд 4/34/3) и Y (1/31/3), цветовые триплеты с массой 1015 ГэВ10^{15}\ \text{ГэВ}.
  2. Нарисуйте X→uuX\to uu: одна вершина, входящий X, два выходящих u. Анализ: согласовано, с цветовым множителем εabc\varepsilon_{abc} в правиле — вершина, нарушающая барионное число.
  3. Расчёт: Γ(X→uu)=4,17×1012 ГэВ\Gamma(X\to uu) = 4{,}17\times10^{12}\ \text{ГэВ}. Нарисуйте X→e+dˉX\to e^+\bar d: та же ширина. Относительные вероятности 12/12\tfrac12/\tfrac12 — это результат Джорджи — Глэшоу; для Y вы получите 12 (ud), 14 (e+uˉ), 14 (νˉdˉ)\tfrac12\ (ud),\ \tfrac14\ (e^+\bar u),\ \tfrac14\ (\bar\nu\bar d) и ΓY=ΓX\Gamma_Y = \Gamma_X, как и подобает двум членам одного дублета.
  4. Нарисуйте uu→e+dˉuu\to e^+\bar d и постройте: три диаграммы (X в ss-канале, Y в tt- и uu-каналах). При s=2×1015 ГэВ\sqrt s = 2\times10^{15}\ \text{ГэВ} σ=7×10−26 пб\sigma = 7\times10^{-26}\ \text{пб}; вкладка «Анализ» даёт полное символьное выражение через s,t,MX,MYs, t, M_X, M_Y.

Файл модели (вкладка Теория) показывает, как записываются такие связи: билинейные формы ψTCγμPLχ\psi^TC\gamma^\mu P_L\chi с матрицей зарядового сопряжения CC. См. страницу Теории.

MSSM: импорт UFO

e+e- → нейтралино нейтралино
  1. Скачайте папку UFO MSSM_SLHA2 из базы моделей FeynRules (или любую модель UFO из поставки MadGraph). В палитре нажмите 📁 Импортировать папку UFO и выберите папку.
  2. Через секунду: «UFO импортирован: 50 полей, 822 вершины». Модель появляется в разделе Мои модели, а палитра перечисляет нейтралино, чарджино, скварки, слептоны и глюино с их кодами PDG (1000022 для χ~10\tilde\chi^0_1…).
  3. Нарисуйте e−e+→χ~10χ~10e^-e^+\to\tilde\chi^0_1\tilde\chi^0_1 (две выходящие линии n1; нейтралино — майорановская частица, античастицы нет), затем Сгенерировать: обмен Z в ss-канале, e~L\tilde e_L и e~R\tilde e_R в tt- и uu-каналах.
  4. Расчёт при s=500 ГэВ\sqrt s = 500\ \text{ГэВ}: σ=269 фб\sigma = 269\ \text{фб} в точке SPS1a — в диапазоне исследований для линейного коллайдера. Роли строки и столбца у внешних спиноров меняются от диаграммы к диаграмме: TreeLevel применяет правила фермионного потока Деннера, и результат не зависит от выбранной ориентации.