Determining the dynamic deformation of $^{140}$Ce by constraining coupled-channels parameters for fusion

기술 요약: 결합 채널 매개변수 제약을 통한 $^{140}$Ce의 동적 변형 결정

문제 정의
핵 구조 매개변수(사중극자 $\beta_2$, 팔중극자 $\beta_3$, 육중극자 $\beta_4$)를 추출하는 것은 핵물리학의 핵심적인 과제이다. 융합 반응은 결합 채널(coupled-channels, CC) 분석을 통해 핵 변형을 탐구하는 민감한 프로브 역할을 하지만, 현재의 방법론은 상당한 한계에 직면해 있다. 구체적으로, 매개변수 튜닝이 개별 반응 시스템에 국한되는 경우가 많아, 이를 다른 투사체-표적 조합에 적용할 때 일관성이 없거나 비물리적인 결과를 초래한다. 또한, 융합 장벽 분포($D$)를 도출하기 위한 전통적인 이중 미분(double-derivative, DD) 방식은 단계 크기(step size)에 매우 민감하며 고에너지 영역에서 정확도가 떨어져, 추출된 구조 정보의 신뢰성을 제한한다. 진동 핵으로 알려진 $^{140}$Ce의 경우, 이전 연구들은 사용된 실험적 프로브와 이론적 모델에 따라 $\beta_2$가 0.003에서 0.15까지 나타나는 등 현저히 다른 변형 매개변수를 보고해 왔다.

방법론
본 연구는 새로운 실험 데이터, 견고한 분석적 장벽 분포 프레임워크, 그리고 여러 반응 시스템에 걸친 엄격한 통계적 분석을 결합한 체계적인 접근 방식을 채택한다:

1. 실험 데이터 획득:

   • $^{16}$O + $^{140}$Ce: IUAC의 중이온 반응 분석기(HIRA)를 사용하여 새로운 융합 단면적($\sigma_{fus}$) 측정을 수행하였다. 데이터는 Bass 장벽(57–76 MeV)보다 약 12.4% 낮은 지점부터 약 17% 높은 지점까지 수집되었다.
   • $^{36}$S + $^{140}$Ce: 기존 문헌의 융합 데이터를 활용하였다.
   • $^{28}$Si + $^{140}$Ce: 기존의 들뜸 함수(excitation functions)를 검증용 테스트베드로 사용하였다.

2. 장벽 분포 분석:

   • 융합 장벽 분포($D$)는 두 가지 방법, 즉 전통적인 DD 방법과 **가우시안 분석 레시피(Gaussian analytic recipe)**를 사용하여 도출되었다. 후자는 단계 크기에 독립적이고 $\sigma_{fus}$의 큰 불확실성에 대한 탄력성이 높으므로, 완전하고 잘 정의된 스펙트럼을 제공하기 위해 선택되었다.

3. 결합 채널(CC) 계산:

   • 계산은 Woods-Saxon 핵 퍼텐셜을 사용하는 ccfull 코드를 사용하여 수행되었다.
   • 표적 모델링: $^{140}$Ce는 일중양자 사중극자($2^+$) 및 팔중극자($3^-$) 들뜸을 포함하는 순수 진동 핵으로 모델링되었다.
   • 매개변수 공간: $\beta_2$와 $\beta_3$를 0.00에서 0.60까지 체계적으로 변화시켰다.
   • 검증: $^{28}$Si + $^{140}$Ce 시스템은 투사체의 고유 변형(진동체 및 회전체 모두로 취급)과 양의 Q-값 2중자(2n) 전이 채널을 포함하는 수정된 결합 스킴을 통해 분석되었다.

