Efficient calculation of two-neutrino double-beta-decay nuclear matrix elements

기술 요약: 두 중성미자 이중 베타 붕괴 핵 행렬 요소의 효율적 계산

문제 정의
신뢰할 수 있는 핵 행렬 요소(NME)는 이중 베타 붕괴($\beta\beta$) 실험을 해석하고 측정된 반감수를 근본적인 약한 상호작용 물리학과 연결하는 데 매우 중요하다. 두 중성미자 모드($2\nu\beta\beta$)는 여러 핵에서 관측된 표준 모델 과정이지만, 무중성미자 모드($0\nu\beta\beta$)는 렙톤 수 위반의 주요 징후로 남아 있다. 최근의 통계적 껍질 모델 분석은 $2\nu\beta\beta$와 $0\nu\beta\beta$ NME 사이의 강한 상관관계를 나타내며, 이는 $2\nu\beta\beta$의 정확한 계산이 후자를 제약할 수 있음을 시사한다.

그러나 $2\nu\beta\beta$ NME($M_{2\nu}$)를 계산하는 것은 상당한 계산적 과제를 안겨준다. 표준적인 접근 방식은 홀수-홀수(odd-odd) 핵의 모든 $1^+$ 중간 상태에 대한 명시적인 합산을 요구한다. 거대 모델 공간(예: 무거운 핵 또는 확장된 가전자 공간인 $sd-pf$)에서는 기여하는 상태의 수가 지나치게 많아지며, 실행 합(running sum)이 느리게 또는 비단조적으로 수렴할 수 있다. 많은-체 해밀토니안의 고유 벡터를 얻기 위한 전체 대각화는 종종 계산적으로 불가능하며, 특히 $10^{28}$년까지의 민감도를 목표로 하는 차세대 실험의 경우 더욱 그러하다.

방법론
본 논문은 $M_{2\nu}$를 계산하기 위해 전체 대각화를 요구하지 않는 란초스 반복(Lanczos iterations) 기반의 강도 함수(strength-function) 방법의 개선된 구현을 제시한다. 이 접근 방식은 문헌 [24]의 형식론을 따르며 문헌 [25]에 상세히 기술된 수치적 처방을 활용한다.

1. 란초스 반복: 이 방법은 초기($0^+_i$) 또는 최종($0^+_f$) 바닥 상태에 감마-지티(GT) 연산자를 적용하여 구축된 GT "도어웨이(doorway)" 상태를 사용하여 란초스 반복을 시작한다.
2. 삼중 대각 행렬: 많은-체 해밀토니안의 전체 고유 벡터를 저장하는 대신, 알고리즘은 훨씬 작은 삼중 대각 란초스 행렬을 구축한다. NME는 이 작은 행렬의 고유값과 고유 벡터를 합산함으로써 근사된다.
3. 계수: 확장 계수($g^\pm_m$)는 란초스 벡터와 중간 상태와 관련된 에너지 분모로부터 유도된다.
4. 효율성: 이 접근 방식은 전체 해밀토니안을 대각화하는 데 드는 저장 및 처리 비용을 피한다. 이는 란초스 벡터 자체를 사용하여 반복적으로 수렴을 추적함으로써, 이전의 강도 함수 구현에 비해 약 2배의 계산 효율성 증가를 제공한다.
5. 고차 항으로의 확장: 본 프레임워크는 위상 공간 인자(phase-space factors) 및 관련 두 중성미자 전자 포획 관측량의 테일러 전개 처리에 필요한 고차 NME($M_{2\nu-3}$ 및 $M_{2\nu-5}$)를 계산하도록 확장되었다. 이러한 항들은 에너지 분모의 더 높은 거듭제곱을 포함하므로 저에너지 상태에 더 민감하다.

