Extracting Mellin moments of double parton distributions from lattice data

本論文は、ハドロン相関関数の歪み依存性を種々のモデルを通じて解析することにより、ユークリッド格子データから二重パートン分布のメリンモーメントを再構成する際の運動学的歪みの影響を調査する。

要約

陽子を固体の大理石ではなく、パートン（クォークとグルーオン）と呼ばれる小さく目に見えない住民で賑わう混沌とした都市として想像してみてください。物理学者たちは何十年にもわたり、この都市の「人口密度」、つまりどの地域に何人の住民が住み、どれほどの速度で移動しているかをマッピングしてきました。この地図が**パートン分布関数（PDF）**と呼ばれます。

しかし、この都市はあまりにも複雑で、時として二人の住民が全く同じ瞬間に外部世界と相互作用することがあります。これを二重パートン散乱と呼びます。これを理解するためには、**二重パートン分布（DPD）**と呼ばれる、はるかに複雑な新しい地図が必要です。この地図は単に一人の住民がどこにいるかを示すだけでなく、特定の場所にいる特定の二人の住民が、特定の速度で、すべて同時に存在する確率を示します。

問題は、この新しい地図を描くことが極めて困難だということです。顕微鏡で陽子を見るだけでは不十分です。量子力学の法則により、すべてを同時に観測することは不可能だからです。

格子「タイムマシン」

これを解決するため、物理学者たちは格子 QCDと呼ばれるスーパーコンピュータ手法を用います。これは陽子の都市の凍結されたスナップショットを連続して撮影するようなものです。これらのスナップショットが「ユークリッド時間」という、私たちの現実世界の数学的な鏡のような時間の中で存在するという仕組みのため、コンピュータは互いに特定の限られた距離にある住民しか見ることができません。

DPD の全体像を得るためには、物理学者たちはこれらのすべてのスナップショットを、単一の連続した映画に組み合わせる必要があります。数学的には、彼らがイオフェ時間（これを「時間シフト」と呼びましょう）と呼ぶ変数にわたって情報を加算（積分）する必要があります。

ここで問題があります。コンピュータは「時間シフト」の短い期間しかスナップショットを撮影できないのです。10 分間の映像しか持っていないのに、2 時間の映画を再構築しようとしているようなものです。欠落した部分については推測するしかありません。

「スキューネス」のひねり

以前の研究では、著者たちは単純で滑らかな曲線（多項式）を仮定することで、欠落した部分を推測しようとしました。彼らはスキューネス（これを「傾き」と呼びましょう）と呼ばれる変数を導入しました。

• 傾き = 0: これが二重パートン散乱に関心を持つ通常の状態です。
• 傾き = 1: これは奇妙で極端な状態であり、二人の住民が陽子の運動量のほとんどすべてを運び、都市の残りの部分には何も残さない状態です。

著者たちは、以前の「滑らかな曲線」という推測に二つの重大な欠陥があることに気づきました。

1. 端の問題: 彼らの滑らかな曲線は、「傾き」が極端（1 に近い）になったときに、十分に速くゼロに落ちませんでした。物理学は、この極端な状態では、そのような配置を見つける確率が緩やかな傾斜ではなく、崖の縁のようにほぼ瞬時に消滅すべきであることを示唆しています。
2. 滑らかさの問題: 彼らは曲線が至る所で完全に滑らかであると仮定しました。しかし、物理学は、中心（傾き = 0）と端（傾き = 1）の両方で、光源が移動したときに影が急激に変化するのと同様に、曲線に「こぶ」や鋭い点がある可能性を示唆しています。

新しいモデル

この論文において、チームはこれらの「崖の縁」や「こぶ」を尊重するように設計された、欠落した映像を推測する四つの新しい方法を試みました。

1. べき乗則: 端で急激に減少する曲線。
2. 積分モデル: 高エネルギー衝突における粒子の分裂の仕方に基づいた形状。
3. 余弦モデル: 鋭い縁または滑らかな縁を持つように調整できる波状の形状。
4. 多項式（旧方式）: 比較のために保持された、以前使用していた滑らかな曲線。

結果：欠落したピースを持つパズル

チームは、どのモデルが最もよく適合するかを確認するために、コンピュータデータをこれらの新しいモデルに投入しました。