4. 통계적 분석:

   • $\chi^2$ 최소화: 분석적 장벽 분포와 이론적 CC 결과 사이의 차이를 최소화함으로써 변형 매개변수를 추출하였다.
   • 베이지안 분석: $^{16}$O 및 $^{36}$S 시스템 모두에서 $\beta_2$와 $\beta_3$의 사후 확률 분포를 결정하기 위해 마르코프 연쇄 몬테카를로(MCMC) 프레임워크를 사용하였다.
   • 모델 평균화: 두 시스템의 결과를 결합하여 최적의 값과 관련 불확실성을 얻기 위해 독립적인 베이지안 모델 평균화(NBMA)를 적용하였다.

주요 결과

• 장벽 분포: 가우시안 분석 레시피는 $^{36}$S + $^{140}$Ce 시스템에 대해 잘 정의된 장벽 분포를 성공적으로 생성하였으며, $E \simeq 104$ MeV에서 주요 피크와 넓은 고에너지 구조를 보여주었다. 이는 spurious(가짜) 3중 피크 구조를 생성했던 전통적인 DD 방법과 대조되며, 고에너지에서의 DD 방식의 한계를 부각시킨다.
• 추출된 변형 매개변수:
  • $^{16}$O + $^{140}$Ce의 경우: $\beta_2 = 0.06^{+0.07}_{-0.06}$ 및 $\beta_3 = 0.18^{+0.02}_{-0.04}$.
  • $^{36}$S + $^{140}$Ce의 경우: $\beta_2 = 0.09^{+0.02}_{-0.02}$ 및 $\beta_3 = 0.18^{+0.01}_{-0.01}$.
  • NBMA 결합 결과: 본 연구는 $^{140}$Ce의 동적 변형 매개변수를 $\beta_2 = 0.09 \pm 0.03$ 및 $\beta_3 = 0.18 \pm 0.02$로 결정하였다. 이 값들은 두 투사체 시스템 모두에서 일관되게 나타난다.
• $^{28}$Si + $^{140}$Ce에서의 검증:
  • 추출된 매개변수($\beta_2=0.09, \beta_3=0.18$)를 $^{28}$Si + $^{140}$Ce 시스템에 적용했을 때, CC 계산은 융합 들뜸 함수와 분석적 장벽 분포를 모두 성공적으로 재현하였다.
  • 양의 Q-값 2n-pickup 채널의 포함이 하부 장벽 융합 증폭(sub-barrier fusion enhancement)과 장벽 분포의 형태를 재현하는 데 필수적임을 확인하였다.
  • $^{28}$Si를 진동체 또는 강체 회전체로 처리하는 것은 결과에 미미한 차이만을 주었으며, 이는 추출된 표적 매개변수의 견고성을 시사한다.
• 문헌과의 비교: 추출된 $\beta_2$ 값은 Pritychenko 등의 값($\beta_2 = 0.10$)과 밀접하게 일치하지만, $\beta_2 \approx 0.15$를 제안했던 이전의 융합 기반 연구들과는 차이가 있다. 추출된 $\beta_3$는 기존 문헌에 보고된 값보다 크다.

의의 및 주장
본 논문은 이 연구가 단일 시스템 분석에서 흔히 발견되는 모호함을 해결하며, 여러 반응 시스템에 의해 제약될 때 추출된 변형 매개변수의 일관성과 견고성을 입증한다고 주장한다. 또한, 본 연구는 융합 장벽 분포를 도출하기 위한 효과적이고 민감한 도구로서 가우시안 분석 레시피를 검증하였으며, 실험적 불확실성을 처리하는 데 있어 기존의 DD 방식보다 우월함을 증명하였다. 또한, 단일 세트의 $^{140}$Ce 변형 매개변수를 사용하여 세 가지 서로 다른 투사체 시스템($^{16}$O, $^{36}$S, $^{28}$Si)의 데이터를 성공적으로 재현한 것은 채택된 방법론의 신뢰성을 강조한다. 저자들은 이러한 다중 시스템 제약이 물리적으로 정당한 핵 구조 매개변수를 도출하는 데 필수적이며, 가우시안 분석 접근법이 고유한 핵 변형을 탐구하는 강력한 프로브 역할을 한다고 상정한다.