Solution to the uncertainty problem of nuclear matrix element for neutrinoless double-$\beta$ decay

기술 요약: 뉴트리노 무중성 중성자 베타 붕괴의 핵 행렬 요소 불확실성 문제에 대한 해결책

문제 정의
뉴트리노 무중성 이중 베타 붕괴($0\nu\beta\beta$)의 검출은 뉴트리노의 마요라나(Majorana) 성질을 확인하고 표준 모형을 확장할 수 있는 뉴트리노 물리학의 주요 목표이다. 그러나 실험적 반감기 한계로부터 유효 뉴트리노 질량($\langle m_\nu \rangle$)을 추출하는 과정은 핵 행렬 요소(Nuclear Matrix Elements, NMEs)의 정확도에 결정적으로 의존한다. 40년 넘게 이론적 NME 계산은 큰 불확실성을 보여왔으며, 서로 다른 핵 구조 모델(예: QRPA, Shell Model, IBA)의 결과값들은 넓은 범위에 걸쳐 분포되어 있다. 이러한 분산은 뉴트리노 물리학의 발전을 저해해 왔으며, NME의 불확실성은 곧 유도된 뉴트리노 질량의 불확실성으로 직결된다.

방법론
저자들은 $0\nu\beta\beta$ 붕괴와 이중 뉴트리노 베타 붕괴($2\nu\beta\beta$) 사이의 관계를 활용하여 이 불확실성을 줄이는 방법을 제안한다. 이들의 접근 방식의 핵심은 다음과 같다:

1. 섭동된 전이 연산자 (Perturbed Transition Operator): 저자들은 NME가 섭동된 전이 연산자를 사용하여 도출되는 계산 프레임워크를 활용한다. 여기에는 고차 항(leading-order인 1차 Rayleigh-Schrödinger 섭동을 넘어선 항)인 정점 수정(vertex corrections) 및 이체 전류 항(two-body current terms, RS two-body current terms이라 지칭됨)이 포함된다.
2. 유효 축방향 벡터 결합 상수 ($g_A^{eff}$): 저자들은 leading-order NME에 완전히 포함되지 않은 다체 효과(many-body effects)를 포괄하는 유효 축방향 벡터 결합 상수 $g_A^{eff}$를 정의한다. 이 파라미터는 leading-order NME와 $g_A^{eff}$를 사용해 계산한 반감기를 bare coupling인 $g_A^{bare}$와 섭동된 NME를 사용해 계산한 반감기와 동일하게 만듦으로써 결정된다.
3. 유사성 관계 (Similarity Relation): 저자들의 이전 연구에서 도출된 핵심 발견은, 두 과정 사이의 거대한 뉴트리노 운동량 차이($\sim 100$ MeV/c 대 $\sim 2-3$ MeV/c)에도 불구하고, $0\nu\beta\beta$와 $2\nu\beta\beta$ 모드에 대한 이론적 $g_A^{eff}$ 값이 매우 유사하다는 것이다. 이러한 유사성은 $g_A^{eff}$를 정의하는 비율의 분자와 분모에 존재하는 뉴트리노 포텐셜이 대략적으로 상쇄되기 때문에 발생한다.
4. 재규격화 절차 (Renormalization Procedure): $2\nu\beta\beta$ 반감기는 실험적으로 측정되므로, 저자들은 $2\nu\beta\beta$ 모드에 대한 $g_A^{eff}$ 값($g_{A,2\nu}^{eff}$)을 실험값과 일치하도록 결정한다. 그 후, 이 값을 $g_A^{eff}$ 값들의 비율로부터 유도된 재규격화 인자를 사용하여 $0\nu\beta\beta$ 계산에 적용한다. 이를 통해 고차 섭동 항의 완전한 수렴 없이도 $0\nu\beta\beta$ NME에 비섭동적 다체 효과를 포함할 수 있다.

적용 및 결과
이 방법은 $^{136}\text{Xe}$, $^{76}\text{Ge}$, $^{130}\text{Te}$, $^{110}\text{Pd}$의 네 가지 동위원소에 적용되었다. 저자들은 다양한 그룹으로부터 다양한 핵 구조 방법(QRPA, Shell Model, IBA)을 사용한 NME를 수집하고 이들의 재규격화 절차를 적용하였다.