주요 기여 및 결과
저자는 표준 유효 해밀토니안(예: GXPF1A, KB3G, JUN45, GCN2850, GCN5082, SVD)과 그에 대응하는 가전자 공간을 사용하여 몇 가지 실험적으로 중요한 동위원소($^{48}\text{Ca}$, $^{76}\text{Ge}$, $^{82}\text{Se}$, $^{128,130}\text{Te}$, $^{136}\text{Xe}$)에 이 방법을 적용하였다.

• 벤치마킹: 란초스 확장을 통해 얻은 "실행 합(running sums)"을 수백 또는 수천 개의 중간 $1^+$ 상태에 대한 명시적 합과 비교하여 방법을 엄격하게 테스트하였다.
  • $^{48}\text{Ca}$ (pf-shell)의 경우, 25회 반복된 란초스 결과는 1000개 상태에 대한 직접 합을 재현하며, 급격한 상승, 강도 영역의 변화, 그리고 최종적인 안정적 평탄부(plateau)를 포착한다.
  • $^{82}\text{Se}$ 및 $^{136}\text{Xe}$와 같이 직접 대각화가 어려운 더 무거운 핵의 경우, 란초스 방법은 직접 대각화(예: $^{136}\text{Xe}$에 대해 2000개 이상의 상태)의 수렴 동작을 재현한다.
  • 이 방법은 GT 거대 공명 영역의 기여를 성공적으로 포착하며, $2\nu\beta\beta$ 행렬 요소에서 나타나는 비자명한 상쇄(cancellations)를 처리한다.
• 고차 항: $^{48}\text{Ca}$에 대한 $M_{2\nu-3}$ 및 $M_{2\nu-5}$의 벤치마크는 명시적 합과 매우 잘 일치함을 보여준다. 이는 주요 NME를 위해 사용된 동일한 란초스 계수가 에너지 분모의 더 높은 역수 거듭제곱에 대해서도 일반화될 수 있음을 확인시켜 준다.
• 쿼칭 인자(Quenching Factors): 본 연구는 경험적 행렬 요소를 재현하기 위한 다양한 해밀토니안에 대한 최적의 GT 쿼칭 인자($q_{opt}$)를 결정한다. 결과는 이러한 인자들이 보편적인 상수가 아니라 특정 상호작용 및 가전자 공간에 의존하는 현상론적 교정임을 강조한다.

의의 및 주장
본 논문은 제안된 란초스 강도 함수 방법이 직접적인 합산이 계산적으로 매우 비싸거나 불가능한 모델 공간에서 $M_{2\nu}$를 평가하기 위한 실용적이고 효율적인 해결책을 제공한다고 주장한다.

• 정확도: 이 방법은 벤치마킹이 가능한 경우, 큰 상쇄와 고에너지 상태의 기여를 포함한 명시적 합의 정확도를 보존한다.
• $0\nu\beta\beta$에 대한 유용성: $2\nu\beta\beta$와 $0\nu\beta\beta$ NME 사이의 강한 통계적 상관관계를 고려할 때, 신뢰할 수 있고 효율적인 $2\nu\beta\beta$ 계산은 $0\nu\beta\beta$ 예측에 유리한 유효 해밀토니안 및 연산자 재규격화를 식별하는 데 도움을 줄 수 있다.
• 확장성: 이 방법의 효율성은 이전에 직접 대각화를 통해 접근할 수 없었던 확장된 모델 공간(예: $2\hbar\omega$ 들뜸 포함)의 탐색을 가능하게 한다.
• 향emerging 전망: 저자는 이 프레임워크가 측정된 반감수로부터 도출된 통계적 샘플과 현상론적 제약을 제공함으로써, 부상하는 ab initio 계산에 대한 보완적인 대규모 벤치마크 역할을 할 수 있다고 제안한다.

본 연구는 새로운 실험적 징후를 제안하거나 $0\nu\beta\beta$ 문제를 직접 해결한다고 주장하는 것이 아니라, 효율적인 $2\nu\beta\beta$ NME 계산을 무중성미자 붕괴 탐색의 이론적 불확실성을 줄이기 위한 필수적인 도구로 위치시킨다.