• 良い知らせ: すべての新しいモデルが、利用可能なコンピュータデータに非常に良く適合しました。それらはすべて物語の「中央部」（中程度の傾きにおける振る舞い）について合意していました。
• 悪い知らせ: これらのモデルを用いて地図の最も重要な部分、すなわち傾き = 0（実際の二重パートン散乱事象）を再構築しようとしたとき、結果は極めて不確実でした。
  • コンピュータデータは「時間シフト」の極端な端（欠落した映像がある場所）で非常に「ノイズ」が多く（ぼやけて）なるため、異なるモデルが中心部に対して非常に異なる答えを出しました。
  • 一部のモデルは、2 という値を予測しました（これは「数則」と呼ばれる基本的な法則が示す値です）。
  • 他のモデルは、はるかに外れた値を予測するか、値自体の数百倍もの大きな誤差範囲（不確実性の範囲）を持っていました。

結論

著者たちは、現在のコンピュータデータのみを用いては、最も重要な点（傾き = 0）における二重パートン分布を完全に再構築することはできないと結論付けました。

これは、隅のピースと中央のピースはすべて持っているが、それらをつなぐピースが欠けているパズルのようなものです。パズルの形状を推測することはできますが、より多くの情報がなければ、中央のピースがどのように正確に組み合わさっているかは確信できません。

これを修正するためには、「時間シフト」の極端な端で「ノイズ」の少ない、より良いコンピュータデータが必要だと言っています。それまでの間、データに任せるのではなく、答えを正しく強制するために、数則のような追加の理論的規則に頼らざるを得ません。

要約すると: 彼らは複雑な量子地図の形状を推測するためのより良いツールを構築しましたが、現在の陽子の「写真」は、最も重要な詳細を明確に見るには鮮明さが不足していることがわかりました。作業を完了するには、より鮮明な写真が必要です。

技術概要

技術的サマリー：格子データからの二重部分子分布のメリンモーメントの抽出

問題定義
二重部分子分布（DPDs）は、陽子内の 2 つの部分子間の量子相関を記述し、LHC などの衝突器における二重部分子散乱（DPS）を含む多部分子相互作用の理解に不可欠である。単一部分子分布関数（PDFs）はよく制約されているのに対し、DPDs はその複雑さと実験的な抽出の困難さにより、未だ十分に理解されていない。格子 QCD は DPDs を制約するための非摂動的アプローチを提供する。具体的には、DPDs のメリンモーメントは、異なる空間位置にある 2 つの局所電流のハドロン行列要素と関連付けられることができる。しかし、これらのモーメントをユークリッド格子データから再構成するには、$-\infty$ から $+\infty$ まで変化する変数 $\omega$（イオフェ時間）にわたって積分を行う必要がある。格子シミュレーションは $\omega$ の有限範囲に制限されている。この課題に対処するため、著者らは以前、$\omega$ とフーリエ共役な「スキューネス」パラメータ $\zeta$ を導入した。物理的な DPDs は $\zeta = 0$ に対応する。課題は、制御されていないバイアスを導入することなく必要な逆フーリエ変換を実行するために、メリンモーメントの $\zeta$ に対する関数依存性を決定することにある。

手法
著者らは、既存の CLS H102 集合体（$N_f=2+1$）からの格子データを用いて、DPDs の最低次メリンモーメント $I_{q_1q_2}(\zeta, y^2)$ の抽出を再検討する。手法は 3 つの主要な段階で進行する。

1. 理論的制約: 本論文は、以前の多項式 Ansatz によって完全には満たされなかった、メリンモーメントの $\zeta$ 依存性に関する 2 つの重要な物理的性質を確立する。

   • モーメントは $|\zeta| \to 1$ として急速に消滅しなければならない。
   • モーメントは、$\zeta = 0$ および $\zeta = 1$ において非解析的な挙動を示す可能性がある。これは、運動量分数 $x=0$ および $x=1$ における PDFs の挙動に類似している。
     これらの性質は、スキュード DPDs の支持領域と、短距離極限における摂動的分裂機構から導き出される。

