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.
논문 ID : 2112.06785제목 : Mass estimates of the SU(2) 0 + + 0^{++} 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 P C = 0 + + J^{PC}=0^{++} J PC = 0 ++ 인 SU(2) 글루볼 스펙트럼을 다양한 격자 간격 데이터를 사용하여 연구했다. 기저 상태 질량 추정은 상관함수의 큰 시간 지수 거동에 기반한 전통적 방법과 일치하며, 스펙트럼 밀도는 여기 상태의 흔적도 포함한다.
해결하려는 문제 : 전통적인 격자 QCD에서 입자 질량 계산은 이점 상관함수의 큰 시간 감쇠 거동 분석에 의존하는데, 본 논문은 Källén-Lehmann 스펙트럼 표현의 직접 역변환을 통해 입자 질량 정보를 획득하는 방법을 제안한다.문제의 중요성 :글루볼은 순수 규범 이론의 결합 상태로서 QCD의 비섭동적 성질을 이해하는 핵심 스칼라 글루볼(0 + + 0^{++} 0 ++ )은 가장 가벼운 글루볼 상태로서 이론적, 실험적 중요성을 가짐 도포 기법에 의존하지 않는 대체 분석 방법 제공 기존 방법의 한계 :전통적 방법은 신호 대 잡음비 개선을 위해 도포 기법 사용 필요 큰 시간 거동의 지수 적합에 의존하여 여기 상태 정보 손실 가능 최대 엔트로피 방법은 광범위하게 사용되지만 특정 경우 안정성 부족 가능 연구 동기 : 동량 공간 전파함수로부터 스펙트럼 정보를 직접 추출하는 방법 개발로 전통적 방법의 기술적 복잡성을 회피하면서 여기 상태를 식별할 수 있는 방법 개발Tikhonov 정규화 기반의 스펙트럼 밀도 역변환 방법 제안 , 음이 아닌 최소제곱(NNLS)과 결합하여 물리적으로 합리적인 양의 스펙트럼 밀도 보장중간자 장난감 모형에서 방법의 유효성 검증 , 최대 엔트로피 방법과 비교 가능한 성능 입증SU(2) 0 + + 0^{++} 0 ++ 글루볼 질량 추정에 성공적 적용 , 다양한 격자 간격에서 전통적 방법과 일치하는 기저 상태 질량 획득여기 상태의 존재 식별 , 글루볼 스펙트럼학에 추가 정보 제공도포 기법이 필요 없는 분석 경로 제공 , 격자 QCD 분석 절차 단순화격자 QCD 이점 상관함수 G ( p ) G(p) G ( p ) 로부터 Källén-Lehmann 스펙트럼 밀도 ρ ( ω ) \rho(\omega) ρ ( ω ) 를 추출하여 글루볼의 질량 스펙트럼을 결정. 입력은 동량 공간의 상관함수 데이터이고 출력은 주파수 영역의 스펙트럼 밀도 함수.
상관함수와 스펙트럼 밀도의 관계:
G ( p ) = ∫ 0 ∞ 2 ω ρ ( ω ) d ω ω 2 + p 2 G(p) = \int_0^{\infty} \frac{2\omega\rho(\omega) d\omega}{\omega^2 + p^2} G ( p ) = ∫ 0 ∞ ω 2 + p 2 2 ω ρ ( ω ) d ω
발산 문제 처리를 위한 감산 도입:
G ( p 2 ) = ∑ k = 0 n − 1 a k ( p 2 − p ˉ 2 ) k + ( − p 2 + p ˉ 2 ) n ∫ 0 ∞ 2 ω ρ ~ ( ω ) d ω ω 2 + p 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} G ( p 2 ) = ∑ k = 0 n − 1 a k ( p 2 − p ˉ 2 ) k + ( − p 2 + p ˉ 2 ) n ∫ 0 ∞ ω 2 + p 2 2 ω ρ ~ ( ω ) d ω
직접 수치 감산을 회피하기 위한 Fourier 변환:
C ( τ ) = F { G ( p ) } ( τ ) = ∫ − ∞ ∞ G ( p 2 ) ∣ p ⃗ = 0 , p 4 ≠ 0 e − i p 4 τ d p 4 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 C ( τ ) = F { G ( p )} ( τ ) = ∫ − ∞ ∞ G ( p 2 ) ∣ p = 0 , p 4 = 0 e − i p 4 τ d p 4
τ ≠ 0 \tau \neq 0 τ = 0 에 대해, C ( τ ) = L { ρ ( ω ) } ( ∣ τ ∣ ) C(\tau) = \mathcal{L}\{\rho(\omega)\}(|\tau|) C ( τ ) = L { ρ ( ω )} ( ∣ τ ∣ ) , 즉 Laplace 변환 관계.
Laplace 역변환 문제를 정규화된 최소제곱 문제로 변환:
J α = ∥ K ρ − C ∥ 2 2 + α 2 ∥ ρ − ρ ∗ ∥ 2 2 J_\alpha = \|\mathbf{K}\boldsymbol{\rho} - \mathbf{C}\|_2^2 + \alpha^2\|\boldsymbol{\rho} - \boldsymbol{\rho}^*\|_2^2 J α = ∥ K ρ − C ∥ 2 2 + α 2 ∥ ρ − ρ ∗ ∥ 2 2
여기서 K k ℓ = e − ω ℓ τ k Δ ω \mathbf{K}_{k\ell} = e^{-\omega_\ell \tau_k}\Delta\omega K k ℓ = e − ω ℓ τ k Δ ω , ρ ∗ \boldsymbol{\rho}^* ρ ∗ 는 사전 추정값.
NNLS 해결기를 통해 ρ ( ω ) ≥ 0 \rho(\omega) \geq 0 ρ ( ω ) ≥ 0 보장:
A = ( K α 1 ) , b = ( C α ρ ∗ ) \mathbf{A} = \begin{pmatrix} \mathbf{K} \\ \alpha\mathbf{1} \end{pmatrix}, \quad \mathbf{b} = \begin{pmatrix} \mathbf{C} \\ \alpha\boldsymbol{\rho}^* \end{pmatrix} A = ( K α 1 ) , b = ( C α ρ ∗ )
J α = ∥ A ρ − b ∥ 2 2 J_\alpha = \|\mathbf{A}\boldsymbol{\rho} - \mathbf{b}\|_2^2 J α = ∥ A ρ − b ∥ 2 2 를 최소화하되 ρ ≥ 0 \boldsymbol{\rho} \geq 0 ρ ≥ 0 제약.
수치 감산 회피 : Schwinger 함수 방법을 통해 직접 수치 감산 연산 회피물리적 제약 : 양성 제약을 강제하여 스펙트럼 밀도의 물리적 합리성 보장매개변수 선택 : Morozov 차이 원리를 사용하여 정규화 매개변수 α \alpha α 자동 선택로그 샘플링 : 주파수 영역에서 로그 균등 샘플링 사용, 광범위한 주파수 범위 [ 10 − 5 , 10 5 ] [10^{-5}, 10^5] [ 1 0 − 5 , 1 0 5 ] GeV 커버Yamanaka 등Phys. Rev. D 102, 054507 (2020) 이 발표한 SU(2) 순수 Yang-Mills 이론 데이터 사용:
β \beta β 격자 체적 배치 수 N N N a σ a\sqrt{\sigma} a σ a [ Λ − 1 ] a[\Lambda^{-1}] a [ Λ − 1 ] 2.1 10 3 × 12 10^3 \times 12 1 0 3 × 12 1,000,000 7 0.608(16) 0.356(27) 2.2 12 4 12^4 1 2 4 9,999,990 7 0.467(10) 0.273(20) 2.3 14 3 × 16 14^3 \times 16 1 4 3 × 16 4,100,000 9 0.3687(22) 0.216(15) 2.4 16 3 × 24 16^3 \times 24 1 6 3 × 24 2,030,000 13 0.2660(21) 0.156(11) 2.5 20 3 × 24 20^3 \times 24 2 0 3 × 24 520,000 13 0.1881(28) 0.110(8)
기저 상태 질량 추정의 정확성(전통적 방법과 비교) 스펙트럼 밀도 함수의 물리적 합리성(양성, 질량 간격 등) 여기 상태 식별 능력 전통적인 큰 시간 지수 감쇠 분석 방법 최대 엔트로피 방법(장난감 모형 검증에서) 주파수 샘플링 점 수: N ω = 1000 N_\omega = 1000 N ω = 1000 주파수 범위: [ 10 − 5 , 10 5 ] [10^{-5}, 10^5] [ 1 0 − 5 , 1 0 5 ] GeV(로그 균등 분포) 정규화 매개변수는 Morozov 차이 원리로 결정 원본 미도포 데이터 사용 중간자 붕괴 모형에서의 검증은 데이터 포인트 증가 및 잡음 감소에 따라 재구성 질량이 개선됨을 보여준다. 이 방법은 최대 엔트로피 방법과 성능이 상당하다.
β \beta β 전통적 방법 m ϕ m_\phi m ϕ /GeV 스펙트럼 방법 m ϕ m_\phi m ϕ /GeV 상대 차이 2.1 1.34 1.382 +3.1% 2.2 1.43 1.444 +1.0% 2.3 1.48 1.401 -5.3% 2.4 1.53 1.505 -1.6% 2.5 1.63 1.392 -14.6%
기저 상태 질량 추정은 전통적 방법과 높은 일치도를 보이며, 대부분의 경우 차이가 5% 이내.
β = 2.4 \beta = 2.4 β = 2.4 를 제외한 모든 조건에서 3.3-3.9 GeV 범위에 첫 번째 여기 상태 존재 표시:
β \beta β 기저 상태 피크 위치 (GeV) 여기 상태 피크 위치 (GeV) 2.1 1.382 3.712 2.2 1.444 3.354 2.3 1.401 3.941 2.4 1.505 2.670, 6.118 2.5 1.392 3.798
양성 제약 없음 : 현저한 음수 진동 및 과적합 현상 발생양성 제약 있음 : 적외선 진동 제거, 물리적으로 합리적인 스펙트럼 밀도 생성질량 간격 : 모든 스펙트럼 밀도는 약 1.4 GeV의 기저 상태에 해당하는 명확한 질량 간격 표시여기 상태 구조 : 여기 상태 피크는 상대적으로 넓으며, 더 나은 신호 품질 필요 가능격자 효과 : β = 2.4 \beta = 2.4 β = 2.4 의 비정상 거동은 유한 격자 간격의 UV 차단과 관련 가능격자 QCD 글루볼 스펙트럼학 : 전통적 상관함수 분석 방법 사용스펙트럼 함수 추출 : 최대 엔트로피 방법, 베이지안 추론 등Källén-Lehmann 표현 : 격자장론에서의 응용격자 QCD 글루볼 연구의 기초 위에 구축, 특히 Yamanaka 등의 연구 스펙트럼 함수 추출 방법 확장, Tikhonov 정규화의 대체 방안 제공 최대 엔트로피 방법과 상호 보완 전통적 방법의 도포 기법 의존성 회피 여기 상태 정보 제공 수치 구현 상대적으로 간단 직접 방법 유효성 : Tikhonov 정규화 기반의 스펙트럼 밀도 추출 방법은 글루볼 질량 추정에서 양호한 성능 발휘결과 일치성 : 전통적 방법의 기저 상태 질량 추정과 높은 일치도여기 상태 정보 : 여기 상태 구조 성공적 식별로 글루볼 스펙트럼학에 추가 통찰력 제공기술적 장점 : 도포 기법 불필요, 분석 절차 단순화데이터 품질 의존성 : 방법은 통계 오차에 민감하여 고품질 데이터 필요여기 상태 불확실성 : 여기 상태 피크는 상대적으로 넓으며 위치에 큰 변화유한 체적 효과 : 유한 격자 체적의 영향 충분히 논의되지 않음UV 차단 효과 : 격자 간격의 영향 추가 연구 필요다른 양자수의 글루볼 상태로 확장 SU(3) 규범 이론에서의 응용 동역학적 페르미온 포함 경우 여기 상태 분해능 개선 기법 방법 혁신 : 전통적 격자 분석의 효과적인 대체 방안 제공기술적 엄밀성 : 수치 역변환의 병적 성질 올바르게 처리충분한 검증 : 장난감 모형과 실제 데이터를 통한 이중 검증물리적 통찰 : 여기 상태 정보 제공으로 글루볼 스펙트럼에 대한 이해 풍부화실용적 가치 : 도포 기법 회피로 분석 절차 단순화오차 분석 부족 : 스펙트럼 밀도의 통계 오차 분석 상대적으로 단순체계적 오차 : 방법 자체가 도입하는 체계적 오차 충분히 논의되지 않음매개변수 민감성 : 정규화 매개변수 선택에 대한 민감성 분석 미흡물리적 해석 : 여기 상태 피크 폭의 물리적 해석 강화 필요이론적 기여 : 격자 QCD 스펙트럼학 분석에 새로운 도구 제공방법적 가치 : 다른 강자 스펙트럼학 연구로 확대 가능기술적 영향 : Tikhonov 정규화의 격자장론에서의 성공적 응용높은 통계 정확도의 격자 QCD 데이터 분석 여기 상태 정보가 필요한 스펙트럼학 연구 도포 기법 사용이 불편한 경우 글루볼 등 순수 규범 상태 연구 본 논문은 주로 다음의 중요한 연구에 기반:
Yamanaka et al., Phys. Rev. D 102, 054507 (2020) - SU(2) 글루볼 데이터 제공 Asakawa et al., Prog. Part. Nucl. Phys. 46, 459 (2001) - 최대 엔트로피 방법 종설 Dudal et al., Nucl. Phys. B 952, 114912 (2020) - 스펙트럼 표현 방법의 선행 연구 이 논문은 격자 QCD 글루볼 스펙트럼학 분야에서 가치 있는 새로운 방법을 제공하며, 일부 기술적 세부사항에서 개선의 여지가 있지만 기본 사상과 초기 결과는 모두 전망이 밝으며, 해당 분야의 추가 발전을 위한 견고한 기초를 마련했다.