LUNAR: a Monte Carlo generator for bound-nucleon decay in liquid argon

기술 요약: LUNAR – 액체 아르곤 내 결합된 핵자의 붕괴를 위한 몬테카를로 생성기

문제 정의
액체 아르곤 시점 투영 챔버(LArTPC), 특히 DUNE(Deep Underground Neutrino Experiment)에서의 핵자 붕decay 탐색에는 결합된 핵의 환경이 붕괴 생성물의 운동학을 어떻게 변화시키는지에 대한 정밀한 정량적 기술이 필요하다. 자유 핵자와 달리, 결합된 핵자는 고유의 운동량(페르미 운동)을 가지며 결합 에너지(제거 에너지)로 인해 자유 질량보다 낮은 상태로 존재한다. 이러한 효과는 딸 입자의 스펙트럼을 넓히고 가용 에너지를 이동시키며, 잠재적으로 특정 붕괴 모드를 불가능하게 만들 수 있다. 또한, 붕괴 메존은 잔류 핵을 통과해야 하며, 검출되기 전에 산란, 전하 교환 또는 흡수와 같은 최종 상태 상호작용(FSI)을 겪어야 한다. GENIE, NuWro, GiBUU와 같은 범용 뉴트리노 이벤트 생성기들은 일반적인 핵 모델을 포함하고 있지만, 이들은 뉴트리노-핵 상호작용을 위해 설계되었으며 전용 핵자 붕괴 연구를 수행하기에는 너무 무겁거나 내부가 불투명할 수 있다. 따라서 결합된 핵자 붕괴의 독특한 운동학에 특화된 가볍고 투명하며 확장 가능한 도구에 대한 구체적인 필요성이 존재한다.

방법론: LUNLUR 생성기
저자들은 $^{40}$Ar(아르곤-40)에 결합된 양성자와 중성자의 이체 붕괴를 위해 특별히 설계된 헤더 전용(header-only) C++17 몬테카를로 생성기인 LUNAR(Lancaster University Nucleon-decay ARgon generator)를 제시한다. 이 생성기는 다음과 같은 모듈형 파이프라인을 통해 작동한다:

1. 초기 상태 샘플링:

   • 운동량: 결합된 핵자는 10가지 선택 가능한 운동량 분포($n(p)$) 중 하나에서 추출된다. 이는 단순한 평균장 페르미 가스(전역 및 국소)부터 단거리 상관관계(SRC), Bodek–Ritchie 가스, 조화 진동자 껍질 모델, 그리고 가우시안 분포를 포함하는 모델까지 다양하다. 두 가지 옵션은 스펙트럴 함수를 사용하는데, 하나는 Jefferson Lab(benhar)의 테이블 그리드를 사용하고 다른 하나는 해석적 구성(ankowski)을 사용한다.
   • 결합/오프-쉘(Off-shell) 운동학: 핵자는 세 가지 제거 에너지 규정(기본값인 운동량 의존 광학 퍼텐셜, 상수 분리 에너지, 또는 껍질 모델 레벨 스킴) 중 하나를 사용하여 오프-쉘 질량 상태로 결합된다(상관된 핵자의 경우, 제거 에너지는 파트너 핵자의 반동에 의해 결정된다).
   • 생성기는 붕괴가 운동학적으로 허용되는지 확인한다. 만약 오프-쉘 불변 질량 $W$가 두 딸 입자 질량의 합보다 낮아지면, 핵자는 재샘플링된다.

2. 붕괴:

   • 이체 붕계($N \to d_1 + d_2$)는 오프-쉘 질량 $W$에 기반한 상대론적 운동학을 사용하여 핵자의 정지 좌표계에서 수행된다.
   • 생성된 4-벡터는 실험실 좌표계로 부스트(boost)된다.
   • LUNAR는 양성자 및 중성자 붕괴를 모두 다루며, 초대칭 선호 모드(예: $p \to K^+\bar{\nu}$)와 대통일 이론(GUT) 모드(예: $p \to e^+\pi^0$)를 포함한 14가지 표준 채널을 지원한다.

3. 최종 상태 상호작용 (FSI):

   • 붕괴 메존은 NuWro의 kaskada 접근 방식을 모델로 한 준고전적 핵 내부 카스케이드(intranuclear cascade)를 통해 핵 밖으로 전파된다.
   • 카스케이드는 Woods-Saxon 밀도를 통해 하드론을 추적하며, 하드론-핵 단면적(파이온/핵자의 경우 Metropolis 테이블, 스트레인지/에타 메존의 경우 보간된 데이터 사용)을 기반으로 상호작용 확률을 계산한다.
   • 과정에는 탄성 산란, 전하 교환, 흡수 및 다중 메존 생성이 포함된다. 나가는 핵자에 대해서는 파울리 차단(Pauli blocking)이 적용된다.
   • "형성 영역(formation zone)" 옵션을 활성화할 수 있으며, 이는 색 단일체(color singlet)가 형성되는 데 걸리는 시간을 시뮬레이션하기 위해 새로 생성된 메존이 상호작용하기 전 일정 거리 $L_f$ 동안 자유롭게 이동하도록 허용한다.

