NNLOJET: a parton-level event generator for jet cross sections at NNLO QCD accuracy

本論文は、電子・陽電子、レプトン・ハドロン、およびハドロン・ハドロン衝突における NNLO QCD 精度でのジェット断面積および関連する観測量を計算するためにアンテナ減算法を実装したオープンソースの NNLOJET パートンレベル事象生成器を導入する。

要約

以下は、この論文を日常的な言葉と比喩を用いて翻訳・解説したものです。

全体像：より優れた水晶玉の構築

あなたが物理学者で、ほぼ光速で衝突する 2 つの微小な粒子が何を起こすかを予測しようとしていると想像してください。あなたには、これらの粒子が「どのように振る舞うべきか」を示す一連の規則（標準模型）があります。しかし、自然は厄介です。粒子が衝突すると、単に跳ね返るだけでなく、大砲から飛び散るコンフェッティのように、無数の新しい粒子が混沌としたシャワーとなって噴き出します。

真に正確な予測を行うためには、単に主要な爆発を推測するだけでは不十分です。爆発に「伴って」起こる、ごく小さくほとんど目に見えないさざ波や二次的な噴出まで計算する必要があります。粒子物理学の世界では、これらの計算をNNLO（次々項）と呼びます。これらは非常に複雑で、ハリケーンの中の砂粒一つ一つの正確な軌道を予測しようとするようなものです。

この論文は、これらの超精密計算を「ジェット（粒子の噴流）」に対して行うために設計された、新しいオープンソースのコンピュータプログラムNNLOJETを紹介しています。これは、古いモデルが見逃していた詳細まで捉えることができる、高解像度の水晶玉を科学者に与えるようなものです。

核心的な問題：「無限」の混乱

物理学者がこれらの衝突を計算しようとすると、「赤外特異点」と呼ばれる数学的な悪夢に直面します。

• 比喩: 部屋にいる人の数を数えようとしているが、誰かがささやくたびに部屋が無限に騒がしくなる状況を想像してください。ノイズを合計しようとすると、数は無限大になってしまいます。
• 現実: 粒子物理学では、粒子が非常にゆっくりと動いたり、非常に近接したりすると、数学が無限大に発散します。現実世界に無限のエネルギーが存在しない以上、これらの無限大は計算手法における単なる誤りです。

これを修正するために、著者らはアンテナ減算法と呼ばれる手法を使用します。

• 比喩: ノイズキャンセリングヘッドフォンだと考えてください。ヘッドフォンは煩わしい背景ノイズ（数学的な無限大）を聞き取り、それを完全に打ち消す「逆位相」の音波を生成します。
• 仕組み: プログラムは、厄介で無限大になる部分を個別に計算し、それからそれらを差し引きます。そうすることで、クリーンで有限の、現実世界に即した答えのみが残ります。

NNLOJET とは何か

NNLOJETは、これらの計算を実行する「エンジン」です。この論文以前、これらの計算の多くは、ある一人のシェフのキッチンに閉じ込められた秘密のレシピのようなものでした。それを使いたければ、シェフに懇願するか、ゼロから自分のキッチンを作らなければなりませんでした。

NNLOJET はオープンソースであることでゲームを変えます。

• 比喩: 全レシピ本と調理器具を一般公開するようなものです。コンピュータを持っている人なら誰でもダウンロードして、料理を作り、材料を確認できます。
• 「ジェット」の部分: このプログラムは「ジェット」に特化しています。粒子物理学において、物質の構成要素であるクォークやグルーオンが放出されると、単独で留まることはありません。それは瞬時に粒子の噴流へと変わります。この噴流を「ジェット」と呼びます。NNLOJET は、これらの噴流がどのくらい大きく、どのくらい速く、どの方向に進むかを正確に予測します。

どのように使うのか（ランカード）

NNLOJET を使うためにコーディングの魔術師である必要はありません。この論文では、ランカードを使ってプログラムを制御すると説明しています。

• 比喩: ランカードをフライトプランやレシピカードだと考えてください。飛行機のエンジンがどのように動くかを知らなくても、パイロットに行き先を指示するだけでよいのです。
• あなたが書くこと: プログラムに以下を指示します。
  • どこで: 陽子（大型ハドロン衝突型加速器 LHC のようなもの）を衝突させるのか、それとも電子を？
  • 何を: Z ボソン、ヒッグス粒子、それとも単なるジェット噴流を探しているのか？
  • 規則: ジェットはどのくらい大きく、どのくらい速くあるべきか？
  • 出力: 最終的にどのようなチャートを見たいか？

この論文は、衝突の種類から特定の数学設定に至るまで、この「フライトプラン」の書き方について詳細なマニュアル（第 5 節と第 6 節）を提供しています。

ワークフロー：組立ライン

これらの衝突を計算するのは非常に重く、1 台のコンピュータだけでは実行できません。何年もかかるでしょう。