2. 格子データ解析: 本研究は、陽子状態における 2 つのベクトル電流の 4 点相関関数を利用する。不変関数 $A_{q_1q_2}(\omega, y^2)$ が格子データから抽出される。解析は、等価なフレーバーの組み合わせと比較して優れた統計的精度を持つため、フレーバー組み合わせ $ud$に焦点を当てる。データは、離散化効果と有限サイズ効果を最小化するために、距離範囲$4a \le y \le 14a$ に制限される。

3. モデリングとフィッティング: 著者らは、$\zeta$ と $y$ の依存性が $I(\zeta, y^2) = J(\zeta)K(y^2)$ と因子分解されると仮定する。彼らは、格子データの正規化された $\omega$ 依存性 $\hat{A}(\omega)$ を、$J(\zeta)$ に対するいくつかの新しい Ansatz のフーリエ変換にフィットさせる。

   • 多項式モデル: 以前の手法（$\zeta^2$ の多項式）。
   • べき乗則モデル: $(1-|\zeta|)^a$ に比例する。
   • 積分モデル: 摂動的分裂に触発された積分表現に基づく。
   • 余弦モデル: $\cos^b(\frac{\pi}{2}|\zeta|^a)$ に比例する。

   フィッティングは、因子分解仮説を検証するために、グローバルに行われるとともに、$y$ の狭いビン内でも行われる。さらに、数総和則（$S_{ud} = 2$）を満たすようにパラメータを固定した制約付きフィッティングも実施される。

主要な結果

• モデル比較: 多項式モデルは、$\omega$ 空間における格子データへの良好なフィットを提供するが、物理的な境界条件（$\zeta=1$ での非消滅と $\zeta=0$ での滑らかさ）を満たすことに失敗する。新しいモデル（べき乗則、積分、余弦）はこれらの条件を満たし、データに対して同様に良好なフィットを提供する。
• 再構成の不確実性: 2 つの自由パラメータでフィットする場合、新しい柔軟なモデルは、物理的な DPD 極限である $\zeta=0$ におけるメリンモーメントに対して極めて大きな不確実性を生じさせる。例えば、積分モデルは、理論的な期待値が 2 であるのに対し、総和則量に対して $4.4 \pm 938.4$ という値を与える。これは、追加の理論的入力なしに、特に大きな $|\omega|$ における格子データが、$\zeta=0$ 近傍の関数の挙動を制約するには不十分であることを示している。
• 制約付きフィッティング: 新しいモデルの 1 つのパラメータを固定して数総和則（$S_{ud}=2$）を強制すると、データと整合する安定した滑らかなメリンモーメントの再構成が得られる。
• 平均スキューネス: すべてのフィットは、一貫して小さな平均二乗スキューネス $\langle \zeta^2 \rangle \approx 0.07$ を生じさせ、小さなスキューネスを持つ部分子配置が支配的であることを示唆している。
• 因子分解: 本研究は、$\zeta$ と $y$ の依存性の因子分解が、現在のデータの統計的精度の範囲内で成り立つことを確認する。

意義と主張
本論文は、格子 QCD が DPDs に関連する不変関数に対する高精度データを提供し得る一方で、ゼロ・スキューネス（$\zeta=0$）における物理的メリンモーメントの再構成は、現在、アクセス可能なイオフェ時間（$\omega$）の有限範囲によって制限されていると主張する。著者らは、特定の関数形式（多項式など）を仮定することが物理的に矛盾した結果をもたらす可能性があるのに対し、物理的に動機付けられたモデル（べき乗則、積分、余弦）は必要であるが、制約を課さない場合、$\zeta=0$ の極限において大きな統計的不確実性を導入することを示している。

主要な結論は控えめである：現在利用可能な格子データでは、追加の理論的入力（例えば、数総和則の強制など）なしに、$\zeta=0$ におけるメリンモーメントを再構成することはできない。 著者らは、この状況を改善するには、アクセス可能な $\omega$ の範囲を拡大するために、はるかに高い精度で大きな陽子運動量および大きな距離における格子シミュレーションが必要であると強調している。しかし、非ゼロ $\zeta$ に対する結果は依然として価値があり、部分子の結合分布に関する情報を提供し、小さなスキューネスを持つ配置の方が確率的に高いことを確認している。本研究はまた、$\zeta$ 依存性および理論的入力の必要性に関する教訓は、準分布または疑似分布を用いた将来の DPDs の格子研究に関連することを示唆している。