2025-11-24T01:49:17.881693

Mass estimates of the SU(2) $0^{++}$ glueball from spectral methods

Dudal, Oliveira, Roelfs

The estimation of the KÃ¤llÃ©n-Lehmann spectral density from gauge invariant lattice QCD two point correlation functions is proposed, and explored via an inversion strategy based on Tikhonov regularisation. We test the method on a mesonic toy model, showing that our methodology is competitive with the traditional Maximum Entropy Method. As proof of concept the SU(2) glueball spectrum for the quantum numbers $J^{PC}=0^{++}$ is investigated, for various values of the lattice spacing, using the published data of arXiv:1910.07756. Our estimates for the ground state mass are in good agreement with the traditional approach, which is based on the large time exponential behaviour of the correlation functions. Furthermore, the spectral density also contains hints of excites states in the spectrum. Spectroscopic analysis of glueball two-point functions therefore provides a straightforward and insightful alternative to the traditional method based on the large time exponential behaviour of the correlation functions.

academic

SU(2) $0^{++}$ グルーボール質量のスペクトル法による推定

基本情報

論文ID: 2112.06785
タイトル: Mass estimates of the SU(2) $0^{++}$ glueball from spectral methods
著者: David Dudal (KU Leuven & Ghent University)、Orlando Oliveira (University of Coimbra)、Martin Roelfs (KU Leuven)
分類: hep-lat(格子場の理論)、hep-ph(現象論)
発表時期/会議: 第38回国際格子場理論シンポジウム、LATTICE2021(2021年7月26-30日)
論文リンク: https://arxiv.org/abs/2112.06785

要旨

本論文は、ゲージ不変な格子QCD二点相関関数からKällén-Lehmannスペクトル密度を推定する方法を提案し、Tikhonov正則化に基づく逆変換戦略を通じて探索している。著者らはメソン玩具モデルでこの方法をテストし、従来の最大エントロピー法と競争力のあることを示した。概念実証として、量子数 $J^{PC}=0^{++}$ を持つSU(2)グルーボールスペクトルを異なる格子間隔のデータを用いて研究した。基底状態質量推定値は相関関数の大時間指数減衰に基づく従来の方法と一致し、スペクトル密度はさらに励起状態の兆候を含んでいる。

研究背景と動機

解決すべき問題: 従来の格子QCDにおける粒子質量計算は二点相関関数の大時間減衰挙動の分析に依存しているが、本論文はKällén-Lehmannスペクトル表現の直接逆変換を通じて粒子質量情報を取得することを提案している。
問題の重要性:
- グルーボールは純ゲージ理論における束縛状態であり、QCDの非摂動的性質を理解するための鍵である
- スカラーグルーボール( $0^{++}$ )は最軽量のグルーボール状態であり、理論的および実験的に重要な意義を持つ
- スメアリング技術に依存しない代替分析方法を提供する
既存方法の限界:
- 従来の方法は信号対雑音比を改善するためにスメアリング技術を使用する必要がある
- 大時間挙動の指数フィッティングに依存しており、励起状態情報を失う可能性がある
- 最大エントロピー法は広く使用されているが、場合によっては十分に安定していない可能性がある
研究動機: 運動量空間伝播関数からスペクトル情報を直接抽出する方法を開発し、従来の方法の技術的複雑性を回避しながら励起状態を識別できる手法を確立する。

核心的貢献

Tikhonov正則化に基づくスペクトル密度逆変換方法を提案し、非負最小二乗法(NNLS)と組み合わせて物理的に妥当な正のスペクトル密度を確保した
メソン玩具モデルで方法の有効性を検証し、最大エントロピー法と同等の性能を示した
SU(2) $0^{++}$ グルーボール質量推定に成功し、複数の格子間隔で従来の方法と一致する基底状態質量を得た
励起状態の存在を識別し、グルーボールスペクトル学に追加情報を提供した
スメアリング技術を必要としない分析経路を提供し、格子QCD分析プロセスを簡素化した

