Proton and neutron electromagnetic form factors from lattice QCD in the continuum limit

기술 요약: 연속체 극한에서의 양성자 및 중성자 전자기 폼 팩터

문제 정의
양성자와 중성자의 전자기 폼 팩터는 강입자의 내부 구조를 설명하는 근본적인 물리량이다. 양성자 폼 팩터는 실험적으로 잘 측정되어 있는 반면, 중성자 폼 팩터는 자유 중성자 표적의 부재로 인해 전자-중수소 또는 전자-헬륨 산란으로부터 간접적으로 추출해야 하므로 정밀도가 상대적으로 낮다. 격자 양자 색역학(Lattice QCD)은 이러한 양들을 계산하기 위한 제일 원리적 접근법을 제공한다. 그러나 기존의 격자 계산들은 종종 한계에 직면해 왔다. 즉, 연속체 외삽(continuum extrapolation)이 결여되었거나, 비물리적인 쿼크 질량으로부터의 카이랄 외삽(chiral extrapolation)에 의존했거나, 혹은 디스커넥티드 기여(isoscalar 항)를 충분한 정밀도로 포함하지 못했다. 본 연구는 물리적인 파이온 질량에서 직접 폼 팩터를 계산하고, 물리적 쿼크 질량을 가진 앙상블을 사용하여 연속체 극한을 취함으로써 이러한 공백을 해결한다.

방법론
저자들은 $N_f=2+1+1$ twisted mass clover-improved 페르미온을 이용한 격자 QCD를 채택하였다. 시뮬레이션은 격자 간격 $a=0.080$ fm, $0.068$ fm, $0.057$ fm인 세 가지 게이지 앙상블을 사용하며, 이들은 모두 물리적 쿼크 질량(the "physical point")에 맞춰져 있다.

• 상관 함수(Correlation Functions): 뉴클레온 행렬 요소는 2점 및 3점 상관 함수로부터 추출된다. 들뜬 상태 오염(excited state contamination)을 억제하기 위해, 저자들은 0.5 fm에서 1.5 fm에 이르는 다양한 싱크-소스 시간 간격을 활용한다.
• 들뜬 상태 분석(Excited State Analysis): 엄격한 다중 상태 피팅(multi-state fit) 분석이 수행된다. 표준적인 플래토 피팅(plateau fits)과 달리, 이 접근 방식은 2점 함수와 3점 함수에서 서로 다른 들뜬 상태 에너지를 허용하며, 이는 전자기 전류가 보간장(interpolating field)과는 다르게 다입자 상태(예: 파이온-뉴클레온 상태)에 결합할 수 있음을 인지한 것이다.
• 디스커넥티드 기여(Disconnected Contributions): 쿼크 루프(disconnected diagrams)를 포함하는 아이소스칼라 기여는 통계적 노이즈를 완화하기 위해 고통계 기법을 사용하여 계산된다. 여기에는 low-mode deflation, color-spin dilution, 그리고 hierarchical probing이 포함된다. 디스커넥티드 기여는 로컬 벡터 전류를 사용하여 계산되며, 이는 재규격화가 필요하다. 반면, 커넥티드 기여는 보존된 벡터 전류를 사용한다.
• 재규격화(Renormalization): 맛(flavor) 비단일 및 단일 벡터 연산자에 대한 재규격화 인자는 RI'-MOM 스킴을 사용하여 결정된다. 분석 결과 싱글렛과 논싱글렛 인자가 일치함을 확인하였으며, 이는 아이소스칼라 전류의 재규격화를 단순화한다.
• 관측량 추출(Extraction of Observables): 폼 팩터는 디폴(dipole) 및 $z$-expansion 안사츠(Ansätze)를 모두 사용하여 추출된다. 저자들은 $Q^2$ 의존성과 연속체 외삽($a^2$에 대하여)을 동시에 피팅하는 "one-step" 접근 방식을 사용한다. 모델 평균화(model averaging)를 사용하여 서로 다른 피팅 범위와 파라미터화의 결과를 결합하며, 최종 통계적 및 계통적 오차를 도출하기 위해 아카이케 정보 기준(AIC)에 따라 가중치를 부여한다.

주요 기여