• 분산 감소: 가장 중요한 결과는 NME의 분산이 극적으로 감소했다는 점이다. $^{136}\text{Xe}$의 경우, 수정된 NME( $g_A^2$ 가 곱해진 값)의 범위가 기존 구간의 4분의 1로 줄어들었다. $^{76}\text{Ge}$의 경우, 변형된 QRPA 결과의 분산이 89% 감소하였으며, 모든 방법론에 대한 전체적인 분산은 약 40~50% 감소하였다.
• NME 값의 하락: 많은 경우, 수정된 NME 값은 이전 계산들에 비해 낮아졌다. 예를 들어, $^{76}\text{Ge}$의 평균값은 2.7 단위 낮아졌다.
• 유효 뉴트리노 질량 ($\langle m_\nu \rangle$): 수정된 NME와 현재의 실험적 반감기 하한치를 사용하여 저자들은 유도된 $\langle m_\nu \rangle$를 계산하였다.
  • $^{136}\text{Xe}$의 경우, 결과는 (가장 가벼운 뉴트리노 질량이 $< 10$ meV라고 가정할 때) 추출된 $\langle m_\nu \rangle$가 아직 역 질량 계층 구조(inverted-mass-hierarchy) 영역에 도달하지 못함을 나타낸다.
  • $^{76}\text{Ge}$ 및 $^{130}\text{Te}$에서도 유사한 경향이 관찰되었다. 수정된 계산은 현재의 실험적 한계가, 일부 이전의 더 높은 NME 추정치가 시사했던 것과는 달리, 역 질량 계층 구조 영역을 탐사하기에 충분하지 않을 수 있음을 시사한다.
• 계통적 특성 (Systematics): 본 연구는 $g_{A,0\nu}^{eff}$와 $g_{A,2\nu}^{eff}$의 비율이 여러 동위원소에 걸쳐 1에 가깝다(1차에서 수렴함)는 것을 확인하여 방법론의 적용 가능성을 입증하였다. 수정 항(정점 및 이체 전류)에 대한 분석은 질량수($A$) 의존성, 특히 이체 전류 항에 대한 강한 의존성을 보여주었으며, 이는 서로 다른 핵에서 관찰되는 상쇄 정도의 차이(예: $^{76}\text{Ge}$에서의 강한 상쇄 vs $^{136}\text{Xe}$에서의 큰 보정)를 설명한다.

의의 및 주장
본 논문은 $0\nu\beta\beta$ NME의 오랜 불확실성 문제에 대한 "가능한 해결책"을 제시한다고 주장한다. 이 연구의 의의는 다음과 같다:

1. 현상론적 보정 (Phenomenological Calibration): 고차 항의 이론적 수렴에만 의존하는 대신, 실험적으로 알려진 $2\nu\beta\beta$ 붕괴를 기준으로 보정함으로써 비섭동적으로 다체 효과를 포함하는 방법을 제공한다.
2. 보편성 (Universality): leading-order NME가 일관되게 계산된다면, 서로 다른 그룹이 사용하는 서로 다른 핵 구조 모델의 결과에 적용 가능한 방법이다.
3. 불확실성 감소: NME 값의 분산을 크게 좁힘으로써, 실험 데이터로부터 유효 뉴트리노 질량을 추출할 때 발생하는 이론적 불확실성을 줄인다.
4. 현실적인 전망: 결과는 이전의 더 높은 NME 추정치들과 비교했을 때 수정된 NME가 더 낮은 붕괴율을 나타내므로, 현재 및 가까운 미래의 실험에서 역 질량 계층 구조의 관측 가능성에 대해 더 보수적인 전망을 제시한다.

저자들은 이 접근 방식이 두 붕괴 모드의 유효 결합 상수 사이의 유사성 관계에 의해 뒷받침되는, NME 불확실성 해결을 향한 "독특하고 새로운 방향"을 제시한다고 결론짓는다.