方法の詳細

タスク定義

格子QCD二点相関関数 $G(p)$ からKällén-Lehmannスペクトル密度 $\rho(\omega)$ を抽出し、その後グルーボールの質量スペクトルを決定する。入力は運動量空間の相関関数データであり、出力は周波数領域のスペクトル密度関数である。

モデルアーキテクチャ

1. Källén-Lehmannスペクトル表現

相関関数とスペクトル密度の関係は以下の通りである: $G(p) = \int_0^{\infty} \frac{2\omega\rho(\omega) d\omega}{\omega^2 + p^2}$

2. 減算スキームの処理

発散問題に対処するため、減算を導入する: $G(p^2) = \sum_{k=0}^{n-1} a_k (p^2 - \bar{p}^2)^k + (-p^2 + \bar{p}^2)^n \int_0^{\infty} \frac{2\omega\tilde{\rho}(\omega) d\omega}{\omega^2 + p^2}$

3. Schwinger関数法

直接的な数値減算を回避するためにFourier変換を使用する: $C(\tau) = \mathcal{F}\{G(p)\}(\tau) = \int_{-\infty}^{\infty} G(p^2)|_{\vec{p}=0, p_4\neq 0} e^{-ip_4\tau} dp_4$

$\tau \neq 0$ に対して、 $C(\tau) = \mathcal{L}\{\rho(\omega)\}(|\tau|)$ 、すなわちLaplace変換関係が成り立つ。

4. Tikhonov正則化逆変換

Laplace逆変換問題を正則化最小二乗問題に変換する: $J_\alpha = \|\mathbf{K}\boldsymbol{\rho} - \mathbf{C}\|_2^2 + \alpha^2\|\boldsymbol{\rho} - \boldsymbol{\rho}^*\|_2^2$

ここで $\mathbf{K}_{k\ell} = e^{-\omega_\ell \tau_k}\Delta\omega$ 、 $\boldsymbol{\rho}^*$ は事前推定値である。

5. 非負制約

NNLS求解器を通じて $\rho(\omega) \geq 0$ を確保する: $\mathbf{A} = \begin{pmatrix} \mathbf{K} \\ \alpha\mathbf{1} \end{pmatrix}, \quad \mathbf{b} = \begin{pmatrix} \mathbf{C} \\ \alpha\boldsymbol{\rho}^* \end{pmatrix}$

$J_\alpha = \|\mathbf{A}\boldsymbol{\rho} - \mathbf{b}\|_2^2$ を最小化し、制約 $\boldsymbol{\rho} \geq 0$ を満たす。

技術的革新点

数値減算の回避: Schwinger関数法を通じて直接的な数値減算操作を巧妙に回避した
物理的制約: 正性制約を強制実装し、スペクトル密度の物理的妥当性を確保した
パラメータ選択: Morozov差異原理を採用して正則化パラメータ $\alpha$ を自動選択した
対数サンプリング: 周波数領域で対数均一サンプリングを使用し、広い周波数範囲 $[10^{-5}, 10^5]$ GeVをカバーした

実験設定

データセット

Yamanaka等Phys. Rev. D 102, 054507 (2020)が発表したSU(2)純Yang-Mills理論データを使用:

$\beta$	格子体積	配置数	$N$	$a\sqrt{\sigma}$	$a[\Lambda^{-1}]$
2.1	$10^3 \times 12$	1,000,000	7	0.608(16)	0.356(27)
2.2	$12^4$	9,999,990	7	0.467(10)	0.273(20)
2.3	$14^3 \times 16$	4,100,000	9	0.3687(22)	0.216(15)
2.4	$16^3 \times 24$	2,030,000	13	0.2660(21)	0.156(11)
2.5	$20^3 \times 24$	520,000	13	0.1881(28)	0.110(8)