1. 물리적 지점에서의 연속체 극한: 본 연구는 물리적 파이온 질량을 가진 앙상블을 사용하여 양성자 전자기 폼 팩터에 대한 연속체 극한을 달성한 최초의 계산 중 하나로, 이를 통해 카이랄 외삽의 필요성을 제거하였다.
2. 포괄적인 들뜬 상태 제어: 다양한 시간 간격을 가진 다중 상태 피팅을 사용하고 2점 및 3점 함수에서 별개의 들뜬 상태 에너지를 허용함으로써, 기존 연구보다 더 견고한 바닥 상태 행렬 요소 추출을 제공한다.
3. 디스커넥티드 기여 포함: 아이소스칼라 폼 팩터에 대한 고통계 디스커넥티드 기여 분석을 포함하며, hierarchical probing 및 deflation과 같은 고급 노이즈 감소 기법을 사용하여 커넥티드 기여와 대등한 정밀도를 달성하였다.
4. 계통 오차 추정: 피팅 범위, 모델 의존성(dipole vs. $z$-expansion), 그리고 연속체 외삽으로부터 발생하는 계통 불확실성에 대한 상세한 분석을 제공하며, 모델 평균화를 통해 이러한 효과를 정량화하였다.

결과
저자들은 연속체 극한에서의 양성자 및 중성자에 대한 전기 및 자기 반경, 그리고 자기 모멘트를 제시한다. 결과는 통계적 및 계통적 오차(후자는 피팅 및 연속체 외사를 포함함)와 함께 인용된다:

• 양성자 전기 반경: $\sqrt{\langle r^2_E \rangle_p} = 0.860(38)(23)$ fm.
• 중성자 전기 반경: $\langle r^2_E \rangle_n = -0.147(48)$ fm$^2$.
• 양성자 자기 반경: $\sqrt{\langle r^2_M \rangle_p} = 0.870(53)(15)$ fm.
• 중성자 자기 반경: $\sqrt{\langle r^2_M \rangle_n} = 0.913(67)(19)$ fm.
• 양성자 자기 모멘트: $\mu_p = 2.849(92)(52)$.
• 중성자 자기 모멘트: $\mu_n = -1.819(76)(29)$.

또한, 저자들은 초미세 구조 분할(hyperfine splitting) 계산에 중요한 제마크(Zemach) 및 프리어(Friar) 반경을 계산한다:

• 양성자 제마크 반경: $r_Z^p = 1.091(34)$ fm.
• 중성자 제마크 반경: $r_Z^n = -0.0405(81)$ fm.

반경 및 자기 모멘트에 대한 결과는 결합된 통계적 및 계통적 오차 범위 내에서 입자 데이터 그룹(PDG) 값과 일치한다. 계산된 양성자 전기 반경은 역사적인 전자-양성자 산란의 더 큰 값들과 뮤온 수소 분광학의 더 작은 값들 모두와 겹치는 불확실성을 가진다.

의의 및 주장
본 논문은 이러한 결과들이 격자 QCD 공식화의 벤치마크를 제공하며, 비섭동적 계산을 실험 데이터와 검증한다고 주장한다. 물리적 지점에서 직접 폼 팩터를 계산하고 연속체 극한을 구함으로써, 본 연구는 특히 중성자와 같이 실험적으로 덜 알려진 강입자 구조 양들을 평가할 수 있는 격자 QCD의 능력을 입증한다.

저자들은 자신들의 결과가 실험 데이터와 일치하지만, 양성자 전기 반경의 정밀도는 현재 통계적 및 계통적 오차에 의해 제한되어 있다고 겸허히 언급한다. 그들은 개별 앙상블의 정확도를 높이고, 더 작은 격자 간격 및 더 큰 물리적 부피(더 낮은 $Q^2=0$ 근처의 운동량을 접근하기 위함)를 가진 앙상블을 포함하는 것이 오차를 더욱 줄이고, 잠재적으로 양성자 반경에 대한 서로 다른 실험적 결정 사이의 불일치를 해결하는 데 필요할 것이라고 제안한다. 이 작업은 또한 격자 QCD에서 디스커넥티드 기여를 포함하고 엄격한 들뜬 상태 분석을 수행하는 방법론의 타당성을 입증하는 역할을 한다.