주요 기여

• 모듈형 설계: 코드는 독립적인 모듈(운동량, 결합, 운동학, 카스케이드)로 구조화되어 있어, 핵심 설계를 변경하지 않고도 사용자가 핵 모델이나 붕괴 채널을 교체할 수 있다. 이 코드는 C++17 표준 라이브러리 외에 외부 의존성 없이 오픈 소스 소프트웨어로 공개된다.
• 핵 효과의 분리: 생성기는 페르미 운동(스펙트럼 폭을 제어)과 핵 결합(스펙트럼 중심값을 제어)을 분리하여 분석할 수 있게 한다.
• 채널별 FSI 정량화: 일반적인 효율 가정에서 벗어나, 다양한 붕괴 모드에 따른 메존 생존 확률에 대한 상세한 내역을 제공한다.

결과
저자들은 LUNAR를 사용하여 핵 효과가 관측 가능한 메존 스펙트럼과 DUNE(20 kt 피두셜 질량 및 400 kt·yr 노출 가정)의 예상 이벤트 수율에 미치는 영향을 분석하였다:

1. 스펙트럼 형태: 관측 가능한 딸 입자 운동량 스펙트럼의 폭은 선택된 운동량 모델(특히 고운동량 꼬리 함량)에 의해 결정되며, 중심값은 결합 모델에 의해 결정된다. 이 두 불확실성은 대체로 직교한다.
2. 최종 상태 상호작용 (FSI):
   • $K^+$ 모드: 핵은 $K^+$ 메존에 대해 거의 투명하다. 약 99%가 상호작용 없이 탈출하며, FSI 효과는 무시할 수 있는 수준이다(운동량의 몇 퍼센트 정도의 이동).
   • 파이온, $\eta$, 및 안티카온 모드: 이 채널들은 상당한 억제를 겪는다. 파이온, $\eta$ 메존, 안티카온의 약 40–60%가 흡수되거나 전하 교환을 일으킨다. 예를 들어, $p \to e^+\pi^0$ 수율은 흡수로 인해 약 50% 감소한다.
   • 형성 영역: 형성 영역을 활성화하면 강하게 상호작용하는 파이온의 생존 수율이 크게 증가한다(예: $p \to e^+\pi^0$의 경우 $L_f=1.0$ fm에서 49%에서 81%로 증가). 이는 비-카온 채널에서 주요 계통 불확실성임을 보여준다.
3. 이벤트 수율: 현재의 Super-Kamiokande 제한을 DUNE 예상 수율로 변환하면, "골든(golden)" 채널인 $p \to K^+\bar{\nu}$는 여전히 깨끗한 상태를 유지하지만, 다른 모드들(예: $n \to \bar{\nu}K^0$, $p \to \bar{\nu}\pi^+$)은 수십에서 수백 개의 이벤트를 볼 수 있다. 그러나 이러한 비-카온 모드들의 경우, FSI 억제가 지배적인 계통 효과이며, 이는 핵 운동량 모델의 선택으로 인한 $\pm 10\%$의 편차보다 더 큰 효과인 약 절반의 비율로 수율을 감소시킨다.

의의 및 주장
본 논문은 LUNAR가 범용 뉴트리노 생성기의 복잡성을 피하면서 핵자 붕괴 커뮤니티를 위한 전용의 투명한 도구를 제공함으로써 중요한 공백을 메운다고 주장한다. 주요 의의는 다음과 같다:

• "골든" 채널의 검증: $p \to K^+\bar{\nu}$가 핵 FSI로부터 거의 면역되어 있어 액체 아르곤 검출기에 가장 깨끗한 탐사 도구임을 재확인한다.
• 계통 불확실성 정량화: 모든 표준 붕괴 모드에 대해, 최종 상태 상호작용이 핵 초기 상태의 불확실성보다 더 큰 영향을 미치는 주요 보정 요인임을 입증한다.
• 유용성: 신호 효율 연구, 계통 불확실성 추정(핵 모델의 확산을 통한), 그리고 핵자 붕괴 물리학의 교육적 탐구를 위한 빠르고 확장 가능한 도구를 제공한다.

저자들은 분석적 운동량 모델이 피팅된 것이 아니라 파라미터화된 것이라는 점, 카스케이드가 대표적인 단면적을 사용한다는 점, 그리고 검출기 반응이 포함되지 않았다는 점을 인정하며 신중한 입장을 유지한다. 그들은 LUNAR를 액체 아르곤 내 핵물리학의 기초적인 도구로서, 커뮤니티에 의해 확장되고 개선될 수 있도록 위치시킨다.