• 比喩: 巨大な壁画を描く必要があると想像してください。1 人なら一生かかります。代わりに、100 人の画家のチームを雇います。
• ツール: この論文では、現場監督として機能するワークフロースクリプト（nnlojet-run）について説明しています。
  1. ウォーミングアップ: 現場監督は数人の画家を送り出し、壁をテストし、最も良い描画方法を見つける（これは「VEGAS グリッド」の適応です）。
  2. 本番: 計画が固まると、現場監督は数百人の画家を同時に実際の作業に送り出します。
  3. 完了: 全員が終了すると、現場監督はすべてのピースを集め、誤り（外れ値）をチェックし、それらを貼り合わせて 1 つの完璧な絵に仕上げます。

何の計算ができるのか

この論文は、現在のレシピブック（バージョン 1.0.0）で利用可能な特定の「レシピ」をリストアップしています。

• 電子・陽電子衝突: 過去の LEP コライダーのようなもの。
• 陽子・陽子衝突: 現在の大型ハドロン衝突型加速器（LHC）のようなもの。
• 特定の事象: 以下の生成を予測できます。
  • Z ボソン（光子の重い親戚）。
  • W ボソン（弱い力を運ぶ粒子）。
  • ヒッグス粒子（他の粒子に質量を与える粒子）。
  • 光子（光）とジェット。
  • これらすべての組み合わせ（例：ヒッグス粒子＋ジェット）。

なぜこれが重要なのか

この論文は、このコードをオープンで使いやすいものにすることで、科学者が理論的な予測を実験データと直接比較できるようになると結論付けています。

• 比喩: モデルカーを構築している場合、それが実際に走行するかどうかを知る必要があります。NNLOJET は、可能な限り最も正確な「風洞実験」を提供します。実データが NNLOJET の予測と一致すれば、宇宙に対する私たちの理解は正しいことになります。一致しない場合、それは何か新しいものを見つけたことを意味するかもしれません！

要約すると: この論文は、物理学者が粒子衝突の厄介な詳細を極めて高い精度で計算できるよう支援する、強力で無料かつ使いやすいツールを公開するものです。これにより、宇宙の地図が可能な限り正確であることが保証されます。

技術概要

NNLOJET の技術的概要：NNLO QCD 精度におけるジェット断面積のためのパートンレベル事象生成器

問題提起
標準模型（SM）の精密検証は、ハドロン終状態（ジェット）の実験的測定値と理論的予測との比較に大きく依存している。これらの予測には、大きな強い結合定数（$\alpha_s$）を扱い、軟・共線放射に起因する複雑な赤外（IR）特異性構造を処理するために、次々々次世代（NNLO）量子色力学（QCD）補正の計算が必要となる。次世代（NLO）計算はほぼ自動化されているが、二重仮想（VV）、混合実 - 仮想（RV）、および二重実（RR）の寄与間の IR 発散を相殺する必要があるため、NNLO 計算は依然として困難である。既存の NNLO コードは多くの場合、特定のプロセスに限定されており、オープンソースとして広く利用可能ではないか、あるいは各観測量ごとに大幅な手動調整を必要とする。電子 - 陽電子、レプトン - ハドロン、およびハドロン - ハドロンといった様々な衝突環境にわたって NNLO 精度でジェット断面積を計算できる、柔軟でオープンソースのパートンレベル事象生成器が必要とされている。

手法：アンテナ減算法
NNLOJET コードは、IR 特異性を処理するためにアンテナ減算法を実装している。この手法は、IR 極限における QCD 実放射振幅の普遍的な因子分解に基づいて減算項を導出する。

• 中核メカニズム: この手法は、未解決の実放射構成を、2 つの硬い放射体パートンのペアとして特定する。それらの間の放射は「アンテナ関数」によって記述される。
• 減算項: 未積分の減算項は、アンテナ関数とより低いパートン多重性の縮小行列要素の積から構成される。これらの項はアンテナ位相空間上で積分され、仮想補正からの$\epsilon$極を明示的に相殺する積分済みの減算項（$d\hat{\sigma}_T$および$d\hat{\sigma}_U$）を生成する。
• 実装: コードは、アンテナ関数の運動学を縮小行列要素から因子分解する位相空間マッピングを利用する。実装は、終 - 終、始 - 終、および始 - 始の放射体構成を網羅している。
• 計算フレームワーク: NNLOJET は、主に現代の Fortran で書かれ、C++ および Python モジュールを備えたパートンレベル事象生成器である。位相空間積分には適応型モンテカルロルーチン VEGAS を使用する。コードには、HPLOG、TDHPL、CHAPLIN を通じて特殊関数（調和多対数関数）を評価するためのライブラリが含まれており、パートン分布関数（PDF）には LHAPDF とインターフェースする。

主要な貢献と特徴
本論文は、GNU General Public License v3.0 の下で配布される NNLOJET コードのオープンソース版（v1.0.0）を発表する。主な技術的貢献は以下の通りである：

1. 広範なプロセス対応: このコードは、以下のジェット生成に対する NNLO 計算をサポートする：

   • 電子 - 陽電子衝突: $e^+e^- \to 2j$ および $3j$。
   • レプトン - ハドロン衝突: 中性カレント（NC）および荷電カレント（CC）DIS プロセス（$ep \to \ell + 2j$、$ep \to \nu + 2j$など）。
   • ハドロン - ハドロン衝突: 包括的ジェット生成（$pp \to j$）、ダイジェット生成（$pp \to 2j$）、ベクトルボソン＋ジェット（$V+1j$）、光子＋ジェット（$\gamma+1j$）、ダイ光子（$\gamma\gamma$）、およびヒッグスボソン＋ジェット（$H+1j$）生成。
   • 派生次数: このフレームワークは、1 つ追加のジェットを持つプロセスに対する NLO 補正と、2 つ追加のジェットを持つプロセスに対する LO を暗黙的に提供する（例：2 ジェット要求を伴う $ZJ$プロセスを通じて$Z+2j$ を NLO で計算するなど）。

2. 柔軟なユーザーインターフェース（ランカード）: コードは、プレーン ASCII の「ランカード」ファイルを通じて制御され、ユーザーは以下を定義できる：

   • プロセスと衝突器: 衝突器タイプ（pp、ep、$e^+e^-$）、ビームエネルギー、およびジェットアルゴリズム（anti-$k_T$、Cambridge-Aachen、$k_T$、Durham、Jade）の選択。
   • 観測量と選択: 忠実度カット、オブジェクトレベルのセレクター（ジェット $p_T$、ラピディティ）、およびイベントレベルのセレクターの定義。
   • ヒストグラム: 柔軟なビニング（線形、対数、カスタムエッジ）およびヒストグラムセレクターを介した多次元微分分布。
   • スケール: 不確実性推定のための、複数の繰り込みスケールおよび因子化スケールの組み合わせの同時評価。

3. 自動化されたワークフロー: 配布物には、ジョブ提出、リソース割り当て、結果の結合を管理するための Python ベースのワークフロースクリプト（nnlojet-run）が含まれている。このワークフローは以下を処理する：

   • ウォームアップおよび生成フェーズ: 適応型グリッド生成に続き、生成実行を行う。
   • リソース管理: 目標相対精度に達するための計算リソースの動的調整。
   • 結果結合: 外れ値の影響を軽減しつつ、独立したジョブからの部分的な結果を結合するための堅牢なアルゴリズム（トリミング、k-scan、重み付き結合を含む）。

4. 検証: 実装には厳格な検証チェックが含まれている：

   • ポイントテスト: 特異極限近傍で、行列要素と減算項の比率が 1 に近づくことを確認する。
   • 技術的カットチェック: 結果が技術的赤外カットオフに依存しないことを確認する。
   • 極テスト: RV および VV 寄与における$\epsilon$極の明示的な相殺。

結果と検証
本論文は、各種プロセスの実装の詳細を記述し、検証の詳細については特定の出版物を参照している。具体的な例として、著者らは LHC（$\sqrt{s} = 7$ TeV）における$Z + 1j$生成の計算を提示する：

• 設定: anti-$k_T$ジェットアルゴリズム（$R=0.4$）、NNPDF3.1 NNLO PDF、およびレプトンとジェットに対する特定のカット。
• 結果: 総 NNLO 忠実度断面積は 68.5(2) pb と計算された。
• 微分分布: 論文は、先頭ジェットの横運動量分布を示し、LO から NLO、そして NNLO への収束を提示する。不確実性バンドは 7 点スケール変換を表す。

意義と主張
著者らは、NNLOJET が、従来は専有であったり特定のプロセスクラスに限定されていたりしたコードが残したギャップを埋める、ジェット観測量に対する NNLO QCD 計算のための普遍的に適用可能なオープンソースフレームワークを提供すると主張している。

• アクセシビリティ: GPLv3.0 の下でコードを配布し、ユーザーフレンドリーなランカード構文とワークフロースクリプトを提供することで、著者らは高精度の NNLO 予測を、より広範な現象論家および実験家のコミュニティが利用できるようにすることを目的としている。
• 柔軟性: コードを再コンパイルすることなく任意の忠実度カットや観測量を定義できる能力により、実験データとの直接比較が可能となる。
• 拡張性: モジュラー設計により、将来のリリースで新しいプロセスや機能の追加が可能である。

本論文は、NNLOJET が精密物理学研究、標準模型のストレステスト、および SM パラメータとパートン分布関数の決定の精緻化のための成熟したツールであると結論付けている。コードは HEPForge で利用可能であり、著者らはユーザーからのフィードバックとバグ報告を歓迎している。